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APPSJDIX 



SIWOPSIS 


Network flow pioblems have received increasing attention in recent 
years due to significant advances in iiiqjlecientation tedanology, coE?)uter 
software and laodelling techniques. In this dissertation, we stu3y the 
role of parametric programming in solving some goieralizations of the 
minimum cost flow problem. It is established that parametric programming 
Can be used to obtain simple and intuitive d-gorithms for several network 
pioblans. In doing so, we preserve some of the essoatial features of 
most network flow problems, namely, the basis triangularity. 

Since our approach is based on parametric programming for network 
flow problene, we first study parametric analysis of some network 
problems to get better insight. We consider the mapciraum flow problKi, 
the minimum spanning tree problea and the shortest path problem and 
develop procedures for their parametric analysis. The analysis is 
based on the linear programming formulation of the network problems. 

We show how this analysis can be used to solve the minimm ratio spanning 
tree problem and the minimum ratio path problem in acyclic networks. We 
also consider capacitated tr^sshipmeit networks in i^diich the suppli^ 
at source nodes, demands at sink nodes and capacities of arcs are linear 
functions of a single parameter, Aa algoritihm b^dd on the parametric 
analysis of the tr^sshipment problem is developed to determine all 
values of the parameter for vMcii a feasible- flow exists. Wiai the 
tr^sshipment network is s~t planar, properties of the dual network are 



used to develop a simpler algoriliim based on the parmetric ^alysis 
of the shortest path problem. 

We next consider the constrained network capacity expansion problan. 

which is essentially a generalization of the constrained mini -mi wn cost 

flow pixiblem and the network edacity expansion problem. In networks, 

where capacities of arcs can be increased by incurring some additxonal 

cost, the constrained netvrork edacity expansion problem is to minifflize 

a linear obj ective fmetion with the additional constraint that the 

cost of flow and capacity expansion does not exceed a prescribed budget 

D. We suggest a parametric approach to solve this prohlaa. We define 

the concept of optimum basis structure for the constrained network 

capacity expansion problem. The optimum basis structure is thoa used 

to parametrize D, The algorithm produces ah optimum solution for all 

values of D greater than the prescribed value, Almost all the 

computations are performed in integers and a near-optimum int^er 

solution is readily available. It is demonstiuted that the algorithm 

can he used to solve a variaht of the constrained network capacily 

expansion problan and a bicriteria network problan. Some special cases 

■of the constrained network edacity problem are also considered and 

restilting simplifications in the algorithm are pointed out, 

■ t 

We also consider capacity expansion of a capacitatai transshipmait 
network where the si^>plies at source nod ^ aad demands at sink nodes 
are linear functions of a parameter. An algorithm is suggested to 
determine the least cost of increasing arc c^aciti^ for all values of 
the parameter for which a feasible flow exists. 



Hie appioach for the constrained netvaork capacity exp^sion problem 
is extended to develop aa algori'Sm for the TninlTnum cost flow problem 
with piecewise linear convex cost fmctions. It is wellknown that 
this convex cost network flow probloa can be transformed to the mint Tmim 
cost flow problem by introducing additional variables, Ihese additiKial 
variables are considered implicitly by defining the concept of optimum 
basis structure for the convex cost network flow problem, Ihe optimum 
basis structure is then used to parametrize v, the flow to be transshipped 
from source to sink, Ihe resulting algorithm irplicitly enumerates 
shortest paths from source to sink and at^ents flow over these paths. 

Ihe algorithm is shorn to be polynomially bounded, 

¥e next show that the proposed algorithm for the convex cost network 
flow problem can also be used to solve the following problous; 

(i) the time-cost tradeoff problem in a OIM netvork vhere time-cost 
tradeoff functions for various arcs are given by piecewise linear con'vex 
ftnctionsj 

(ii) the maximum flow problem in a network wi1h piecewise linear concave 
gain functions} and 

(iii) optiiaum allocation of a giv«i budget to increase capacities of 
various arcs in order to maximize the flow in a network where costs of 
eapacily expansion are glvm by piecewise linear convex ftinctions. 

The algorithm for the convex cost network flow problem has bear 
coded and tested on a nxmher of problons. OoEputational results 
indicate that the algorithm can solve large-sized network problems in 
reasonable amount of time. 



GHAPIES I 


HfTHSDTJCTIOU 

1, 1 INTRODUCTION 

Network flow problems are well-inown for their diverse ^plications 
as well as efficient solution procedure. -Apart from 1iie obvious trans- 
portation and coimrmnication ^plications, a wide variety of problems 
arising in facility location [l3, 23, 76, 77 ] , producticm planning [25, 
32, 91, 95 ] , project manggemeit [40, 58] , machine loading [l5, 29] , 
operations scheduling [82, 87] , cash managemmt [44 , 88 ] can be 
formulated as network flow problene, Recent developments in solution 
methodology, in^ilementation technology and inproved computer software 
have enabled the network flow problems to be handled more effectively 
and conveniently, 

The minimum cost flow probloa is peihfig>s "aie most popular network 
flow problem. Some of the reasons for its populari"^ ares (i) simplicity 
and intuitive appeal of the problem as well as its solution procedures; 
(ii) integer valued optimum solution; (iii) little computer memory 
requiiemait; (iv) efficiency of conputer cedes in solving very large- 
sized problems; and (v) the abili^ to foimilate several complex problCTS 
as a sequence of the TniniTOiim cost flow problaB, 

Eractical considerations such as losses and gains associated 
with flows, additional constraints arising due to br^getary and time 
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constraints, fixed charges associated with new routes, nonlinear costs 
Tfihich are oftoa oacountered in practice etc, have led the researchers 
to investigate gmeralizations of the mninBm cost flow problem. Such 
generalizations have been sought with a view to enhance the scope of 
applicability of the miniinum cost flow problem. It is obvious that 
the structure of the original problem is affected, to a certain extmt, 
by these gmeralizations, Ihe essential effort has been to consider 
only those goieralizations for which it is possible to retain coinputa- 
tional efficiency of the solution procedures. 

In this dissertation, we explore the vise of parametric programming 
to exploit special structures of some generalizations of the minimum cost 
flow probloii. Specifically, we consider the following problems: 

(i) the minimum cost flow problem with an additional linear constraint; 

(ii) the network edacity expansion problems; 

(iii) the minimum cost flow p3x>blem with piecewise linear convex cost 
functions; 

(iv) the time-cost tradeoff problem in CJM networks; and 

(t) the maximum flow problem in networks with piecewise linear 
concave gain functions, 

fhe constrained minimum cost flow problem and the network capacity 
expansion problem have been studied by several researchers. We consider 
these two problens in a unified ncdel, which we term as the constrained 
network capacity expulsion problem, ilhe approach commonly used to solve 
a constrained minimum cost flow problem is to develop the sinplex 
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adaptation for it. If the aspect of capacity expansion is incorporated 
in the constrained mnimim cost flow problem, the siuplecs adaptation 
does not exploit the structure effectively^ We ^ow that if parametric 
pjoDgraHming is used to solve the constrained network capaci'ty expansion 
problem, the structure can be exploited vepy effectively. 03he main 
advantage of this a^iproach is that the basis triangularity is preserved, I 
liirtheimore, a near-optimum integer solution is readily available. 

Some variations and special cases of the constrained network capacity 
expansion problems are considered. We also consider capacity expaasion 
of a transshipment network where supplies and demands are linear function 
of a parameter, M algorithm is outlined to deteraine the least cost 
of increasing arc capacities for all values of the parameter for which 
a fesisible flow exists. 

Ihe parametric approach for the constrained network capacity 
expansion problem is extended to develop an effLciait algorithm for the 
convex cost network flow problem. 3?he algorithm iirplicitly eaumeratos 
shortest paths from source to sink and augments flow over these paths. 

It is observed that this algorithm can be used to obtain the project 
cost curve of a CPU network where time-cost tradeoff functions for 
various ares are piecewise linear convex functions, Ihe algorithm, 
thoi, ImpT icitly Qiumerates a cutset at each iteration and modifies the 
activity times of all activities belonging to this cutset. We also 
^ow that the algorithm for the convex cost network flow problem can 
be used (i) to obtain the maximum flow in a network with piecewise 
linear concave gain flmctionsj and (ii) to optimally allocate a 
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prescribed budget to increase the cs^acities of various arcs to maximize 
the flow in a network vdiere cost of capacity esp^sion is given by 
piecevsise linear convex functions. 

Since our approach is essentially based on parametric progismming, 
a better insight into it is gained by studying the parametric analysis 
of some network probloas. Specifically, we consider the maximum flow 
problem, the mi n imum spanning tree problem and the shortest path problem i 
and develop efficimt solution procedures for their parametric aialysis. 

We show how this analysis can be used to solve some TniniTmiTn ratio 
network problems and a feasibilily problem arising in capacitated 
transshipmmt networks. 

1,2 FBELim^ ABIES 1 

Some notations and w^Ll-known conc^ts of graph theory are used 

!| 

throu^out the thesis. Ibr the sake of completeness, they are given 

I 

below, 

A directed graph G == (ur,A), consists of a finite set F of 
elemoits, called nodes, aid a set A of ordered pair^ of nodes called arcis 

V ■ f 

A directed network is a directed graph in vdiich numerical values are 
attached to the nodes and arcs of the graph. Let n = InI and 
m = IaI • Sie two specified nodes s and t are called the source 
and the sink respectively, 

4n arc (i,3) has two end points, i and and it is said to be 
incident from node i and incidait to node j . Let l(i) and 0(i) 
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dmote, respectively, the sets of arcs incident "to ^hd incident from 
node i, !I3ie dpgree of a node i is the number of arc^ incident to 
or incidoat from that node, 

A Path in G = (N,a) is a sequaice i^yi^, ..«,i^ of distinct 
nodes of If such that either (i^, e A or ® ^ 

each k=1 , , , , , r-1 . A directed pato is defined similarly, except that * 
for each k=1,...,r-1, A cycle is a path together 
with an arc (i , i ) or (i , i ) , A directed cycle is a directed 

r I II* 

path together with the arc (i^, i^), 

A graph G' = (lf», A«) is a sribgraph of G = (1 T,a) if N' £1 i 
and A’ c A. A graph G* = (IT' , A’ ) is a spanning subgraph of 
G = (H,a) if N* =H and A» c A. ; 

Two nodes i and j are said to be connected if idiere is 
atleast one path between them. A graph is said to be connected if 
all pairs of nodes are connected; otherwise it is called disconnected . 

A set Q £A such that the grs^h G' = (ir,A-Q) is disconnected and 
no subset of Q has this property, is called a cocycle of G, A cocycle 
a cutset if it disconnects source and sink. 

A graph is acyclic if it does not contain any <ycle, A tree 
is a connected eicyclic graph. A subtree of a tree T is a subgraph 
of T as well as a tree, A tree T is said to be a spanning tree 
of G if T is a spuming subgrg^ih of G, Arcs belonging to a span n i n g 
tree T are called tree-arcs ^d arcs not belonging to T are called 


nontree-arcs . A spanning tree of G = (K,a) has exactly (n-l) tree-arcs 
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In a spanning tree, there is a unique path between any two nodes, 
iddition of any non-tree are to a spanning tree creates exatfJtly caie 
cycle. If any arc in this cycle is dropped, the resulting graph is 
again a sp s a inin g tree. If any tree-arc is dropped from a spaaning tree, 
two subtrees are formed. Ares v&xxch have their aid points belonging 
to the differeaat subtrees constitute a cocycle. Ibis coqrele is a 
cutset if the dropped arc b^ongs to the unique path in the spanning 
tree from source to sink. If soy arc belonging to this cocycle is 
added to the subtrees, the resulting gr^h is again a spanning tree, 

1,3 OUTLim OF TEE THESIS 

In Chapter II, we study the parametiic ^alysis of some network 
problems. The following problems are considered: 

(i) the maximum flow problem; 

(ii) the minimum spanning tree problem; and 

(iii) the shoirbest path problem. 

In the maximum flow problon, the lower and upper bounds on the 
flow in each airc are considered to be linear functions of a parameter X . 

In the minimum spanning tree gaid short^t path problems, arc lengths 
are parametrized, Sie parametric analysis is based on the linear 
programming formulations of the probleas aad the structure embedded in 
the problems is used to simbli^ the computaticais. We show how this 
analysis can be used to soliiil ri) the minimum ratio spanning tree problcaa; 
jaid (ii) the minimum ratic^||th probloa in acyclic networks. 
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We also consider capacitated transshipment networks in inhic^ the 
supplies at source nodes, demands at sink nodes aid capacities of arcs 
are linear functions of the parameter X. An algorithm based on the 
parametric analysis of the transshipment pixiblem is developed to 
determine all values of X for which a feasible flow exists. A simpler 
algorithm is outlined for s-t planar networks, Yaiious special Cases 
of this problem have been considered by researchers. Our contribution 
is in proposing simpler algorithms for more goieral classes of networks. 

In C2iapter III, we consider a generalization of the constrained 
mi n i m u m cost flow problem and the network capacity expansion problem, 
Tihich we term as Idie Constrained Network Capacity Expansion (CNOE) 
problem, llhe mathematical statement of the GNGE problem is as follows; 


Minimize Z = I c. . x. . , 

( 1 , 3 ) E A 


( 1 . 1 ) 


subject to 


1 X.. 

(j,i)el(i) 


l^-v, if i = s, 

I X. . 0, if i / s,t, y ie N, 


(i,j)eo(i) 


V, if i = t, 


0<.y^<a^^-b.., y (i,j)eA, 

I ®l-l *ii ^ ^i1 ^ii 

(i,j) eA {i,j)cA 


( 1 . 2 ) 

(1.3) 

(1.4) 

(1.5) 


vhere x. . is the flow aod y. . is the capacity expansion of arc (i, 3 ). 
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?ix)blems related to the constrained TniniTnum cost flow pioblOE 
have beoa considered by researchers and their anphasis is on developing 
primal sinplex adaptations [43* 54, 62 j • Shis approach, however, 
does not appear to be promising for the C3JCE problem, She constraints 
(1.3) are not ordinary upper bound restrictions and may not be 
considered inplicitly in the sinplex ad^tation. Hence the need for 
an alternate approach arises. 

We suggest a parametric approach to solve the C3JCE problem. 
Properties of the optimum solution of the CHGE problem and extensions 
of the concurs of bounded variable linear programs are combined to 
evolve the concept of optimum basis structure for the OTCE problem, 

Ihe opti m um basis structure considers the variables y. . and the 
constraints (1,3) and (1.4) iiif>licitly, ther^y reducing the size of 
the optimum basis considerably. lUrthermore, the constraint (l,5) 
is considered in a manner vfaich does not affect the triangularity of 
the basis, 

Ihe optimum basis structure is used to parametrize D.. Initially, 
the algorithm obtains an optimum basis structure for a large value of D, 
Ohe value of D is thm decreased continuously gjad an opti mu m basis 
structure is maintained at each step imtil either the prescribed value 
of D is reached or infeasibility of the problan is indicated. Each 
iteration of the algorithm is similar to the primal siii¥>lex iteration 
of the TniniTmgn cost flow probloa. AlEKjst all the computations are 
performed in integers ^ad a near-optimum integer solution is r*eadily 


available. 
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Yfe show that the OTCE algorithm can be used (i) to solve a vaixation 
of the caJCE problem? and (ii) to obtain efficioat ®>lutions of a bicriteri; 
network problem. Some special cases of the dJGE probloa are also consider 
and the resulting sin?)lifications in the GECE algorithm are pointed out. 

We also consider the capacity expansion of a capacitated transshipmca: 
network where supplies at source nodes and demisids at sink nodes are linear 
functions of the parameter X , lEhe problem con^dei«d is to determine 
the least cost of increasing arc edacities so that the changing demmds 
are met by the cdignging sl:^)plies res\ilting from the changes in X . in 
algorithm, essentially similar to the C3JCE algoiithm, is developed to 
solve this capacity expansion problem for all values of the parameter for 
which a feasible flow exists, ! 

!Qie ideas developed in Chapter III are further refined and extaaded 

in Chapter IV to develop an efficient algorithm for the Convex Cost j 

' i 

Network Plow (CCNE) problem. | 

The CCNP problem is to determine the minimum cost flow in a network | 
when cost of flow over each are is given by a piecewise linear convex 
function. iSie CCNP pixiblem has been studied by several researchers and 
following approaches have bean siiggested to handle such problems: 

(i) the primal-dual approach [ 64, 70, 71 3 } 

(ii) the shortest path approach [ 51 ] ; and 

■' ! 

(iii) the negative cycle approach [l1, 35, 36, 59y 67, 93] • 

It is wUll-known that the CCNP problem caa be transformed to the 
•minimum enst: flow probleo, by introducing one variable for each linear 
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segmaat [72] , We note that in the optinun solution, all the aiditional 
vaiiables be considered ingjlicitly, !Ihis observation allows us to 
define the optiauia basis struotxare for the ecu' problen. Dbe optiiaun 
basis stsnicture is then used to parametrize v, the flow to be transshipped 
from source to sink. Initially, the algorithm obtains ah optimum baisis 

structvire for v = 0, The value of v is thm increased and an 
optimum basis structure is maintained at eveiy step until either the 
desired flow is established or infeasibilily of the problem is indicated. 
The resulting algorithm successively augments the flow on "tiie shortest 
paths from source to sink ^ihich are implicitly enumerated by the algorithm, 
Haace the algorithm essentially belongs to the shortest path approach. 

The corputational corplexity of the algorithm is shown be o(mv) , 

Hu's algorithm [51] for the CCaJP problem is also based on the 
shortest path approach. His algorithm solves a shortest path problem 
at each iteration and augments flow over that path. Our algorithm 
generates these paths without solving shortest path problems and hence 
is more efficient than Hu’ s algorithm. 

We show that the GCSTP algorithm cah be used to obtain the project 
cost curve of a CIM network whoi time-cost tradeoff functions for various 
activities are given by piecewise linear convex functions. 53ie 
algorithm, thm, idaxtities a cutset at each iteration ^d activity times 
of all the activities belonging to this cutset are daanged so that the 
project duration is decreased using minimum additional cost, A similar 
algorithm to obtain the project cost curve is suggested by Haillips and 
Dessotxki [ 78 } , Hie basic differoace betwem the two algorithms is 
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that the cutsets are obtained by applying a cut search procedure at each 
iteration in Shillips.and Dessouliy‘s algoilthn, ^ereeis th^ ax^ 
inplicitly enumerated by otir algorithn. The algorithm of HaiUips and 
Dessouiy ae well as those of Fulkerson [4o] and Kelly [58] are 
primarily suggested for networks with linear tine-cost tradeoff Emotions 
but can handle piecewise linear convex functions by introducing one arc 
for each linear segment. This results in ^ilarging the network 
substantially. Our algorithm does not require any additional arc to be 
introduced. 

We also establish a relationship betweoa the CCKF problem said 
the maximum flow problem with piecewise linear concave gain functions. 
Modifications in the CGNF algorithm are suggested to solve this problem, 

!I3ie problem of optimally allocating a given budget to increase 
the capacities of various arcs to maximize the flow in a network, is 
considered by several researchers for differmt -capacity Kcpansion costs, 
Yfe show that whoa capacity eocpansion costs are givai by piecewise linear 
convex functions, this problem caa be solved by the GOJF algorithm. 

To test the computational performance of the CGN’F algorithm, a 
program was writtaa and tested on a nimiber of network problems, 0o3m>tu 
tational results demonstrate the efficioicy of GGtSTF algorithm in solving 
large-sized problQHS, 

cair-^ter Y cont^s the concluding remarks. 

Since the problems considered in various ch^ters are quite 
different, we survqy the literature diapterwise. 



CHAPTER II 


PAEAMETEIG AMAiYSIS OP SOME lEEWOBK PBDBLSiS 

2.2 INTRODUCTION 

Algoiithms for several optimLzatioa problems solve the shortest 
path problem or the maximum flow problem as their subproblans [ 1 , 28, 
51, 78 j . Goaerally, data of the subproblems gmerated in two 
consecutive iterations is different only for a few arcs. Thus, if 
instead of solving the subproblans afresh, parametric analysis is used 
to talce care of the changes in data, significait con?)utational savings 
can be achieved. 

Parametric programming for network problems was initiated by 
Srinivasaa and Qlhompson [85, 86] , wtio developed an operator theory 
of parametric programming for the transportation problem, !I3ais work 
was subsequently extended to the gaaeralized transportation problem by 
Balachandran and Thompson [ 5, 6,7,8 ] . Some results are also available 
on the sensitivity analysis of network problems [ 28,58, 85, 94 ] . 

In this chapter, we develop efficient methods to perform 
parametric analysis of the following network problems; 

(i) the maximum flow problem; 

(ii) the nnnimxm spanning tree problem; and 

(iii) the shortest path problem. 
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In the mggiTfi u m flow piobleia, the lower and upper bomds on the 
flow in each arc are considered to be linear functions of the parameter x * 
In the miniimaa spanning tree and shortest path problems, arc lengths 
are parametrized. We show how liiis analysis can be xised to solve 
(i) m i n i mu m ratio spaming tree problem? and (ii) the minimum ratio 
path problem for acylic networks. 

We also consider multipleusource, multiple-sink: directed networks 
in which the supply at each sotirce node, the daaand at each sink node 
and the capacity of each arc is a linear function of the parameter X , 

!Qie problem considered is to determine all the values of X for which 
a feasible flow exists, Ihis problem has been considered by several 
researchers for certain restricted classes of networks. Doulliez ^d 
Eao [27] considered single:.sotirce, multiple-sink networks T/daeie only 
the demand at each sink node increased linearly with tine. In [28] , 
they studied the same problem for s-t plaiar networks sard suggested an 
algorithm which obtained the optimum solution by solving a sequoace of 
shortest path problms. Minieka[68] extended the problai of Doulliez 
and Rao [ 27] by considering networks with multipleusourees and multiple- 
sinks. In this chapter, we propose simple and efficient algorithms 
for laore goaeral networks in vdoich the capacities of arcs are also 
considered to be linear functions of the pai^eter X . 

We con^der the linear programming formulations of the netwjrk 
problems to perform the parametric analysis and the structure embedded 
in the po^blens is used to simplify the corputations. Initially, an 
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optimun basis for soae value of X is obtained, Optiraality criteria 
of this basis is then used to detenaine all values of X for vdiich this 
basis continues to remain optimm. Haese values of X constitute an 
interval which is known as the (diaract eristic intearval associated with 
that basis. At the end points of the chairaot eristic intervalf a 
siinplex iteration or a dual simplex iteration is perfoiaaed to obtain an 
alternate optimum basis which may allow further increase or decrease 
in X , Ihe feasibility pioblem cited above is also solved by this 
approach. 


2.2 PARJiMETRIC ANALYSIS OF MAXIMUM Fim PROBLEM 


In this section, we suggest a method to perform the parametric 
analysis of the maximum flow pioblen vitoaa the lower and tipper botinds on 
the flow in each arc are linear functions of the parameter X. 

let four numbers a?., al?'.* b?, and b*. be associated with 

each arc (i,3)eA, let a°. + X . and b?. + Xb*. be the lower 

and upper bounds, respectively, on the flow in each arc (i,3) e A, We 

assume that the network contains the fictitious airc ("tts) vdaich 

= a* = b* 

^:s °t£ 

the maximum flow problem can be stated as 

(2.1) 


a. = ac = =0 and b . = «• , Then, for a fixed value of X » 

ts ts ts ts 


Minimize 
subject to 


‘ts 


I X - I X = 0, ie U, (2.2) 

(3,i)el(i) (i,3)e0(i) 




(2.3) 
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Panwalkar [ 75] has studied the parametric analysis of bounded 
Variable linear programs* !Ehe approach suggested by us for the 
parametric analysis of the maximum flow problem is similar to that of 
Panv/allcar [ 75 ] . 

2*2,1 Initial Optimum Basis Structure 

A basic feasible solution of the maximum flow problem consists 
of the basic variable saad the nonbasic variables vdrLch are at their 
lower or upper bounds. Let B,L and U represent the sets of arcs 
corresponding to the basic variables, and the nonbasic variables at 
their lower and upper bounds respectively. We refer to B as a basis 
^d arcs belonging to B as basic arcs , !Che triplet (B,L,TJ) is 
referred to ais a basis structure ^ 

!I3ae basis of the maximum flow problem is a spanning tree, which 
consists of two subtrees joined by the arc (t,s) [22] . Arcs vdaida 
have their aid points belonging to the differmt sxibtr.ees constitute a 
cutset. In the optimum basic feasible solution, this cutset is the 
miniinum cutset. We assume that a value of X » say X^, is knowi for 
vdaich the maximum flow problem is feasible. Let x. . be the optimiaa 
flow and (B,L,U) be idle optimum basis structure for X= X^, 

2.2.2 Characteristic Interval 

It is easy to see that if X is varied, the dual feasibility of 
the optimum basis structure is not affected, but the prim^ feasibility 
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may be violated due to changing lower and tqjper boundsi Ihusj liie 
characteristic interval of (B|IijTl) consists of all the values of X 
for vshich this basis structiire continues to faiiaia feasiblea 

The characteristic interval is detennined from the follbvjing 
consid erations { 

(i) the flow in each aye (i||)eB does not violate its lower or 
v^per bound} 

(ii) the flow in each arc (i,j)el/ is kept at its changing lower 
bound} and 

(iii) the flow in each arc (i»3) e U is kept at its changing upper 
bound* 


It follows f]x>m these considerations that as long as the current 
basis structure remains feasible, flow in each arc (ijj) eB is of 
tho form x. . + (^-X ) z. The following theorem derives m eaepression 
for z. .* 

Theorem 2*1 1 let T. and T. be the resulting subtrees, containing 
' 10 

node i and node ^ respectively, whm any arc (i,3)eB is dropped 

from 1016 basis* let ^1. . and Q. . be the sets of arcs from T 

10 -T-O i 

to T. and T. to T. respectively* Thai 
0 0 1 


T b* T b* + T 

^ kr -v Tv kfc 


(k,il) nu {k,£) eQ.^ nu (k,£) ni* 




I 

(fc,£)e4.ni. 




(2.4) 


z. . 
10 



17 


Proof ; Addition of any arc (k, Jt) e L Uu to the basis creates exactly 
one cycle consisting of basic arcs. Let W, . and W be the sets of 

JC Af "jC Jo 

arcs in this cycle along and opposite to the orientation of arc (k, £,), 
The flow in each arc (k, 5.) e F increases at the rate of b^j_ as X is 
increased. Keeping in view the flow conservation at all the nodes, 
it amounts to increasing flow in the cycle at the rate of If 

the arc ( 1 , 3 ) ^ then flow over it increases at the 3rate of ; 

and if (i,j) then flow over it decreases at the ra.te of * 

We further note that iff (k, «-) ^d (i,3) e 

iff (k, s,) e 5. Hius the first two terms in (2.4) represoat the 
effect of changing tipper bounds of all the ao?cs in U on the flow in 
arc (ijj). Similarly, the last two terms in (2,4) represent the 
effect of changing lower bounds of all the arcs in L on the flow in 
arc (i, 3 ). 

An efficient method to determine z for each (ijj) eB is 
stated below. It is easy to see that the method terminates in o(m) 
iterations, 

Stq> 1 ; Set TT^ = 0, ¥ i e K, Perform the following operations for 
each (k, £) eLU U and go to St^ 2 , 

'k ' ’k "n = 'It - =kji > “ 

'k * 'k ' ‘’kl '1 “ 'll - '’kt ' (l£,lt)sD. 

Step 2: In tile basis select a node of degree 1 ^ bbs ' if ^d find tiie 
basic arc (PfC[) snch that p = i or = i. If 
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p = i, then, set z = - ir. and laadiiy it to it + ir , 

» pq 1 q P q 

If q = i, thaa set z = ir. and loodi:]^ ir. to ir "*■ tTq* 

pq ^ p p q 

Delete the arc (p,q) from liie basis and repeat this step 
■until no arc ranains in the basis, 

Ihe justification of this method is as follows; 

The flow in each arc (k, 2.) eU is kept at its changing upper 
bound, vhich increases at Ihe rate of b*^. To do so, additional 
flow is required at node k and getting available at node 2 at the 
rate of b*^. Similarly, on account of aay arc (k,Jl) eL, aiditional 
flow is required at node k and getting a'vailable at node £ at the 
rate of tt^, at the end of Step 1, denotes the cumulative 

requirement of flow rate at node i. The tree structure of the basis 
is thm lised to determine unique values of z. . which satisfy the 
flow reqToirement of all nodes. 


Having determined z. ., the characteristic interval is obtained 
by using the following inequalities; 




(2.5) 


which csn be restated as 


a..+ (X-Xj (X-X„) (X_XJ b*,, ¥(i,j)eB, 

( 2 . 6 ) 


®ij^ ®ij -^ij '"'"O' ij-aj ' o' ij' 


vtiere a. . = a?. + X a*, aid b = b° + X b* . 
ij ij o xj ij 13 o X3 


I 

! 
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let us drflue the uui*ers and ^ 

for each ( 1 , 3 ) eB as follows; ^ 


AX . = < 
ij ^ 


_ f ‘^13- « ^d3 < . 

A1 * - 

(_ “ , otheiwise, 

. 'If 

^ » , Otherwise, 


AX^-. . = i 


li ii " ' ~i -i ““t -5 » •^■'- Z . . > b* . 

«' -^J ij 13 i3 ±3 ^ 


(2.7) 


(2.8) 


z. - > a*. 


AX^ 


ij 


otherwise, 


(2.9) 


AX^ 


■ij 


I <*’13-^13 < b* , 

^ l-». 


otherwise. 


(2,10) 


let 


Using (2,7)-(2,10), (2,6) reduces to 

^ (i,a)eB. 


( 2 . 11 ) 


and 


Xg = X^ + non. { 


(i» 3 )e B 


-r^l —2 , 

AX 


ij ■ 13 


(2.12) 




(2.13) 
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The conditions in (2.1 1) now ffim pi 1 -fy to 






(2.14) 


Bans ( Xg) is the characteristic interval of (S,Jj,U), 

Bor all values of X e (l^, Xg), (B,L,II) is tiie optiiaum to^Ls structure 
and the optimum flow x. . is givei by 


, ¥(i,j) gB, 

S.. + (X-X^) aj^ , ^i±,2) £ L, (2.15) 

^ hi ^ ‘’13 > 


If it is required to solve the maximum flow problem for X > 

( ^ “^ ^Lg) t thai a dual simplex iteration is performed to obtain 
alternate optimum basis structure for X= X-{\=X^), 

2,2,5 Dual Sirpleg Iteration 

let X. . denote the optimum flow for X=X^* AtX=X_ 

3-3 a n 

— -t - ■ '2 

the flow in an arc (p,q) e B, for whidi X^ = X^ + min. { AX AX 

equals its lower or xapper bound. If X = X +‘ax^ , thm 

B o pq' 

a^. and x = a?. + X^ a?.. However, if X„ = X + AX^ , 

3-3 3.3 13 13 B ^3 » B o pq' 

then z. . > b*. and x. . = b?. + X*_ b*. , If X ig increased 

13 3.3 13 13 B 13 

further, without chaiging the basis structure, ttie condition (2.5) is 
violated. Thus a dual siiiqjl ex iteration ±f, performed to obtain an 
alternate optimum basis structure fbr X = X^ viiich may allow further 


increase in the value of X , 
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iEhe dual siuplex iteration, is performed by dropping the arc (p, q.) 
from the basis «d selecting a nonbasic arc to mter the basis. When 
the arc (p,q) is dropped from the basis, tvso sidbtrees are formed, let 
f and f be the resulting subtrees containing node p said rujde q 

P St 

respectively, ill the arcs vdrich have their one aid point in 2 and 

p 

another in 2^ constitute a coqycle Define the orientation of 

the coqrcle along (p,q) if (p,q) leaves the basis at its upper 

bound, and opposite to (p,q) if it leaves the basis at its lower bound, 
let Q and Q be the sets of arcs in the cocycle along ®d opposite 

irH. irTL 

to its orientation respectively. 

Choice of the entering arc depoads upon tbe entries in the 
associated aurplex tableau in the row corresponding to ^d the 

relative cost coefficioats of the nonbasic variables. We now state 
two linown results, which may be fomd in [l0, 22, 2?] . 


Property 2,1 : A nonbasic variable has an mtiy ~1 in the row 

corresponding to if aid only if (lc,S.) e where 

E = (Q ni) u (Q nu). 

pq Vpq r ^q 

Property 2,2 t let ioiote the relative cost coefficioat of 

tbe variable 2hm 




0 , ¥ (fc,£) 

-1 , ¥ (k, A) e , 


(2.16) 


where 




is the set of arcs in the mininaaii cutset. 
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Let us define the sets E aod E as follovra: 

pq -pq 

Ep, - e <!jg,)= 0} , (2.17) 

s : (aj^_ ej^u) = _1 } . (2.18) 

If E is eii53ty» then all the mtiles in the row ooirespoading 

ir M 

to are nonnegative. Ihis indicates that the dual sin^lex iteration 

cm not be performed md the maximum flow problem is infeasible for 
1 > Xg» If E is not Qi?)igr, thm my arc belonging to B cm 
enter the basis. It is easy to see that the basis chmges by the 
d\ial simplex iteration but the minimum cutset remains unchanged. 

However, if E is enmty, then my arc belonging to E cm alter 
the basis. In this case, miniTmitn cutset also chmges, Bie same 
procedure cm be r^eated to determine optimum soluHons of the maximum 
flow problaa as X increases further. 

To perform the parametric malysis for X < ^ arrc (p,q)cB 

satisfying X.^, = X + max. {dX^ } leaves the basis, 5toe 

" — B 0 — pq^ — pq 

arc (p,q) leaves the basis at its lower bound if X^ = X^ + A X^ ^^ 

2 

and at its upper bound if X.^ = X + AX , The entering are is 

— 3 o pq. 

selected in the same manner as it was selected for X = X^, 

2,2,4 Numerical Example 

We now solve a numerical exaB¥>le to illustrate various steps 
of the parametric malysis of the maximraa flow problm. !i3ie netirork 
is shown in ilg. 2,1, Eor the sake of siD 9 )lici 1 ?y, the lower bound 
on the flow in each arc is eissumed to be zero. Nodes 1 md 6 are the 
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souice acd. sink respectiTelyi Hie paraffifiiiiic gaiaiysis is perfoimed. 
for all feasible values of X > 0| 


Hie st^s of the parametild maiysis are tabtilated in Table 2<1| 
In the tablej the symbol + indidates the basic ard leaving the basis 
said the symbol + indicates "ttie nonbasic arc entering the basis# Since 
the parametric analysis is performed for incr^sing values of X , the 
values of ^nd required* In the first two 

iterations| basis as well as the minimum dutsets diange* In Hie Hiird 
iteration, only the basis changes and the mlDiiniBn cutset roBains 
unchanged* Hie basis obtained in the fourth iteration is optiimm for 
all values of X ^ 4 . Hie basis in various iterations eure shown in 
iig. 2,2. 


2*2,5 A Special Oasei ELanar Uetwirks 

A network is planar if all the arcs connecting the nodes om be 
drawn in the plane such that no two arcs cross each other, A network 
is s-t pl^ar if by adding ah arc from node t to node s, the 
resulting network raoains planar* Associated with each s-t plmar 
network is a laiique dual network, which csn be constructed by the meHwd 
described in f26 ] . We denote the dual of the pl^ar network 

ff = (K,A) by ff* = (#, A*). 

In Gr* = (N*, a*), there is a source node as well as a sink node, 
let s* aid t*' daiote the a^urce and sink in &* respectively. Hirtherj 
there is one to one correspondence betwem ares in A and A*. It is 
wellloaown [37 ]that when the lower boiaad on the flow in each are is zero| 




25 


!I?able 2,1 : Solution of the networfc flow piDblem 


SS"" i-B \ 

Hunijer 


Basic Atos 


Nonbasic Arcs 


1. 


0 1 16+6X 


(i»0) 

X. . 

z. . 
IJ 


(2,4) (3,5) (4,6) (5,6) (6,l) 
6 10 16 


10 

2 


1 


2 6 


CO CO 


CO CD 


(1,2) (1,3) (3,2) (4,5) 
6 10 0 0 


1 2 22+4X 


(i,d) 


x. . 
ij 

2 . . 
13 


“ij 

AX^. . 
IJ 


(1,2) (3,5) (4,6) 
10 12 10 

2 2 2 

00 CO CO 

00 -j oo 


(5,6) (6,1) 
12 22 
2 4 


(2,4) (1,3) (3,2) (4,5) 
10 12 0 0 


(i,3) 


2 4 26 +3 X 


ij 


ij 




(1,2) (3,2) (4,6) 
12 0 12 

1 1 2 


2 

4- 


(5,6) (6,1) 
14 26 

1 3 


(2,4) (1,3) (3,5) (4,5) 
12 14 14 0 


4. 


32+3X 


(i,3) 


X. . 

13 


13 

_1 
4^ ij 


(1,2) (1,3) (4,6) 
14 18 16 
2 12 


AX 


13 


(5,6) (6,1) 
16 32 
1 3 

00 CO 

OD CO 


(2,4) (3,2) (3,5) (4,5) 
16 2 16 0 












Basis in various iterations 
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the minimum cutset in & coiresponds to the ^Kjrtest pat^ from s*' 
to t* in G- , In such cases, parametric aiaalysis of the maximum flow 
problem reduces to the parametric analysis of the shortest path probloa, 
for vdiich a simpler aid more efficiait algorithm is d^ciibed in 
Section 2,4* 

2.3 FARAMETBIC MAMSIS OF MINIMUM SPANNING TREE PROBLEM 


In this section, we consider the parametric analysis of 1die 
minimum spanning tree (MSI) problsa. Associated with eacii airc (i,j) eA 
are two numbers e?. and c*. and the loogth of the arc (i,3) is 
c?. + Xc* . . The direction of arcs is ignored in this section. 

let T be the minimum spanning ' tree for some value of x • 
let the set P, consist of arcs in the unique patia from node k to 
node £, The necessary ^d sufficient conditions for I to be a 
MSI [38 ] ares 




(2.19) 


where I = A-T, 


Define the nunhers snd 2L]j£ each (k, Ji) e I as follows: 


k£ 


= ^ 


min. ^ (of — c?.)/(c*.— c* ) } , 
/. .N Z kZ ±y ^ ij kr * 

oa , if is empty, 


( 2 . 20 ) 


28 


r Uc® -c° )/(c| -c| )} , 

I (ijj) ^ ^ ^ ^ 

\:il " I 

- 00 , if is en^ily, 


( 2 . 21 ) 


*er3 = {(i,j) e Pj^^i o*^ > e» ) and 4^- {(1,3) e P^.^: 
Using (2.20) ^ (2.21), (2.19) reduces to 


. <<#-}. 

13 kSL. 


(2.22) 


#iieh further simplifies to 

At i. ^ 1 f (2.23) 

lEhus ( ^T» characteristic interval for T, Star 

all values of X e (At» ^ is a MST aad its loigth is 

( I c° + I c|), 

(i,3)el (i,j)el 

To perfoim the parametiic analysis for X > X^, identify an 

arc (u.v) e T for which X = X aud an arc (p.q) e 1 for 

vdiich (c° - c° )/(c* - c* ) = X , ^ alternate MST for X= X_ 

pq UTT ' uv pq uv T 

is T* = T U { (u,v)} - {(p,q)} • In the similar laanner, "ttie paiumetric 
snalysis for T* can be performed, 

r 

SimilarLy, to perform the parametric analysis for X < ^ , • 

A. 

identify an arc (u,v) ^ T for wiadh = X.^ and an arc 

(pjQ.) ^ ® for which (c° -c° )/(e* -<# ) = X , . 

^ ' — uv P9 nTT uv pq ~uv 

MST for X = ^ is T» = T U {(u,v) } - {(p,q)}. 


M alternate 



Hie diaractenstic intervals of MST*S obtained in the course 


of parametric analysis do not overlap. It is saLso ea^ to see that 
once X leaves the characteristic interval of a MSI, that ti^ee never 
reappears, Hie total number of MSl’s is clearly finite. Moreover, 
an upper bound on the nuniier of MSf’s can be stated as follows: 

Hieorem 2.2 : jto upper bourui on the number of MSf’s obtained, as X 

varies in the interval -« < x <»» is m(ia-2)/4 if m is evaa, 

2 

and (m-l) /4 if m is odd. 

Proof . We refer to c*. as the incremental laigth of arc (i,3). 

Let us number the arcs ais a^,a2,...,ajj^ in the increasing order of 
their incremental laagths aid consider the parametric analysis for 
increasing values of X . It follows from (2.20) that udaenever the MSI 
changes, the oateiing arc has smaller incremental length than the 
leaving airc. Consider an arc a^. SChere are atmost (i-l) arcs which 
have incremmtal lengths less than that of a^^ and atmost (m-i) arcs 
which have incremoital lengths greater than that of If a^j^ is a 

treeu-arc in a MSI, then it may become a nonti«e-arc by being replaced 
by a nontree-arc of smaller incremental laigth vdiich can not ha^spen 
more thaa (i-l) times. Similarly, if is a nontree-airc, thaa 

it become a tree.,arc in a MSI by replacing m. arc of larger incremental 
loagth and this would not happoa more than (m-i) times. Thus 
may cause MSI to chgnge atmost min. {i—1, m-i} times. When this value 
is summed for all the arcs, we get iiie desired expi^sxons. 
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2,4 PARAMBTEIC ANALISIS OF SEORHEST PAHS PRCBIEM 

In this section, we consider the parajnetric analysis of the 
shortest paith problem. Associated with each arc (i#3) ^ A are 
nuinbers c° and c* and the length of the are (1,3) considered 
to be c + X e* , It is assumed -that the network does iK>t ccmtain 

XJ XJ 

negative cycles for ary value of liie parameter x in the laoge of the 
analysis. 

Ihe shortest path pixjblQn, for a fixed value of X , can be 
stated as 

Mi nimi ae I (c?, + \<^.) i. . , (2.24) 

(i,3)eA 43 43 A2 ’ 


subject to 


1 X. . - ^ X. . =■<( 

(j»i)el{i) (i,j)eO(i) 


-(n-l), if i=s, 
,¥ ie F, 
1, otherwise, 


(2.25) 


^ (ijj) ^ A. (2.26) 

It is wellbnom [22] that ay optimtm basis of liie dbortest 
path problem is a spanning tree. In this speaining tree, there is a 
unique directed palii from source to eveiy other node and liiis pato is 
the shortest path from source to Idaat node. Bier^ore, this spanning 
tree is toown as the ti*ee of ^rtest pabhs rooted at source. Let iP 
be the tree of shortest paths for some value of X . Ihe necessary and 
sufficient conditions for T to be a tree of shortest paths are: 

(-r° cJ, 


(2.27) 
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= A^f, and 

IT® 

d 

and IT* 

are 

0 

II 

0 m 

and 

.r° 

tr . — IT . = 

3 1 

0 
c. . 

TT* =0 

s 

^d 

II 

yn 

1 

c*. 

2.3 


( 2 . 28 ) 

(2.29) 


— kJl numbers defined, for eacii arc (k,fi,) e i, 

as follows; 


- ’k - ("* -”k - “kP ' “■ 

CO ^ otherwise, 


(Tr° _ TT? + c° )/(n* - TT* - cto)»i^ « -< 


^k. M 

Li “ , otherwise. 

Using (2.30) and (2.31), (2.2?) reduces to 

•( ^^(k,t)eT, 

which further simplifies to 


Z k 


“E; <0. 


(2.30) 


(2.31) 


( 2 . 32 ) 


h - ^ » 


(2.33) 


wiaere = min, and ^ = max. * {X,J. 

^ (k,£)ea} ^ (k,£)ef 

^Qius (^, T^) is the characteristic intesrval of f, for all 

values of X e (^, Xj), 5 is a tree of shortest paths and length of the 

shortest path from source to any node 3 e N is + X ir*- , 

3 3 

!Bo perfoim the parametric analysis for X > , a simplex 

iteration is performed to obtain an alternate tree of shortest paths 
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for X = X let (u,v) e T be an arc for isfcich X^^ = X^, and 
(PjQ.) ^ ^ be liie unique arc for which q=v'» Bie tree f' = f U{(ujv) } 

- {(p»q.)} ia an altemate tree of shortest paths for X= X^,. Reseating 
the same procedure, subsequent tree of shortest paths can be obtained 
as X increases further. 

Similarly, to perform the parametric analysis for X < select 

an arc (u,t) e T for #iich ^ " iim unique arc (p,q) 

'IIV* «*» 

such that q = v, fhe tree I* = 1 U {(u,t)} — {(pjQ,)} is an 
alternate tree of shortest paths for X = 

2. 5 PARAMETRIC NETWORK FEASIBILITY PROBIEM 

In this section, we consider multiple-source, multiple-sink 
directed networks in vshich supplies at source nodes, demaids at sink 
nodes and capacities of arcs are linear functions of a paraaeter X • 

An algoiiliim, based on the parametric analysis of liie transshipment 
problem, is dereloped to determine all the values of X for a 

feasible flow exists. A simpler and more efficieat algorithm is 
suggested for s-t planar networks. 

Let SczR denote the set of all source nodes and I c: R denote 

the set of all sink nodes, where S fl 1 ^ , Let the supply of source 

node i e S, the demaad at sink node i e I, capacity of arc ( 1 , 3 ) e A 

be -.r° -Xi^, r° +Xi^ snd b°^ respectively. Here r?, 

b? and b*. are all fixed and X is a parameter ^^ch is permitted 
ia ^ 

to vary . If for eadi i e H — S U f we put r? = ^ then without 

loss of generality we can assume •Kiat ^ == I r* = 0 . 

ie N ieH 




We 3ref er to the problem (2.34)-(2.36) as the Parametric 
HetTOrk feasibility (MP) problem. Although, it is assimed in ■&© 
problem that the lower bomd on the flow in each arc is zero, it 
is by no means restrictive. Problems with nonzero or paxajaetric lower 
bound on the flow in each arc cgai be converted to the PSP problem by a 
simple transformation of variables. It should further be noted that 
the problem to determine ^ can be solved by rolving the HfP 
problem in vfaich b* is replaced by -b* for each i s H and c^^ 
is replaced by -c* for each (i,3) ^ A. 

In principle, the MP problem can be solved as a linear piograa 
by considering X as one of the variaibles. But in doing so, "tiie 
structure of the constraint matrix is signiflc^tly affected. We, 
therefore, consider X as a parameter and ret^n the structure of the 
problem. We then use "tiie concepts of parametric linear prograBHoing 
to determine the range of feeisible values of X. 
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2.5.1 Developmgat of the Algorithm 

!Sie MP problem is a bounded vaii^le linear program. In a 
basic feasible solution of the INP problem, let B,Ii and U r^r^ait 
the sets of asxjs corresponding to the basic variables and the nonb^ic 
variables vdiLch are at their lower and upper bounds respectively. We 

refer to B as the basis and -fee triplet (B,I,U) as the basis 
structure . It is kiom that the b^is of the I?IP problem is a spanning 
tree. 

We assime that a value of X , say is laaom for which the 

BTP problaa is feasible. A basic feasible solution for can 

be obtained by the method described in [ 10 j . In this solution, 
let X. . represent the flow on arc (i,a) e A and (B,1,U) be the 
associated basis structure. If X is increased, the feasibility of 
this basis structure may be violated because of changing supply-demand 
requiiremoits and the changing arc edacities. We first determine the^ 
largest value of X upto whidi the current basis structure continues to 
remain feasible. Ihis largest value of X , denoted by x , is 
deteimined from the following considerations s 

(i) !I3ae flow in each arc (i|3) e B remains witiiin its lower end 
upper bounds} 

(ii) the flow in eanh arc ' (i,o) e L is k^t at zero value} and 

(iii) the flow in each arc (i,3) e U is k^t at its dianging r^per 
bound. 
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It is easy to see that as long as the current basis structure 
remains feasible, the flow in each are (ijj) ^ B is of the foim 

+ (1-^) 2 _ . !13ae following theorem derives liie egressions for z_, 

Oiaeorem 2,5 : Let I. and T. be the resulting subtree, containing 

^ tJ 

node i and node 3 respectively, whaa ar^ arc (i,3) e B is dropped 

from the basis. Let and be the sets of arcs from to 

T. and T. to T. respectively, Ihen 
i3 3 ^ 


! = y 1* 

•i ^ V 


b* - 

FT 


(2.37) 


(k, £) n U (k,£) e HU 


Proof : The first term in (2.37), i.e., ^ i* , represoats the 


k e T. 


effect of changing si?)ply-denieaid requireomts^on the flow in the 
aro (i,3). It is the rate at which the net demand of nodes in T. 
increases with X , Since ^ ^ 1 = 0, it is also the 


k el. 
i 


k eT. 
3 


rate at which the net supply of nodes in increases with X. 

Since (i,3) is the only basic arc connecting I. to I., idae flow 

over it must daange at "Hae rate of ^ i* in order to satis:^ Hie 

kei. 

3 

changing demaads by the changing supplies. 

The last tro teams in (2,37) indicate the effect of changing 
arc Capacities on the flow in aro (i,3) . Proof of Haeorem 2.1 
includes the justification for these terms. 


Hhe following is the procedure to determine z^ . The justification 
of this procedtare is similar to the one used for confuting z^^ in the 


maximum flow problem, 
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Step 1 , 


Step 2 . 


Set Perfoim the following operations 

for each (k,2.) e U and go to Step 2: 

In the basis select a node of degree 1, s^ i, aad find the 

unique basic arc (p,q) such that p=d or q=i. If p=i, 

thoi set z = -ir. and modi:^ ir to tr + ir . If q=i, 

pq 1 q P q 

thm set z = tr. and nodify ir to tr + ir . D^ete 

pq 1 P P q 

the arc (p,q) from the basis aid repeat this step until no 

arc remains in the basis. 


Having determined z. x is determined from the following 
inequalities; 

0 < + (x-x) - (X-X^) b|,, ? (i,j) e B, (2.38) 

where b. . = b?. + x 6*. 

10 13 0 13 

Define the nunibers m AX . AX'...; and ax . ior eada (i,3) e B 

as follows; 


AX^ 


-X. ./z. . , if z. . < 0, 

13 13 3.3 


13 


AX 


OO f otherwise, 

, .-X. .)/(z. .-b* . 

p 13 13^ 13 13 


iTb. .-X. .)/(z. ._b*.) , if z > b* , 
ii x:\.^ 13 13' 


3-3 




CO I otherwise, 


i-x. ./z. . , if z >0, 

13 13 3.3 


L 03 , othervdse. 


( 2 . 39 ) 


(2.40) 


(2,41) 
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AX 




(b. -X .)/(z, if z. . <b*. 

^ 13 13^ ' 13 X3^’ 13 13 


«>, otiieiwisei 


Using (2.59)-(2.42), (2,38) reduces to 


(2.42) 


It is ob-vious from (2.43) that 

X = 1 + niio. {AxV, ) AX^^^) • (2.44) 

° (i,j)=B ^ 

A feasible flow x. for X = X is as follows: 

XJ 

^i/ "'ij » 

= < 0 ¥ (i,3) e I , .(2.45) 

[fir *>« • '' A-3) = ” • 

At X = x» flow in an arc (p,q) for which X^+ min. 

*1 ““2 — 

{AX , AX }= X reaches either zero value or its upper bomd, 

PI PI 

It follows from (2.39) and (2.40) that if flow reaches zero value 

then z 0: aid if flow reaches its upper bomd, thai z > b* . 
pq ’ pq pq 

If X is increased further without changing the basis structure, flow 
in arc (p,q) violates its bomds. 0!hxis to find a feasible flow for 
X > X , a dual siii5)lex iteration is perfoimed to obtain ai alternate 
basis structure for X= X, 

Ihe dual simp lex iteration is perfomed by dropping the arc 
(p,q) from the basis and sheeting a nohbasie arc 1» enter "fee basis. 
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Dropping the airc (p,q) yields the cocycle Q t Define the orleatation 

of the cocycle along (p,q) if (p,q) leaves at; its upper bomd 

opposite to (p,q) if (p,q) leaves at zero value, let and 

be the sets of axes in Q along and opposite to ixs oiioatation 

pq 

respectively. 


Ihe structure of the constraint matrices for the MF problm 
and the masimum flow problem are same and, therefore, Property 2,1 
remains valid for the BJP problem. The relative cost eoeffici^ts for 
all nonbasic variables are obviously zero. Hence, if the set 
®pq ' '^4 n I) u n n) is nonemply, then any aio in it can be 

selected to enter the basis, Bie same procedure can be r^eated 
for the new basis struetTore. However, if E is en^ty, thoi the dual 
conplex iteration can not be performed indicating that X = 

It may be noted that the choice of the arc in E , selected to 

enter the basis, is crucial for the number of iterations reqirLred to 

reach X , Several selection rules were tried computationally, 
max 

finally, the following selection rule which reduced the number of 
iterations substantially was selected. It is discussed for two 
differmt cases. 

Case I : The arc (p,q) leaves the basis at zero value. If any 

arc (k,2.) e alters the basis, then z^j^ in the 

alternate basis structure is givai by 
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2 




~z 


P'1 

1.. 


, if 


z + 0 , if 

pq k £ » 


n !<» 

(kjfi.) e n U 


(2.46) 


let denote the largest increase in the value of X , 

that the ajrc (k,2,) can tolerate without going out of the basis, 
ass um i ng that all the other basic ares continue to remain basic. !0he 
expression for AX’^^^ is as follows: 


AX 


kjt 


kji '' pq ki'* pq k£ * 
\ “ » 


(2.47) 


otherwise. 


where x b*, . 


Case II : !Ehe arc (p,q) leaves the basis at its t;q>per bound. In 
this Case, the exp3?essions for and are as 

follows : 


kJl 


“ 0=.'^) b Op, n 1 , 
‘-Sp, n" > 


(2.48) 


AX' 


kl 




b /(z _b* -bf" ) , if 2 > b* + bf , 

kji^^ pq pq kji^ ^ pq pq k£ » 

CO I otherwise. 


(2.49) 


Ihe number AX' can be takoa as a measure of increase in X 
k£ 

permitted by the arc (k, j) without becoming nonbasic. Hence, if 
at every iteration an arc in hating the largest value of 

is selected to enter the basis, larger increase in X per iteration is 
Thi a conclusion waS verified computationally. 
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2,5.2 Description of the Algorit^im 

A stepwise description of the algorithm is given helow. 

Step 1 . Obtain a feasible basis st3?ucture (b,L,U) for X=X^ by the 
method described in [10 ]. let x. . be the flow on arc 

13 

(i,3) £ A, Set J. = \ , Let b. , = b?. + 1 b^. for each 

0 13 13 i3 

(i,3) e A. 60 to St^ 2, 

Step 2 , Compute z for each (1,3) e B. Let .AX^ . . aad &X^. . 
be the numbers defined for each (1,3) e B as follows: 


-1 

AX 


-oc. ./z. . , if z. . < 0, 

13 13 ’ 3.3 ’ 


i3 


< 


00 y otherwise, 


-2 

id 


(b -X. .)/(z. .-b*.) , if z . . > b*. , 
13 i3^^^ 13 13^ » i3 i3 * 


00 j otherwise. 


Let 

E \ 


-1 _2 

“in. { AX . • » AX } 


Update X, b.. and x. . as follows; 
X = X + ^ , 


^13 “ ‘i 3 *‘^ "id ’ 


x.^+ AX z^^ , ¥ (i,3) , 


3-3 


<0 


b . . 

xj 


, f (i,3) e L , 
, ¥ (i,3) ^ ^ . 



If X - CO , go to Step 4; oldaeiwise identi:^ m. arc (p,q) e B for 


min, {A X 


1 


AX' 


}= AX gjad go to St^ 3. 


p<i ' pq' 

step ? , Dropping the arc (p,q) from the basis yields ^e coorcle 

Define the ori^tation of Q along (p,q) if x. . = b 

pQ. IJ xj 

and opposite to (p,q) if x^^ =0. Det and 

be the sets of arcs in along and opposite to its 

orimtation respectively. let H D) U (Op,^ fl U). 

If Ep^ is go to Step 4; otiierwise define a number 

A X* , . for e£U3h (k.il) e E as follows: 

K Ji ^ ^ ' pq 


h, /(-2 -b*.), if X =0 and — z > h* , 

kr pq pq pq ki » 

AX‘ , „ = ^ «) > if ^ ^ ^ 2 > b* + b* 

k£ ! kJc ' pq pq icy* pq pq pq pq k£ 

» t otherwise. 


Identify m arc (u,v) e E for which AX* = max, 

V » / pq u-y V „ 

_ , (^>^)eSpq 

{ AX', } . Bie arc (u,t) enters the basis, 

•k£ 

Update "toe basis structvire and go to Step 2. 


Step 4: Set X = X aid SIOP. 

— max 

2,5,3 A Special Case ; ELaaar Networks 


A singleusource, single-sink network can al^^s be associated 
with eveiy multiplex-source, multiple-sink network in the manner 
descilbed as follows: The network G = (U,a) is augmented by a 

super-source node i end a st 4 )er— sink node t . An axe (s,i) is 

introduced for each i e S wi'^ b° = — r? aid b* = — i^. Similarly, 

si ^ 0 o 

an arc (i,t) is introduced for eadi i e T with b _ = r. 

it _ ^ 

b* = zf, let the augmoited network be denoted by § = (K,a) . 
it ^ 
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Let us assume that fr is s-t pl^ar. It can be easily seen 
that G is s-t planar if and only if G is s-t planar for eTreiy 
coEibination of sourca-sink pair, The dual network of i = (N, 1) cm 
be obtained by the method d^cribed in [ 2 ^] . Let Q* = (ir*,A*) 
denote the dual network. Let s*’ and t* be the source and sink 
in G . Since there is one to one correspond oace betweai arcs in G 
and G*, let (i*,j*) denote the arc in corresponding to the arc 

(i, 3 ) in G. Let M = {(i,t) e I: i e T } and e A*:(i,3) e M}. 


It is easy to see that the HTP problem is feausible for some X if 
and only if flow in arcs belonging to M equals their respective vqjper 
bound. This ia 5 )lies that the HTP problem is feasible for aame X if 
and only if M is the ndnimum cutset with arc capacities as b?.+ lb*.. 

Since the m i n i mu m cutset in 5 corresponds to the shortest path from 

' *" • 

s* to t* in G*, it follows that the IKP problem is feasible for some A 

if and only if is the shortest path from s* to t* in G* with 

arc laigths as b?. + X b*., 

11 11 

Hie ETP algoilthm for s-t planar networks is now obvious. Hie 

parametric analysis of the shortest path from s^ to t*' in G^ is 

performed for increasing values of X wi-fti arc lengths as b?. + X b^.. 

The highest value of X , upto vdii<di continue to remain the shortest 

path from s* to t*", is -file vaiue of X 

. max 

2,6 MINIMUM RATIO NETWOSK PWBmiS 


Parametric programming is used extensive to solve fractional 
programming probloois £24 j 55 2 • Pbllowing the same approach^ we show 
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that liie parametric gjialysis of netirork probloaa be vtaed to solve 
(i) the nanimum ratio spanning tree probloai and (ii) the miniTOBn ratio 
shortest path problem for acyclic networks, 

2,6,1 Minim-imi Batio Spanning free axiblem 


let two nimibers, c. . ghd d. be associated with each arc 
(ijj) e A, !Ihe MiniiBum Ratio Spainlng IPree (MBST) problem is to 
determine a tree 0* for which the ratio 


I 

(i,3)ea?^ 


c. . 


/ 


I 

(1,3) 


id 


(2.50) 


is minimm. It is assimed that for every spanning tree I, 

I d ^ 0. 

(1,3) e T 

Chandrasekaran [l4] has provided the following characterization 
of the MEST: 

niieorem 2,4 ( LhI ) ; A spanning tree T* is a MEST with m i nimum 
ratio as k?^ if and only if iP*’ is a MST with "tiie length of eeich 
arc (i,j) e A as 

Using this characterization, the following method determines 

a MESS: 

let k*^ be an arbitraiy number. Determine the MSI wilh the 
length of each arc (i,3) ^ A as (c^.-.k°d^^). let be the lEI 

thus obtained. let M° = I Jc -k°d ). If if = 0, then 1° 

(i,3)e3:° ^ 

is a MBSI. If M° < 0, then treating -tiie lengtix of each arc (i,3) 
as (c. .-k d. ), perform the pareffaetric analy^s of the MSI problem for 

13 13 
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decreasing values of k until a spanning T* is obtained for ^sbich 

I (c -k d )=o. Hie tree T* is a MEST* If > 0, 

(1,3) e a* 

thm the parametric cnalysis is performed for increasing values of k. 

It follows from Haeorem 2.2 that this metlKjd is polynomially 
bounded. 

2*6.2 Minimum Batio ?a1ii Problem 

The minimum ratio path problem is to determine a directed pa'Ki 
^ from source to sink for which the ratio ^ c. ./ ^ d. . 

(i,3) e (i,4) e ^3 

is m i nimu m. A characterization similar to Hieorem 2.4, caG be 
obtained for the MHP, The paran^tiie analysis of the shortest path 
problem can be iis^ to identify a MRP, but one serious difficulty 
arises. The parametric analysis assumes that the network does not 
contain any n^ative <ycle for soy value of the parameter in the range 
of the analysis. If the network is acyclic, this assxanption is 
satisfied and the MRP is obtained by the parametric analy^s. However, 
if the network is cyclic, n^ative cycles may arise during psirametric 
analysis and the MRP may not be obtained. In fact, it can be shown 
that the MRP problem for cyclic networks is NP-conplete, 

Hieorem 2,5 1 Oiie minimum ratio path problem is HP-complete, 

Proof ; We show liiat the longest path problem, which is a wellknown 
HP-coE^lete problem, is a special case of the MRP problem. 

Associate a ntmber . with eadh (i»3) ® -4. Hie longest path 

il 

ia a directed path £** froB source to sink for wilich i . 

(i,3)cP« 
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is maxifflum. Define c. . and d for each (i,3) e A in the 

following manner; 


c. 


1 , if (ijj) e 0(s), 


0 , otherwise, 


d..= Z.. 


(2.51) 

(2.52) 


dearly, for eveiy directed path P from source to sink 
I c = 1 , Ihus, the MEP is a directed path from souiee 

(ijj) ep 3 

to sink for which the ratio l/ ][ . is minimum. She ratio 

(i,3)ep^^ 

1/ y jl. . is minimum if aid only if Y Jt. . is maximm. 

(i,j)eP* (1,3) eP* 

Hence the theorem follows# 




GHAMER III 


C^NSIRADTE!) RiSWOEK GA.?AOia!I EXPAHSIOU ESOBiaf 

3.1 INTBODUCTION 

Hae m l nlm jia cost flow problem is to find the optiiaiaii distribution 
of flow thiou^ Idle arcs of a givai network froiB souixie to sink, Ihe 
m i n i m um cost flow probloa sxibsinaes a number of network problems, viz,, 
the loaximum flow problem, the i^iortest patii problem, the assignment 
problem, the transportation problem s^ni the transehipmait problsn. 

In this chapter, we consider tdie Tninitmuti coat flow problem wildi 
an additional linear constraint. Such situations arise often in 
practice, Dantzig [22] points out tdiat in the transportation problem^ 
the additions^, constraint appears in the form of a lower or upper 
bounded partial siaa of some of the vari^les, diaxa^ and Cooper [l5] 
suggest that in mar^ ^plications, "aie additional constraint may singly 
express the proportionalily or equalily of particular variables to force 
two routes to have equal flow. Other applications of the constrained 
imLnimum cost flow problem are the war^ouse funds-flow model discussed 
by Chames and Cooper [15] and a constrained version of Wagner's [ 92 ] 
enployment scheduling problem, 

!I3ie nr! ni.Titiim cost flow problem is a specially structured linear 
program and this structure permits to achieve many computational 
simplifications while impl«aaiting siim>les algorithm or its variants. 



47 


She additional linear constraint affects this struetiire aod maiy of 
the nice feature of the TniniTmm cost flow problQE are lost. ftjr 
instance, the basis becomes nontziaagular aid the optimum solution may 
be nonintegr^. However, the constraint matrix still poss^ses a 
structure and the atteipt is to exploit this structure as f ar gya 
possible. 

If the additional constraint is of a particular form, atten^jts 
have beaa made to tr^sfom the problem into an equivalent Tainimum 
cost flow problem of larger size [l 5, 16, 17, 22,. 60, 92 J. However, 
these transfoimations are not possible with m arbitrary ^ditional 
constraint. KLingman ^d Russell [62] developed a specialization of 
the primal simplex method for the constrained tr^sportation problem, 
which was subsequmtly ext^ided to the constrained transshipment 
problem by Glover et, al [43]»snd to the constrained generalized 
network problm by Hultz and KLingman [53] • 411 these methods 

consist of efficient ways to determine relative cost coefficients, 
representation of the entering arc and iqjdating the basis. 23he 
efjEicdency of these methods have bear ^tablished idirough extensive 
Gon^intational experience. Hie primal sin^lex approach is also extaaded 
to dev^op methods for multiconstrained (i) traisportation problem [61 ] j 
capacitated transshipment problem [isX S a*id (iii) gaieralized network 
problem [53j . Eecaatly, Masch [65] has developed a cyclic method 
of solving the transshipment problem with an additional linear constraint, 
Hit a method is not a simplex adaptation, but exploits the structure more 
efficiently, 

• reserve collect®#?? 
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Several researchers have considered the capacity expansion of a 
network vshere the bMgetaiy constraint spears in the form of 
additional linear constraint. 

A common feature of physical distribution systenssuch as electricity 
or gas supply, road or rail netvjorks is that the supply ^d demaad 
requirements change with time, Hien the existing network can not 
accomodate the increased flow due to increasing deoaa:^ requisreiiaats, 
expansion of the network is required. Sometimes, even if the siipply 
and demand reqxiiremoits do not change, expansion may be called for to 
reduce the transportation coats. In these situations, it is necessary 
to determine the optimum expansion of the network. 

One way to expand the network is by increasing the capacities of 
various aorcs. Several authors have considered the problem of optimally 
allocating a givai btuiget to increase capacities of various arcs so 
that the flow from source to sink is maximized. 5\dkerson [59 } 
considers linear capacily expansicn costs fsad proposes a labelling 
method to solve the problem. Hu [52j solves this problem by solving 
a sequQice of shortest path problems, Bansal anA Jacobscaa [ 9 ] consider 
concave e3q3ansion costs at^ approach the problem by Bender’ s decon?)osition 
procedure. When the cost functions are continuous and piecewise linear, 
the branch and botard algorithm discussed by Tomlin [ 89 ] oan be applied. 
Price [79] has also dev^oped a br^ch and bound algorithm to consider 
step ftmctions. 

An alternate w^ to expand -aie network is by introducing new arcs, 
ibllowing this ao^Pteach, Hammer [48^ solves the problem of maximizing 
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flow from source to sink sxibject to the biMgetaiy constraint bj 
pseudo-boole,^ progransming, Christofides and Brooker p9jdescil.be 
an efficioit tree seardi algoritiim using boimds calcxilated by a 
dynaEuc programming procedure. 

Some of Hxe other problene, related to the optimum expansion 
of a network, are studied by Gomoiy and Hu [ 45 J and Oettli and 

Pr^er [73 J . 


In this chapter, we consider the constrained mini mum cost flow 
problem and the network capacity expansion problsn i» a single model, 
yhich we tem as the Constrained Fetwprk CaJacity Expansion (ca?CB) 
problem, ^Phe mathematical statement of the caiCE problem is as follows: 


Minimize Z = T e. . x. . , 

(i,j)=A 


(3.1) 


subject to 


I 

h 


V, if i=s. 


1 X I X =^o, if ii^Sjt, VieH, (3.2) 

(j,i)si(i) (i.pEod)^^ 


V, if i=t. 


0 < X. . < b. .+ y. . , ¥ (i,j) e A, 
0 < J. . <a. .-b.., ¥ (i,i ) e A, 

T e. .X. . + y f. .y. . < D, 

(i,3)cA (i,j)eA 


(3.3) 

(3.4) 

(3.5) 


where x. . doaot^ the flow cm arc ( 1 , 3 ) and y denotes the amount 
by whi( 2 h the cg(pa3ity of arc ( 1 , 3 ) is increased. We assume that 
f . > 0 for eacdi ( 1 , 3 ) e A- 



50 


As an application of "ttie ®0E pixjblem, consiiier a trsaisportation 
network ere c and e denote, respectively, ttie trav^ time aid 

J-J Xj 

travel cost of transporting one unit of a conmodity over arc {i,3)» 


Each arc (i,3) has an existing c^acily which may be increased by 
incurring some additional cost. Let f > 0 denote the cost of 

XJ 

increasing the capacity of arc (i,j) by one unit, !Ehe increase in 
the capacity is further limited by an amount (a. .-b. .), Ihen, the 
problem of minimizing the total time of transportation (or, equiva^ 
iQitly, the average time of transportation) subject to the consideration 
that the cost of transportation and capacity exp^sion does not exceed 
a given biadget D, is an instgaice of the ONCE problem. 


We develop a parametric algorithm for the C3JCE problem. Properties 

of the optimum solution of the C2TGE problem aid extensions of the 

concepts of bomded variable linear programs are used to evolve the 

concept of optimmi basis structtire for liie (MCE problem. Ihe optimum 

basis structure considers the variables y. . ^d the constraints (5.3) 

1 j 

^d ( 5 . 4 ) implicitly, thesreby reducing the size of the optimum bspis 
considerably. Eurther, the constraint (3.5) is considered in a 
manner that preserves the basis triaigularity. 

Ifae optimum basis structure is xised to parametrize D, Initially, 

^ optimtmi basis structure for a large value of D is obtained. Hie 
value of B is then subsequmtly decreased and an optimum basis 
structure is laaintalned at each step. It is obvious that a given 
basis structure is optimum for values of B contained in an interval. 
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^fcSERVE coLLECimr:! 

ISie algorithm, in fact, deteimines &uch intervals and terminates #ien 
either the desired value of D is reached or infeasibili1^7 of the 
QTCE problem is discovered. 

Althoiagh, the optimum solution of the CJICE problem m^ be non- 
integer, the algorithm performs almost aH the confutations in 
int^ers. Often practical considerations require int^r^ilgr of the 
solution. !Bo take care of such situations, we suggest a sli^t 
modification in the algorlihm to obtain a near-optimum integer solution* 

We show that the CSKTCE algorithm can be used to solve (i) a 
variation of the CWCE probl«aj and (ii) a bicriteiia netvwjrk pioblaa. 
Several special cases of the CHGB problaa are also considered and 
resulting simplifications in the ONCE algosAthm are pointed out. 

lastly, we consider the capacity expansion of a capacitated 
transshipmaat network where sifplies at source nodes ah danaids at 
sink nodes are linear functions of a single paran^ter 1 . Ihe problea 
considered is to determine the least cost of increasing arc edacities 
so that the changing demands are met by the changing supplies resulting 
from changes in X . We develop an algorithm to solve this capenily 
expansion problan for all values of X for hiich a f edible flow exists* 
A special case of this probloa is also considered. 

3,2 OPTIMALITY CONDITIONS uVlIW 

Let us modi:Qr ■Qie OiOl problem by replacing the constraint (3*5) 


by the following ajnstraint: 
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I e X . T f. .y. . = _ D» . 
(i,3)eA x/x3 


(3.6) 


We refer to the modified problem as the Mbdifi^ Constrained 
Network Capacity Expansion (mGNCE) problem. 

Let us partition A into the sets X, Y, L, U ^d V such that 
X U Y consists of a spaaning tree augmented by a nontree-arc. If we set 

(3.7) 

(3.8) 

(3.9) 


=■ 0 , V (i,l) £ 1 , 


then there exists a uniq.ue solution, x. . and y. satisfying (3.2) 

Ij 

and (3.6). If this solution further satisfies 




(3.10) 

(3.11) 


then we refer to X U Y as a featsible basis af^ (X U Y, L,U,V) as a 
feasible basis structure of 1he MOTCE problem. We refer to m arc 
(i,j) as an unsaturated basic arc if (i,j) e X and a saturated 
basic arc if (i,j) e Y. We refer to the arcs belonging to L,U ai^ 

V as the nonbasic arcs at their lower, mi d dle and upper bounds 
respectively. we define an optimum basis structure of ■fiie MCiNCE probl^ 
as a feasible basis structure for ^diich Ihe associated flow is asi 
optimum solution of the MCSCTCE problem. In this section, we derive the 
necessaay and sufficient conditions for a feasible basis structuare 
to be an optimum basis structure. 
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The dual of the MGUGE pioblem is 


Maximize W = v(it -w ) - 7 b. .p. 7 (a. .-b. .) pD , 


subject to 

V"^ir « 

Pij iO a»d \j-°» ^ (i»a) eA, 


(3.12) 


(3.13) 

(3.14) 

(3.15) 


where it , p.^j 5 aad p are the dual variables s^sociated with the 

X XJ Xj 

constraints (3.2), (3.3), (3.4) and (3.6) respectively. 

The con^ilemaitaiy slackness conditions for the MGKCE problem, are 


p. . 


( 3 . 16 ) 

S 

ij 

”10- (i> 3 ) ' ^ 

(3.17) 



(3.18) 

^id 


(3.19) 


It is easy to see that Idle complaaeatary slackness conditions 
for the C31TGE problQH are (3,16)-(3»19) and the following condition: 


v(D - I e X » I 

(i,j)eA ^ (i,3)eA "-3 


( 3 . 20 ) 


Hiis aeiditional condition immediately yields the following 
r^ationship between the optimum solutions of the ®CE and MGUCE problem* 


TToperty 3.1 : The optimum solution of the MCKGE problem for D* = D 
is the optimtun solution of the ®CE problem. 



We now deriTe an altemate r^r^oataii.'On of coioplcmaaliaiy 
slackness conditions for the MGUGE problem. In the deiiTafiony it is 
assumed that p ^ Of which is maintainod throughout the algorithm. 

Yfe first prove the following result: 

Eieorem 3.1 : If ]x > o, thai any optimuffli solution, md' 

of the MOTGE problem satisfies 

- ^13) = “ • ''(1.3) = A. (3.21) 

?roof : Suppose y.. > 0 and x. . < b.. +y.. for an arc (i,j). 

It follows from (3.16) and (3.19) that 


X. . .< b. .+ y. . =^> P. . = 0 , 
13 13 *^13 13 * 




=*■ p. . = «. . + uf. . . 
X3 X3 '^13 


(5.22) 

(3.23) 


Hence 


X < b. .+ y. . and y. . > 0 =>6. .+ pf. . = 0. (3.24) 

13 . 13 ''13 •'13 13 *^13 / 

Since f . . > 0 and 6. . >0, (3.24) iiiplies that p < 0, v^ch 

contradicts our assuaqjtion that p »■ 0. 


iSieorem 3.1 inplies liiaT in the optimm soluld-on of the MGNCE 

problem if y. . >0, thm x. . =b..+ y. ^d if x. .<b. .+ y. thaa 
_13 13 13 ''13’ 13x3 •^X3* 

y. . = 0. Using these properties, the conditions (3.16)-(3.19) can be 
X3 

^own equivalent to 
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'r’i" “ij* ='±3' °> 2^i3^> * (3.27) 

°i3"""i3 ' V’'i“=13*''(\3^13> “• "i3^13’^i3'°’ 

lt(i,j)eA, (3.28) 

I 

TT^W. >c +u(e +f ) — > X =a =£» h 

j i ij ^ ij \r ±3^ 

7(i,j)eA. (3.29) 

It is clear from (3.25)-(3.29) that every optimum solution 

satisfies 

Ji^ = max. { 0 ,Xi_bi^ } , ¥ (i,j) e A. ( 3 . 3 o) 

Elis relationship indicates that in the optimum solution, y. . 
can be considered iiqjlicitly. In view of this r^ationi^ip, we 
restrict our attention to x. . only. It further follows from (3.25)- 
(3.29) that the necessary and sufficient conditions for a feasible 
basis structure (X U Y, L,U,V) to be an optimm basis structure are 


that there exist numbers ir. satisfying the following cwaditions; 

3 


’3 - *i ‘ °i3 * " ®i3 ■ ’ ^ ’ 

’3 -’i “ °id *"(^03 ■^13^’ ’ (^’^> ' 

, ¥ (i,j) e L, 


TT . - V. < c. . + p e. . 
0 1 — 13 


c.. + p e. . < Tr.«Tr. < c. .+ uCe. .+f . .), ¥ (i,i) e U, 

ij 3-1 — J X- 11 13 13^* V ’ 

V. - ir . > c. .+ p(e. .+f. .) , ¥ (i,3) e ¥. 

3 1 - 13 13 13^ ’ ^ 


(3.31) 

(3.32) 

(3.33) 

(5.34) 

(3.35) 


Ihe set of basic ai«s, X U Y, consists of a spanning tree, 
s^ I, and an additional arc, say (p,<l). X = X TIT and Y = Y fl T, 
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Gonsequoitly, we may r^resaat XyY asXUYU ^(p»‘l) }• refer 
"to X u Y as "the basic tree and arcs belonging to it g® tree— arcs « 

Arcs not belonging to the bgisio tree are referred to as nontr eeu.arcs , 
ilirther, we refer to (p,q) as the pj-ypt arc . It may be i*®aeE4)ered 
that the pivot arc cm be an unsaturated or saturated basic arc. 

We now use the special structure of the basis to sii!«>lif:y (3.31 )-(3.35) , 


c d 

Let IT , ir . ir. be the nitdaers satis^Ting 

3 0 3 

= 0 and Tr° - it° = c , ¥ (ijj) e X y Y, 


0 and IT . « ir = 4 


ij 


, ¥ (i,j) e X, 




IT =0 aid TT._ tr, = 

s 3 1 






(3.36) 


(3.37) 


(3.38) 


where p is a positive real nunher. 


It is clearly seen that 


TT . = tr? +p TT*^ , ¥ 3 e H, 

3 3 3 


(3.39) 


Let X = 0 or b or a^„. Let I , U and ¥ be the 
P<J P<1 1><1 

sets of nontree arcs having x. . = 0, b . and a. . respectively. 
Ibr each arc (k,l) ® L y 5, define liie numbera Cj^ and 3^ as 
follovras 

(k,l) e L , 

’^l - ’^k - ^ki » ^ ^ » 


°KL “ ^ 


(3.40) 
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} _d 


1 , 


5^ » if (k,l) e L , 


’'l “ ""k ■ ®]jl » ^ ^ 


(3.41) 


Similarly, for each arc (k,i) e U U V, define the nmnbeis e 


kl 


and as follows; 




*k-*i “kl ' “ “ 5, 


(3.42) 




(3.43) 


’k - *? * “kl * ' 5, 

7l “ ’^k " “kl " ^kl* (4^1) ® V, 

It immediately follows from (3.4o)-(3.43) that 

^kL= -®kL ^ ^kl» ^ ^ ^ ^5.44) 

Substituting (3.39) in (3.3l)-(3.35), md then using (3.4o)-(3.43), 
we obtain the following equivalait conditions: 

Cp^+v5p^= 0, if (p,q) is ai msatuiated basic arc, 


c „+Vd = 0, if (p,q) is a saturated basic arc, 

~pq “Pq 7 7 


(3.46) 

(3.47) 


We ref er to the conditions (3.45 )-(3. 47) as the optimality 
conditions for the MCNCE problem. It is clear from (3.45) that given 
the basic tree, -Hie value of p d^ends upon liie arc (p,q). We refer 
u as t^e pivot ratio. 
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3.3 DEVELOEMEST OF TEE ALGOBIUm 

In tMs section, we describe in detail the deTdopment of the 
algorithm for the ®CE problem, Hie algoiiliim treats D* as a 
parameter aid solves the MCSfOB problem for continuously decreasing 
values of D’ , It first obtains eo. optimum basis structure of the 

MGffGB problem for a sufficiently large value of B’, say B , It thm 

o 

determines B^ such that this b^is structure continues to remain 
optimum for all B* • e (B^ , B^), A dual simplex iteration is performed 
for B* = B^ to obtain an alternate optimum basis structure which m^ 
allow further decrease in B*, Ihe algorithm proceeds in this manner 
until either an optimum soluticai of the CSfCE problem is obtained or its 
infeasibilily is indicated. 

Given an optimum basis structure (X U Y U {(p»<l) ) > B,U,V), let 


us classify the nontree-arcs, E U B U V, as follows; 

S =t(kjl) e i U U : dj^ < 0 aad c ^ 0 } , (3.48) 

S ={(k,l) e U (J ? : dj^ < 0 aid Cj^ > 0 } , (3.49) 

G ={(k,l) e i U ^ ^0 and c^^ <o } , (3.50) 

G ={(k,l) e U U Y ; d < 0 end <.0 > , (3.51) 

H ={(!£:, 1) e B UU : 5^ > 0 } U{(k,l) e tj U Y : d^ ^0 } (3.52) 

We refer to -Hie sets S U S, G U G and H as 1die sets 
of active, critical yard passive arcs irespectively , Hdese sets are 
mutually exclusive well as exhaustive. 



3.3.1 Characteilstlc Interval 

let the optimum basis stmctu3re of the MGNOE problem for D’ = D 
be (X U Y U {(p,q) } ,L,U,V)i Here D is su<ii that (X U Y U {(p,<l) } , 
1,TJ,¥) is not 3h optimum basis structure for D' >5* Let 5. . be tiae 

X.. J 

optimum flow Z = ^ c. . x. . , It will be shown, in Sectaons 

3*5.2 and 3.3.3 that if (p,q) is an unsaturated baoic arc, "Bien 

(p.q) e S, and x = 0 or b , It will also be iSaowu there that if 

* pq pq 

(p»q) is a Saturated basic arc, "aien (p,<l) e J, and x = b or 

a • We now determine the interval (D-5) such that "aie current basis 

pq ^ ' 

structure continues to remain optimum for all D s (1,5) , Ihis intervsl 
is Ijnown as the characteristic interval associated with (X U Y U {(p,<l)}, 
L,U,V)* 

Since the numbers Cj^, d^^, aad d^ are uniquely determined 

for a givm b asis structure, the optimality conditions are not affected 
by change in the value of D' . Itowever, in order to satisiy (3.6), flow 
must be changed keying the basis structure intact. 

In the basis, there is a cycle consisting of basic arcs vhach is 
formed vhen (p,q) is aided to the basic tree, let be the set of 

arcs in this cycle, Ohe flow can be changed by circulating some 
additional flow in Bae direction of additional flnw in Wp^ 

d^sads upon the status of arc (Pf<l) iihich could be either saturated 
or unsaturated. If (p,q) is an unsatuiated b^ic arc, then 
(p,q) e S £(1 U U). Haace to ke^ the arc (p,q) unsaturated, the 
additional flow must be circulated along (p,q) if (p»<l) e i* and opposite 
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to (p,<i) if (p,q) e U. Similarly, if (p,q) is a Saturated basic arc, 

(p»<l) e ^ £,(U U y), ®iis iiEplies that the additional flow 

must be circulated along (p,q) if (p,q) e U and opposite to (p,q) If 

(p»'l) ^ V. After defining the orientation of W accordingly, let 

pq 

^PQ. "^<3. sets of arcs in along and opposite to its 

ori station re^ectively« !Ehe maccimum amount of flow, denoted by w, 
that can be saat in W is calculated by xxsing the bomd restrictions 
of the basic aros, i.e,, (3»lo) aid ( 5 .II). Let us define the nundier 
w . for each (i,j) e W as follows: 


w. .= 
11 


‘’id - \d ’ “ ^ ^ *pq • 

. if (l,d) ' ^ 

%d ■ ^id ’ ^ ^ *pq' 

.^id ' '’id ’ ^ "A-q’ 


( 3 . 53 ) 


•vdiere X gaad Y are the sets of unsatvur'ated al^ saturated basic arcs 
in X UY U {(p,q)} respectively. 


!Ehe w. . deletes the maximum ^ditional flow peimitted by 
11 

any arc (i,l) e W without violating its bomd r^trictions. Hence 
' pq 


w = min. { w. . } . 

(i,d) e w„„ 


(3.54) 


5 + w d , if (p,q) is an imsaturated b£eie arc. 


D + w d • if ip%^} iB a sa:turafBd bgpic arc* 

— 


( 3 . 55 ) 
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Ihe follov?ing theorem establi^es the optimum flow of the MCSIGE 
p3x>bl©Qi for all I>» e (B,5). 

3.2 : ibr all D« e (D,D), the optimim flow of l^e M®GE 
problem is given by 


X. . = 


8 w , y {i, 3 ) e w_ , 


P<1 


V ^3 - 8w , ^r(i,3) , 


X. . 

13 


» ^ (1,3) 9^ Wp^ , 


aid the value of the objective function is 


(3.56) 


where 


Z = 


2 + 0 Cp w j if (p,q) is an msaturated 

basic ai*c, 


jz. 


0.2pq ^ , 11 (P»l) Is a Saturated 

basic arc. 


6= (C- B')/(5 - D). 


(3.57) 


(3.58) 


^luof : Let (p,q.) be an unsaturated basic arc. It can be easily 
shorn using (3.36) and (3 .40) that 


c = I c ^ c 


(3.59) 


aimilarly, using (3.37) aid (3.41 ) it can be shorn that 

d = ( T e. .» T e. .) + ( J f 

(^»3) ^ (i,3)elpq (1,3)^ Wp^ n Y 

”, J ^5.60) 

Hoice c and d denote, respectively- the increase in 
pq pq ' —VI 

the value of Z aid D' if unit amount of additional flow is 
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circulated in along its orientation. Since (p,q) e S, we 

have 0 aid c^^ >0. 53iei«fore Z increases aid D' 

decreases as additional flow is soit in W . Since u = -c /S > 0, 

pq pq P<1 ’ 

it denotes the increase in Z per vaiit: decrease in D* , 


ibr D* £ we see that 0 5. 8 5_1 and haace the flow 

^ij* Sivaa by (3.56), satisfies the flow requirODiaats of 'fee optimm 
basis structure. We now show that the condition (3.6) rosiains 
Satisfied in the interval (2»^) • Substituting x. . from (3.56) in (3*6)’ 

JLJ 

R 4 H.S. of (3.6) = T e. . X. 0 w ( J e. . 

(1,3) L 


I e. .) »ew( I f . I _f. .) 

V,i)eSp,nY 


= - D « 0w d (using (3.6o)) 

pq 

= _ D* (using ( 3 . 55 ) and (3.58)) 


(3.61) 

(3.62) 

(3.63) 


fhis proves the optimality of prove (3.57) consider 

E.H.S. of (3.57) = I o.. !„+«»( L °ii-, ,l “n) 


= Z + 0 w c (using (3 .59) ) 
P9 


(3.64) 

(3.65) 


The proof of the theorem is given for the case Mshai (p,q) is 

an msaturated b^ic arc. If (p>q) is a saturated basic arc, thai 

c aad d dmote, respectively, the increase in the value of Z 

~pq "P<l 

and D* vtien raait amount of additional flow is circulated in 
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aiong its oiioataliion. !I!he iaterpa^taliion of p rsaains the saEis* 
Further, optimality of x aad Z csiP. be prored in the similar 

Xj 

majaaer* 

Let ^ be the optimum flow for D’ = D. Further decrease 
in the value of D’ is blocked by an arc («,6) for which - w. 

If (a, 6) e X, then g =0 or b^^ g j and if (o,6) £ X, then 
— a 6” ^a8 ®'a 3* decreased without changing the 

optimum basis structure, flow in arc (®»6) violates its bomda . Sius 
a dual simplex iteration is performed to obtain an alternate optimum 
basis structure for = D, which may allow furcher decease in liie 
value of D* , 

3»3.2 Dual Sin^lex Iteration 

Ihe alternate optimum basis structure for D* = D is obtained 
by performing a dual singles iteration, Hi tbe pai^eters ^d sets 
associated with the optimum basis structure (X U X U {(Pj<l) } t 
are indicated with a prime for the alternate optimum basis structure. 

The du^ simplex iteration is perfom^ by dropping the arc 

(a, 6) from the basis md sheeting a nonbasic arc to enter the basis. 

When the arc (ot,0) is dropped from Ihe basis, the rmaining basic 

arcs constitute a sp^ning tree, 3Ms spanning tree becomes the basic 

tree in the alternate optimum basis structure. If (o, 3) = (p»<l)» "tlien 

X* = X aid X' = X. In -aiis case, the values of c^, 5^, and 

d» raaain unchanged for all nmtree-arcs other thaa Idle are (p,q), 
“id. 

The values of c , d , c and d are multiplied by -1 to obtain 
pq pq “Pt 
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-^q / (p,q) , then Idae basic tree 

changes. Hie new basic tree is givm by X* uY‘ = X U T U ^Cp, 9.)J- 
{(a, 3) } . Shis basic tree consists of two subtrees 2’ (ccmtaining 
source) and 2^, and the arc (p,q) wfaicii joins these tvuo subtrees, 

2he values of c^, d^, c_^ and also get chenged. 3b coii?>ute 

these values j let us classify the nontree arcs L* UU* U ^ as 
follows; 

B° = {(k,l) e I> U 7» : ke2» and 1 e 2|; or k e 2^ and 1 e 2p, 

( 3 . 66 ) 

s'* = {(k,l) e L» UU« U V': k e 2‘, 1 e 2^ and (k,l) e £« U V; 


or ke2^, le2;j aid (k,l) e U* > (3.67) 

= {(k,l) e L» U U» U V»; k e 2», 1 e 2^ and (k,l) c U»i 

or k e 2^, 1 e2» and (k,l) e L» U ^ > . (3.68) 

Using (3.36)-(3.37) md (3.40)-(3.43) it can be easily ^wn 

that c* , d* . £* and _d' can be expr^sed in teians of c » 
iCJL iCL "iCL Ai. 

d^, and djj^» depending upon the status of the arc (p,q). 33iese 

ejqiressions are as follows: 


Case I ; (p,q) is an msaturated basic eiro 




Cj^, ¥ (k,l) e #, 
°lcL'^°pq> * 


( 3 . 69 ) 
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^ (k,l) £ 




(5.70) 


^kL = 


c^ , ¥ (k,l) e iP, 

^la^pq- * 


(3.71) 


^k 


-kL» 


¥ (k,l) e ^ , 


^kl%q’ ® » 


(3.72) 


Case II ; (p,q) is a saturated basic arc. In (3.69)-(3.72), r^lacing 

c and d by c and d , respectively, we get the 
pq pq / -pq -pq» ^ ^ 

desired expressions for c « , d« , c' and d* . 

KL* kL’ — kL — ki 

Let S’ U S’ doaote the set of active arcs. If S* U is 
empty, it follows from IheorOE 3.6 that the MQfCE problem is infeanible 
for all D’ < D. ifowever, if S’ (J is noneirply, define a nxariber 
^k each (k,l) e S’ U _S’ as follows; 


^k~ ^ 


-5k^k- i'' 0=-^) ' ®'- 


-ck/^k- 


(3.75) 


Clearly, > 0 for eadi (k,l) e §*Us’. Let (g,h) € S’ U S’ 


be an arc satisfying 


I = ndn. { p* - >• 

^ (k,l) e S’ U S’ ^ 


(3.74) 
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iHrrtiierj let 1* j IJ* aid ip be the sets of arcs in 
E' U U ^ — {(g>h)} at tbeir lower, middle and i:q>per bounds 
respectively , We now state sojne liieorenis along with liieir proofs 

to show that (X’ u Y» u Ug,b)} , is m alternate optimum 

basis structure for D* = D with (g,h) as gai unsaturated baisic arc 
if (g,h) e S' j and as a saturated basic arc if (g,h) e _S' . In 
theorems 3.3, 3.4 ^d 3.5, it is asstimed that (p,<l) is an msaturated 
basic arc. On similar lines, results can be proved with {p,<i) as a 
saturated basic arc. 

theorem 3.3 . No arc (k,l) e L» U U» U ^ is critical. 


Proof . If 
v/e get 

If (k,l) e E^, thoi substracting (3.45) from (5.46) aad using (5.69) 

£hd ( 3 . 70 ) we get ( 3 . 75 ). If (k,l) e th«i adding (3.45) to (3*46) 
and using ( 3 . 69 ) and (3.70) we again get (3.75). We have thus obtained 

c^+ Vd^ t.0, ¥ (k,l) e Lt u ij« . (3.76) 


(k,l) e eP, then substituting (3.69) aid (3.70) in (3*46) 


o? + V d» >0. 
kL kl - * 


(3.75) 


Similarly, it can be shown that 




(3.77) 


If an arc (k,l) e &» and is critical, then (k,l) e (£' U H') 
and d^ < 0, c^^ £0. Since P > 0, we hasre e*^^+ n d«^^ <0 


kL 


kL 


kL 


which contradicts (3.76). likewise, if an arc (k,l) ^ 5' and is 
critical, then (k,l) e (U» U and d^ < 0, c^ <.0. Hius, 
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~K1 ^ ^ -KL ^ ° contraflicts (3,77). 

Let y* denote the value of the pivot ratio for the altieinate 
optimum basis structure. It follows from (3.45) that 


h' = S 


-c^/d^, if (g,h) is sn msaturated basic arc, 
“—gh^'^gh’ (s»^) a satiirated basic arc. 


(3.78) 


ConiJaiing (3.78) with (3.73) we find that 


= '■V- 


(3.79) 


Ihus the pivot arc is m active arc which yields the lowest value of 
the pivot ratio. We now show that the value of pivot ratio is 
increased by the dual simplex iteration. 

theorem 3.4 ; p' ^ p • (3.80) 

Proof : Consider any arc (k,l) e S'U^*. If (k,l) e S*, thoa 
(3.76) can be written as 


u < -c* / d’ = p’ . 
— Id. Id. a. 


If (k,l) e s», then (3.77) can be written as 
V ± " b] 


‘k- 


We have thus obtained 

V 1 % , V (k,l) e S‘U s» . 

Since (g,h) e s’ U ;S’, we have 



(3.81) 


(3.82) 


(3.83) 


(3.84) 
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^h^:rgDi_5.;.5 : ae basis structure (X* U T’ U {(g,b) } , X’ , U', 7») 
satisfies tbe optimal i ty concLiticMis » 


dual sii^plex iteration is essentially a degerate iteration 
and flow on any arc does joot change, Therefbre the value of D< is 
not affected by this iteration, She optimality condition (5 *45) is 
satisfied in view of (3 ♦78), 5b prove "ttiat liie remaining conditions 
are also satisfied, consider aty arc (k,l) e S’ y^* , Qlearly 

^ s«, 

v' £ vk " ^ ^5.85) 

Ij-k'^-k » (k,l)eS’, 

which after rearrangemoat reduces to (3,46) sad (3,47), Now ojnsider 
any arc (k,l) U ^ . Since there is no critical arc, (k,l) must 
be a passive are. Hoace d' ^0 if (k,l) e L’UU», aid d» ^0 
if (k,l) e U' u V*. Multiplying p’ ^ p by d|^ aid adding c^ to 
both the sides, we get 




( 3 . 86 ) 


which reduces to (3.46) using (3,76). 


Similarly, multiplying y' >_ p by d^ and adding c^ to bctti 
the sides, we get 



^k^- 


h d* . 


(3,87) 


which reduces to (3.47) using (3.17). 53ais completes tiie proof. 

Sheorem 3.6 t If S’ U is ea^ty, laie MCNCB problem is infeasible 
for all D' < D. 
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• I'®'* be -ttie ndnimm value of D» for #iich the MOJGE 

problem is feasible. dearly, D . is the optimum objective 
function Value of the following linear program: 


Minimize ^ e x + T f v 

subject to 


( 3 . 88 ) 


^ X ^ X = < 

(j,i)ei(i) 3 ^ (i,j)e o(i) 


-V , if i=B, 

0 , if i/s,t, ¥ ie M, 
T , if i=t, 


0 i^ii . » (i,3) ^A, 


“-5^13 -\3 -‘ij ’ 


(3.89) 

(3.90) 

(3.91) 


The dual of (3.88)_(3.9l) is 

Maximize (ir -ir )v _ T b. .P _ Y (gi ._b. .) , (3.92) 

* = (1,3) c A (i.3)> i^ 


subject to 


V’i-"i3 -^13* ^ ’ 

^13 - *13 ^^13 ’ ^ > 

Pij - 0» ^ (i»3) ® A, 


(3.95) 

(3.94) 

(3.95) 


and the complemaataiy slaclsness conditions are 


'’I3^\3"^i3-='i3’ ' 

^13(^13" 


(3.96) 

(3.97) 

(5.98) 

(3.99) 
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We now show that if S' U^* is esrpty, the solution x.. md 

- 0-3 

^3 ~ niax. Co, } solves (3.88)-(3.9'5) and hence D = . 

let Tj- be defined by (3.37). Set w . = y ^ e 1. Hence 

"■ ®i3 " ° ,V(i,3)eX', (3.100) 

’'i**''3 ^ ®i3 " -% * ^ e Y«, (3.101) 

Since the set of active arcs is eaiply s^d there is no critical 
arcy all the nontree—arcs are passive arcs# H^ce 

% 10» ^ ^ U U', (3.102) 

10, ¥ (k,l) e ii' U 1^. (3.103) 


Using (3.41) and (3.43), these conditions become 


Vd ■" ®i3 ’ ^^*3) e 

‘ Vi-’i ®i3 ^ ® 

ir . -ir. + e. . < _f. ., ¥ (i.i) e F . 
13 13 - i 3 ’ 


let us further define the varLdales P. . and d . as follows: 

X3 13 


0 , ¥ (i,3) e X* y i', 
p =V-Tr.-e W(i,5)e U' U ¥», 


13 


3 X X3 

Eli 


, V (i,3) e X’, 


ff-ir.-e. .-f. :, V (i,3) e ¥> , 
jj 1 13 X3^ ^ 

^0 ,¥(i,j)rF. 


(3.104) 

(3.105) 

(3.106) 


(3.107) 


(3.108) 
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It cao be easily seen that tt . p. . and ^ » as defined above, 

3 13 i3 

give a feasible solution of tiie dual pioblem (3.92)-(3,95), v*idi 
together with and satisfy the conplementaiy slackness 

conditions ( 5 , 96 )-(3. 99 ). Hmce and is an optimum solutitm 

of the p3x>blem (3.88)-(3.9l) and 2) = D . , 


iron 


5 . 3.3 Initial Optimum Basis Structure 


I'St be the highest value of D' i^jto vdaich the a>nstiaint 

( 3 . 6 ) affects the optimum solution of the MdTCE problem. le now outline 
a procedtire to obtain an optimum basis structure of the MGSGE problem 
for D’ = D , Justification of the procedure is givKi subsequently. 

Step 1 ; Solve the following mini mum cost flow problon by any primal 
sin?)lex method [ 10 ] : 


Minimi ze 


subject to 


I 


C. , X. 


(i,3)sA 


( 3 . 109 ) 


—v , if i=s 


I l ^ 0 . if » i ' H. 

( 1 , 3 )^ 0 (i) ^ U,if i=t, , 

* * (3.110) 


° i='i3 - “i3 ’ 


(3.111) 


Let X. . be the optimum solution. In this solution^ let 
i3 

and tP be the sets of basic arcs and nonbasic arcs having 

x?. = 0 , b.. aid a>j respectively, Ihe B° consists of a spanning 
X 3 * 10 ij ^ 

tree, 13° is eiii 5 )ty and U consists of nontree-arcs , Let 

X° = {( 1 , 4 ) e B° Y° = {(i, 3 ) eB°:b^^ ^ ^ij - \j 
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lalojig Z y y as a basic tree, confute ^ saad txsijag (3.40)- 

■ K . L iK.1., 

(3.43). Go to Step 2. 

^ * I*®"*! G° U G° be the set of critical arcs. If G° U G° is 
eu^ty, go to step 4; otherwise select m arc (p,q) ® 5° U G° 
and go to step 3, 


Step 3- ‘ bst be the cycle formed when, (p,q) is added to the 

basic tree. Depending upon the status of arc Cp,q), define 

the orioatation of W ^ as described in Section 3.3.1. 

Determine i using (3.53) and (3.54) md send additional 

flow w in W along its oiientation. Drop aa arc 
pq 

(a, 6) e W for ishich w„d= w. !I3ae new basic tree is 
Z^U Y° U {(p,q)> - {(a,e)> . Update if, and 

recompute c^, d^, c^ and d^. Qo to St^ 2. 


Step 4 : Let x. . be the existing flow, Sien D = T e. .x, . 

(i,3) e A 

is Idle optimvan flow for the C3JCE problem. Otherwise let 

‘•"•O O CS' 

S U ^ be the set of active arcs. If S U ^ is em^jty, 

thm the CSliCE problem is infeasible; otherwise define a nunher 

V° for each (k,l) ® U as 

KlL 


u° = 

Id. 


(k,l)e r , 

-4/4 , if (^1) ^ S°, 


( 3 . 112 ) 


and select an arc (g,h) e S° U S*^ for which 4” ®^*_o o'^KL 

^ (k,l) e S U S 
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jOhe initial optimum basis stiuetuie for D* = D is (X° U U 
1 ,U ,7 ), where are the sets of arcs in U 5° U ^ 

" { (Sji^) } having flow equal to Oj siid i:*e^ectively* 

In Step 1, an optimum solution of the MC3NGE problem wi'SKJut the 
constraint (3 •6) is obtained. According to our notations, this solution 
satisfies the following conditions; 

^ (is:.!) e U If and c^LO, ¥ (k,l) e ^ U ^ . (3.113) 

In Steps 2 and 3, m optimum solution of (3.109)-(3.111) requiiing 

minimum value of D' is obtained, Siis is done by selecting an arc 

(p,q) e 5° for which c^^ = 0 and d^^ < 0, or an arc (p,q) e for 

which _c^^ = 0 and < o, and performing the primal a imp l ex pi vat 

iteration. In view of our interpretation for c° , d° , e° and d° , 

pq» pq» -pq -pq» 

the value of Z does not change but the value of decreases by the 

pivot iteration, Since c° = c° = 0, the conditions in (3.113) are 

pq ~p<l * 

maintained. At the end of Step 2, a solution is obtained in i/4d.ch 

can not be decreased without increasing Z, Haace this value of D’ is 

the value of D . It is obvious that if D > D , this ajlution 
max max' 

optimizes the OTOE problKa, 

The optimum basis structure for B* = B is constructed in 

max 

Stqj 4, If the set of active arcs is "aien it follows from 

Theorem 3,6 that the MOTCE problem is infeasible for all B' < B . 

... 

Hence the OTOE probloa is infeasible if B <B . If the set of active 
apcs is non empty I then an gyctiire axo vAmch yi^ds the imnlinuia value of 
the pivot ratio is selected a3 the pivot arc^ We now ^ow that this 
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Erom ( 3 . 45 ) we get Goasider arc (k,l) e S° {j _s°. 


Clearly 


V® = 11° < 

gh -V 




a' ki 


- 4/^1 » ^ S °* 


(3.114) 


-Id -kl 


Alter reatrsaagemetit (3.114) reduce to (5.46) and (3.47). Now 
consider any arc (k,l) ^ (J _S°, Tenuination of Step 2 guareatees 
that there is no critical arc. Hence (kjl) is a passive arc. 3his 
implies that 

d^^O, if (k,l) e L° U 0^5 and d^>0, if (k,l) eD®uiP. (3.115) 


In view of (3.113) ^ad (3.115), the optimalily ccfflditions (3.46) 
and ( 3 . 47 ) are satisfied by the arc (k,l). 

3.3*4 Hiysical InterpretaticHX 

iEhe algorithm first obtains ah optimm solution of the MCSFGE 
problem for a large value of D* . SMs value of D* is then conti- 
nuously decreased and optimality is maintained. At each stqp, the 
algorithm keeps a basic tree and the pivot arc is used to decrease D* 
by performing a simplex-like iteration, ill e value of Z increases 
continuously and the pivot ratio indicates the rate of increase in Z 
per unit decrease in B* , ihe algozitiHa proceeds in this manner until 
B' = B, vhaa the optimum solution of the MC2fCE problaa is the optimum 
solutim of the C3TGE problem. 

The algori-liim is essaitially based on the trade-off of B* vdth 
Z, The' pivot, arc plays the caatral role in this trade-off. The pivot 
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arc is selected among the nontreewaxcs which are cl^sified in 1 » three 
cat^oiles, viZij ciltical arcs| active arcs aad passive arcs# Critical 
arcs are those nontree-arcs which lead to decrease in, D’ wilhout 
increasing Z if any one of those arcs bec»mes /^e pivot arc. Active 
arcs may also lead to decrease in D* but by increasing Z, ISilike 
critical md active arcs, passive arcs do not leai t© decrease in D* | 
HxLle obtaining ttie initial optimum basis structure, critical arcs are 
used to decrease D* , Since critical arcs do not ^pear in later steps, 
active arcs are then used to decrease D* , If there is no active arc, 
the value of B’ can not be decreased aid the algorithm terminates. 
However, if there are active arcs, thai ^ active arc which yields m ni TTBim 
value of the pivot ratio is sheeted as the pivot arc, !Ihis selection 
criteria is quite intuitive as it iu^ilies that at every st^ D' is 
decreeised by incurring Tnin-intiiTn increase in Z, In srhsec^mt iterations, 
value of the pivot ratio increases indicating that decreasing D* 
becomes costlier, 

Z,4 DESCRIPTION OF TSE ALGORIWi 

A 

Step 1 : 


formal description of the GKCE algorithm is given below. 

Obtain the initial opIhjiBim basis structure (X U Y U {(p,q)} t 

L.U.V) for D* = D by the procedure described in 

' ' ' ■ 

Section 5,5.3, let x. . be the optimum flow. Set 

- Z = 1 2^4 and 5 = B . If S ^5, go to 

( 1 , 3 ) e A 

S 1 #ep 5f otherwise go to Stop 2# 
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St ep 2 


\q cycle consisting of basic arcs* Define 

the orientation of W along fn ol i-p n j / \ 

pg axong xf 0 md (p,q) 

is an unsaturated basic arc, or x 


' ~P<i ^pq 

a Saturated basic arc. Ottieiwise define -Sje orientation 

of V opposite to (p,,). lo, s ^ ^ ^ 

pq — pq 

arcs in along aad opposite to the orientation 

^pq respectively. Define a nunijer i. . for eac^ 
Vijj; e as follows: 

^ij”^ij * (^»^) e X , 

w. .= ^ ^ ^-^q » 

^ij"*ij » ^ ^ ^ » 

e i , 

da ere X and Y are the sets of unsaturated aid saturated basic 
arcs respectively, 

Oonpute 


aid 


w 


min. { w. . ) , 

(i.j) t 


? = min. {w, (DJD)/ d } 
^ pq 

I^date 2 and D as follows: 


ii 


r.. + ? , ? Ci,3) e W 


? Cm) e 


2J 

X. . « f 
13 

X. . 

_13 


pq 


”pq 


» ? Cm) f % 


pq 
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2 + r c 


2 = 


pq' 


. Z + r e 4 
L. -^q» 


D = < 


I) + r d 9 
PT 


S + r d 9 

-pq’ 


(p»<l) is an msatnirated basic axcj 
ii (p»<l) is a Saturated basic arc, 
if (p,q) is an unsaturated basic arc, 
if (p,!!) is a saturated basic arc. 


If C — D, go to Step 5; otherwise go to Step 3^ 

Step 3 : Idoitify an ai^c (a,g) g for which ^^^0= w. The 

new basic tree is X LJY U{(p,q)} - {(a,e) } * RecoHpute 
^Id.’ ~kL 1‘et S U_S be 1±ie set of ^tive 

arcs. If 5 U_S is enply, go to Step 4; otherwise define 
a nunher for each (k,l) e S U ^ as follows: 

-V4i ’ “ 

and select an arc (g,h) e S U_S for vshich ndn. 

^ (lE,l) e s U ^ 

{Vj^} . Set (p,q) =® (g»ii), vdiicdi is ah rmsaturated basic 
arc if (p,q) e S, gtid a saturated basic arc if (p,q) e S. 
Update the basis structure and go to St^ 2. 


Step 4 : Ihe OTOE problem is infeasible, STOP. 

Step 5 : Ihe solution and mar. {0,x^^-b^j } is the optimum 

solution of Ihe GtrCE problem with the value of objective 
function 2. STOP. 
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3.4.1 The Txadeott Curve 

Hie CSJGE algorithm cssi be iised to solve the MCNCE problem for 

all feasible valu^ of B’ , If the optimum obj ©itive function value 

of the MCKGE problem is plotted gainst D* , a piecewise linear convex 

function is obtained. In this curve, the slope of a linear segmait 

is n^ative of the value of the pi -rot ratio in the associated optimum 

bahis structure, Olie slope of this curve is aero for B* > B , 

max 

Ehe slope decreases as B decreases and finally becomes — » fer 

B = B . , Since this cuirve daicts the tradeoff between B* and Z, 
mm ^ 

we refer to it as the tradeoff curve . 

3.4.2 ETear-Optimum Integer Solution 

!Phe optimum solutions of the MCKGE problraa obtained by the 
algorithm at the end points of the characteidstd-c intervals are integral. 

If the prescribed value of B is not one of Ihese aid points, the 
optimum solution of l^e OTCE problem provided by the algorithm may be 
nonint^ral. In such cases, the following modification in the 
algorithm yields a nearuoptiaam int^er solution. In the last iteration, 
whm the pr^cribed value of B is met exactly, r may be less than w 
and m^ be noninteger. Hiei, increm si tin g r to "ttie neap:‘est integer, 
i.e., f r + 1 3 t yi^ds an int^ral ais well as a feasible flow. Ihe 
difference in the objective fmetion values of 'ttris integral solution 
and the optimum solution is ( [ r + 1 } — r ) if (p»4) is an ira- 

saturated basic arc as^ ( ^ 3 - ^ (p»^) d.s a 

saturated basic arc. If fee rabie of Cp^ or Cp^ is saaall, the int^ral 


solution is near-optimum, 
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SiS NUMERICAL EXAMPIE 

In this section, a nimeiical exaaple is solved to illustrate 
various steps of the CHGE algoiiltoa, She network is shown in ELg. 
Ilte data associated with arcs is given in fable 3.1. ®ie nodes 1 
and 6 respectively rep reseat soiirce and sink. She flow to be 
transshipped from source to sink is 12 units^ It is desired to 
obtain the tradeoff curve of the numerical exaa^le. 

Table 3.1 i Data of the nmerical exaa^le 


(i>3) 

"ii 

b. . 

c. , 

3-3 

e. . 

3-3 

f . . 

3-3 

(1,2) 

12 

6 

2 

0 

1 

(1,3) 

8 

8 

3 

0 

0 

(2,3) 

5 

3 

2 

0 

1 

(2,4) 

10 

5 

4 

0 

2 

(2,5) 

<x> 

5 

1 

0 

3 

(3,5) 

10 

9 

2 

0 

4 

(4,6) 

12 

8 

2 

0 

3 

(5,4) 

5 

4 

5 

0 

2 

(5,6) 

15 

6 

2 

0 

2 
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!Qie ®GE algoii'ttaa obtains the tradeoff curve in sdLx iterations, 
iiae coii?>utations involved in these iterations are summarized in fable 
5.2. In the table, the te3mi 'status* indicates the set to which sn 
arc belongs. fhe s3mabol + denotes the pivot arc and + denotes the 
basic arc leaving the basis. Basis in various iterations are shown in 
Jig. 3 *2, In this figure, the pivot arcs are indicated by dashed 
lines and saturated baisic arcs are drawn in double lines, Hie tradeoff 
curve is d^icted in ilg. 3*3* 

d,6 BELATED NETWOBK PROBimS 

In this section, we consider (i) a variation of the ®CE problaaj 
and (ii) a bicriteria network problen. We ^low how the ®CB algoriltoi 
can be used to solve "feese problems, 

5.6,1 A Variation of the GICE BrobloB 

In the application of the ®G1 problem cited in Section 3,1j 
our objective was to Tninlinize the total time of transportation subject 
to the consideration that the cost of transportation aid capacity 
expansion does not exceed a prescribed budget D. Our objective caauld 
ah well be to minimize the cost of trgnsportation and edacity expansion 
subject to the consideration tbat the Iptal time of transportation does 
not exceed a prescribed value. -Problems of this kind can be formulate 


as follows: 
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Table 3.2 : Solution of the numeiical exsoijle 


Ite- 

Z D 

ion 

No. 


w 


Sree-arcs 


Hontree-a3?cs 


1 . 60 59 0.5 


(k,l) 


°kL 

-Id 

^Id 

Sta- 

tus 


^kl 

’kl 


( 1 , 3 ) ( 1 , 2 ) ( 2 , 5 ) ( 5 , 6 ) ( 4 , 6 ) 
0 12 12 12 0 


Y X 


8 


6 7 

4 


( 2 , 4 ) ( 2 , 3 ) ( 5 , 5 ) ( 5 , 4 : 
0 0 0 0 


3 

-5 


2 

-4 


5 

.2 


H 


0,6 


0.5 

9 

+ 


2.5 


2 . 72 15 0.6 1 


(k,l) 


_ld 

'kl 

-kl 

-Id 

Sta- 

tus 


( 1 , 3 ) ( 3 , 5 ) ( 2 , 5 ) ( 5 , 6 ) ( 4 , 6 ) 
6 6 6 12 0 


( 2 , 4 ) ( 2 , 3 ) ( 1 , 2 ) ( 5 , 4 | 


0 

5 

-5 


0 

5 

-3 


Y Y X 

1 6 8 


6 

2 

-3 

-2 

2 

H 


0 

5 

-2 


0.6 1 .0 

5 


S 

2.5 


(k,l) 


5 


75 10 1.0 


"kl 

'Id 

kl 

^kl 


( 1 , 3 ) ( 3 , 5 ) ( 2 , 4 ) ( 5 , 6 ) ( 4 , 6 ) 
6 6 1 11 1 


( 2 , 5 ) ( 2 , 3 ) ( 1 , 2 ) ( 5 , 4 ) 


-Id 

Sta ^ 

tus 


Id 

kl 


w. 


X X X Y X 


6 6 4 5 7 


* 


5 

3 

-2 

-3 

5 

S 

1.5 


0 

0 

2 


H 


0 

5 

-2 


6 
-1 
2 
1 

_1 

S S 

1 .0 2.5 







(M) 1(1,3) (3,5) (1,2) (5,6) (4,6)1 (2,5) (2,3) (2,4) (5,4) 


(It,!) (1,3) (3,5) (2,5) (5,6) (4,6) (l,2) (2,3) (2,4) (5,4) 


(k,l) (1,3) (3,5) (2,5) (5,4) (4,6) 1(1,2) (2,3) (2,4) (5,6) 













(3.116) 


MLnimize • £ e. x + T f v 

(i,d)eA ^3 i3 i/id ’ 


siibject to 


f 

= < 0 


^ if i=s, 


I I = ^ 0 , if i/8,t, ¥ie]J, (3.117) 

V , if i=i;. 


(d,i) €1(1)3^ (i,j) co(i)« I 


° -S'ld "ij* » (i.J) ' *. 


y c X < 2* 
(i,d) 'a - 


(3.118) 

(3.119) 

(3.120) 


This problem, vihich is referred to as ^ is clos^ related 
to the CSUCB problea. The only differoace between the two is '&at 
one of the constraints snd the objective functi<Ma are interchanged. 
The following theorem establishes a relationship between the optimum 
solutions of the two problons. 

Theorem 3.7 : Let x?. be the optimum solution of the OTOE problem 
13 

for D = D° < I) with Z° as the value of objective function. 

— max « 

Then, x?. is the optimum solution of H for Z = Z° with as 

the value of objective function. 


Proof : Suppose that the theorem is not true aid there exists a 

solution x^^ of PI for Z = Z° with the value of objective 

function D* < D°, Let 2 ^ = I c. . x*.. We note that 

0 ^ e A ^ 

Z = Z because x*. is also a feasible solution of the OTGE problem 

■ .r-' ' 

for D = D° and Z* < Z° contradicts -flie optimality of the solution z°^. 
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!>*■ < D in^lies that in the optimum solution of the OTC® 
problem for D - D , the constraint {3i5) is not a binding constraint 
and, hence, p = o. This contradicts oiir assunption that y > 0 for 

^ — ^max* Z* = Z , is -jiie optimum solution 

of PI for Z = Z°. 

It is now obvious that the ®CE algorithm solv^ PI * Ihe CKGB 
algorithm can also be viasfed as if it treats Z as a parameter md 
increases its value until either the prescribed value of Z is reached 
or infeasibility of PI is indicated^ 


3 *2 A Biciiteiia Network Problem 


Consider the following bicriteria networi: problaas 

C I c X 1 


lex + I f. y 

(i,5)eA (i,3)e A^3 


subject to 


f-v , if i=s, 


(3.121) 


I I 0 , if i^,t, ¥ieir, (3.122) 

(j,i)ei(i) 3 (i,j)eo(i) ^ ^ 


° -'‘a 


( 3 . 125 ) 


° -J'ij -hi- ‘ij’ ’’ (^>3^ 


(3.124) 


!nhe parametric programming tedkiniques have be«i used extoisively 
to solve multicriteria math^iatical programs, Bor bicriteiia mathe^ 
matical programs, Bacopoulus and Singer [4 J have suggested the 
•constraint czlterla ^proach*, l.e,^ maxi. mi, zing (or minimizing) one 
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cid.i;:eiia keeping the oliaer criteria ae^ a coastraiiit. Biey daowed 
how aXl the efficimt elutions c^ be gmerated by paraiictrically 
varying the level of the constraint# 

Ihe ONCE algorithm treats the criteria ( 1 e. . x. . 

T . X . (i,3)eA 

L asa constraint varies its level parametrically. 

(i,j)eA^J 

fflius the solutions of the CMCE prohlean, whicii generate the tradeoff 
curve, axe the efficimt solutions of (3.121) - (3.124) and the points, 
where slope of the ciirve changes, correi^ond to extr*^e efficient 
solutions, 

3, ? SPECIAL CASES 

In this section, we consider some special cases of the CSCE 
problem and point out the resulting siE¥>lifLcations in the CSCE 
algorithm, 

3.7.1 Constrained ffinimum Post Blow liublem 

Hie mathematical statonait of the Constrained IfiLnimum Ciost Blow 
(CMCE) problem is ae follows: 

Minimize I (3.125) 

{i,3) t A ^ 

subject to 

y X. J X. . 

(j,i)e i(x) (i,j)eo(i) ^ 


Cv , if i=s, 
i 0 , if i/Sftf ¥-ieir, 

I T , if i=t. 


(3.126) 
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0 1»„ , »(i,3) tJl, 


I e. . X. . < 

(ijj) eA 


(5.127) 

(5.128) 


In the GtfCE pitihlem, if we put f. .= 0 ^id a_ = b. . for 

15 ID 

each (ijj) ^ A, then it reduces to the (MGF problem, ibllovm^ 


sinipli£L Cations result in the OTGE alg02?ithEi it is applied to 


the CMGE problems Ihe sets of saturated basic arcs (y) ^d the iioiib^ic 
arcs at their upper bounds (v) are ^ooply^ Ihe optiiauii basis structure 

of the GMCF problem is (X U {(pj^q) } ^ere (p,q) is always ^ 

unsatuiuted basic arc. Only two nuEobers, c, ^ and d, need to be 
defined for each nontree-.arc (k,l). Hence 1he sets ^ md Q are 
empty , Further, the optimality conditions for 1±ie (HDF probloi are 


c + p d = 0, 
pq ^ pq ^ 

°kl (k,l)^iUU, 


(5-129) 

(5.150) 


On accomt of tiiese siii?3lifications, obtaining liie characteiistic 
inteiTral ^d performing the dual simples: iteration is easier. However, 
interpretation of the algorithm does not change. 


5.7.2 Bicriteiia Trapsportation Problem 


Hie Bicriteria Irsnsportarfon (BT) problem has been considered 
by Anej a and Hair [ 5 ] . suggested a method, i*ich obtains all 

the efficioit solutions of the BI problon by solrfng a sequence of 
transportation problems. In Ihe OJCE problem, if we put b^^.= ® 

and f. .= 0 for each (i,3) e A, ^ the underlying network is that of a 



transportation problem, then the CHCE algorilto obtains ail -fee 
efficiaat solutions of the BT problem, Tte CSCE algoritim, then, 
solves one tr^isportation problem in the beginning in each 
subsequoat iteration an effort equivalent to that of performing one 
simplex iteration is, required, Hoace our algoiiihm is coiii)utationally 
superior to the method of Meja ^d Hair [5 ] for solving the BI 
problem. 

3.7.3 Constrained Shortest Path Eroblem 

Ihe Constrained ^loitest Path (CSP) problem is to determine the 
shortest path from source to sink satis:^ing an additional linear 
constraint. Ihe CSP problem can be foimulated as tbe CMCP problem 
with additional integrality requirKneats on vaii^les. Eecently, Idae 
CSP problem is proved to be HP-.coB?)lete by dbowing that the in^sack 
problon is a special case of tbe CSP problem [4-9 ] , Since the CSP 
problem is HP-conplete, efficiaat algorithE® yi^ding near-optimum 
solutions are highly desirable. Ihe CMCE ^ori-aim incorporating the 
integrality modification suggested in Section 5.4,2 provides such a 
solution. 

3.7.4 Other Special Cases 

The problem to allocate a given budget D to increase the 
capacities of various arcs so that the cost of flow in the network is 
minimized, cen be formulated eis followss 
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Mardmize Z = ^ c. x 

(i,3) e A * 


subject to 


(5.131) 


I 2C. 

(j,i) el(i) (i 


I 

»j)«: 0(i) 


X. . 

13 


j-v , if i=s, 

j 0 , if i^,t^ ¥ieN, (3.132) 
1 V , if i=t. 


° -"'ij - \3-\3’ ''<"- 3 ) 

I f. . y. . < D. 

(i,j) s A - 


(3.133) 

(3.134) 

(3.135) 


¥e refer to the problem (3.13l)-(3.135) ae the GaPaeity 
Expansion to Minimize Goat of Slow (GBJGE) problem. In the ®CE 
problem, if we put e. = 0 for each (ijij) e A, it reduces to the 
CEMCP problem, !Hiere is no (diange in the OICE algorithm vidiea it is 
applied to the GMUP problem. 


Consider a variation of the C1M3E probl^ in which no fixed 
budget is prescribed and capacities of arcs can be increased as long 
as the total cost of flow and capacily expansion is ndnimized. !Ilie 
formulation of this problem caa. be given by (3.88)-(3.91 ). !Ehe solution 
provided by the GUCE algorithm, vdien the set of active ares is empty» 
is the optimum, solution of this problem. 

!l!he problem of allocating a prescribed budget D to increase 
edacities of various arcs so '^lat ■&e flow in the networtc is maximized, 
is a special case of the CMCP problem. 3MS problem can be solved 
by the CNCE algorithm but the problaa is too simple to be solved by the 
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CKCE algorithm# Instead^ we consider this problem as a special case 
of the parametric netvsiork feasibility problem considered in the next 
section* Similar commaats apply to the constrained ma Ti Tmim flow 
problem vdiich is a special case of the GMOP problaa# We ccxisider the 
constrained maximum flow probloa in Chapter IV* 

3, 8 FMAmiRlC NEFWOm CAPACITY EKPANSKM FBOBIEM 


In this section j we consider c^adLtated traassfaipment networks 
in which supplies at source nodes aad demands at s ink nodes are linear 
functions of a single parameteri aiwi the capacities of arcs csaa be 
increased by incurring additional cost* Hie problem considered is to 
determine the least cost of increaising arc capacities so that the 
changing demands are met by the changing suppliesi resultnng from the 


changes in X , We term this problem as the Parametric Hetworfc Cag)acil 


(HJCE) problQn* Hie mathematical stataaent of the HJCB 


problsan is as follows J 



0 < Xy < yy, » (i,j) C J, 

° -^13 -\ 3 - ” 13 ’ ’ ^ 


(3.156) 

(3.137) 

(3.138) 

(3.139) 


where 


G. . > 0 for each (if^) c 
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2h8 HTGE piDbl^ is « -t.* 

ai'i of too parametilo newioii feasibility 
Pfobl^ 0003110.01 l„ adopts, n tbo oopadty eap^alon p.obloa. 

0003110.01 in tblo Coptot. .e nao tbo oonoopts le^p^ m both of 

these chaptejrs to obtain an effloiait i. ■*,. * 

iicimt algoELtiffli to solve the PNCE 

pioblem for all v^uea n-p \ -p , 

vaLues of X for which a feasible flow exists. 

^timalil^ Condltinno 

let us partition A into the sets X,Y,1,D and T such that 

a spanning tree, She optimm basis structure of the HICE 

problem IS define! as (x u Y, *!(* Is similar to tho optimum 

basis structure of the MfflOE pioblem. The only difference between 

the two is that the basis of the HfOE piobla. consists of a spanning 

tree instead of a spanning tree aogmentol by a nontroa-arc for the 

MOUCE problem. Tbls difference arises because tho alditional linear 

constraint is absent in the mas problem. Wo now dorivo tho nooossaiy 

aud sui^ficiBiil/ condiljioiis ‘to b© saI:! Ktr -*-* i. 

sauxsxiea by Idi© op'fcxsiiii basis sijinclnxr© 

of the HTGE problem. ae dual of the HfOE piDblaii is 

Maximize I (r°+Xi^)n ^ b. .p. _ T b 

(3.t40) 


sub j ect to 

- ° ^ e A, 

*^ij“ *^i3 — ^ij* ^ e A» 

Pi^ > 0 aad 5^^ >0, ¥(i, 3 ) e i, 


(3.141) 

(3.142) 

(3.143) 
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where 6^^ are ihe dual Taxiables associated with the 

constraints (3.137), (3.138) and (3.139) respectiTely. 


33a e con^ilemaataiy slaclaaess conditions for the MCE problem ai« 



, ^ (i,3) e A, 

(3.144) 


, ¥ (i,d) e i. 

(3.145) 

X. .(p. =0 

13 13 1 3 

, ^ (i,3) e A, 

(5.146) 

y. .(c. .+ 6. ._p. .) = 0 

, y (i,3) e A. 

(3.147) 


She following theorem, whidi is similar to Iheorem 3.1, is used 
to deiive an alternate rqpresontation of "Sie conplaaaataiy slackness 
conditions. 


33aeorem 3.8 : Eveiy optimum aalution of "aie MCE problem satisfies 

Proof : Suppose y. . >0 and x. . < b. .+ y for an arc (i,j). 

It follows from (3.144 ) and (3.147) that 


. . < b. . 
3-3 





c . .+ 5 ,. p . . = 0. 
13 13 13 


(5.149) 

(3.150) 


Hence 

X. . < b. .-fy. .and y. . > 0 ==> c + 5. . = 0. (3.151) 

X 3 xj xj ij 

Since 5 >0, (3.151) implies that c. . < 0 vhich contracts our 

id — ’ . 13 — 

assumption that c. . > 0. 



Using aeorem 3.8, the conditims (3.144)-(3.14-7) can be shorn 
eg^uivaloit to the following conditions: 


0 . <b. . => 

13 13 

1 

0 

u 

1 

y^^=0, ¥(i,3) ® A, 

( 3 . 152 ) 

b. . < X. . < a. . => 
13 13 13 

TL-ir = c. .and 
3 1 ij 


(3.153) 

IT.-. TT. < 0 =*• 

3 1 

^ 3 “ 0 and y^.= 0 , ¥ (i, 3 ) e A, 

(3.154) 

0 < c => 

3 1 i3 

X = b. . aad 

13 13 

yy' 0, y (1,3) ^ A, 

( 3 . 555 ) 

TT .-IT. > C. . => 

3 1 13 

X. .= a. . and 

13 13 

=^13“ %3-'’i3’ 

(3.156) 


It is seoa from these conditions liiat in the optimum solution 


~ max, (o, J (ijj) ® A. 3his relationship 

permits us to consider y . inplicitly. It further follows from 

(3.152)-.(3.156) that the necessary and sufficient conditicms for a 

feasible basis structure of the HFCE problem to be s^i optimum basis 

structure are that there exist mrabers ir. satis^lng the following 

3 

conditions: 


Vi " ° ^ 

® X, 

(3.157) 

TT .-ir. = c. . ,¥ (i.j) 
3 1 13 * V 


(3.158) 

»^(ij3) 

el. 

( 3 . 159 ) 

O lt-^ -\-°i 3 

® U, 

(3.160) 

•tr.^TT. > c. . ,¥ ( 1 , 3 ) 
3 1 13 

® ¥. 

(3.161) 


We refer to the conditions (3*t57)-(3.16l) as the optimality 
conditions for the HfGE probloa. In the subsequent discussion in 
this section, the optimum basis struccure and the optimality conditions 



are used in referoice to the MCI problem. We refer to ir. satis:^ing 

3 

the optimality conditions as 1iie optiman dual varlsfeles g^sociated with 

the optimum basis structure (X y Y, L,C,t), It may be noted that given 

any feasible basis structure (X UX, 1,11,?), k. cm be taaiquely 

3 

determined by setting = 0 and -aioa using (5.157)-(3.158). If 

9 3 

thus obtained, satisfy (3.159)-(5.16l ), then tf. are the optimum dual 
variables and (X U Y,I,U,T) is an optimm basis structure. 

3.8.2 Development of the Algorithm 

In this section, we presQit in detail Idse development of an 

algorithm for the MCE problaa. It is assumed that a value of X , 

say X t is hnown fbr vdiich the MCE problem vdth y. .= 0 for each 
o j-3 

(i,3) e A has a feasible solution. Initially, the al^rithm obtains 

m optimum basis structure for X = X^. The value of X is then 

increased and an optimum basis structure is maintained at ea<ri step, 

The algorithm thus produces an optimum solu-cion for all feasible values 

of X >X , lo obtain ah optimum solution for all feasible values of 
— o 

X < X , r^ is replaced by -r* for eac^ i e I and the same algorithm 
— o^ i 1 

is applied, 

A baisic feasible solution for X = doioted by 

obtained by the method described in [l0 ] , In this basic feasible 
solution, let W^, and be the sets of basic arcs and nonbasic 

arcs at their lower and middle bounds respectively. It can be easily 
verified that the ba^s structure (W^ U ^ , Wg, W^, 4>) with the 
associated flow x°., together with l4ie numbers ir = 0 for each 



j e Uj salJisfies the optimalily conditions sndy beace| is en optiiauia 

basis structiire for X = i , 

o 


Let (X U y,L,U,?) be £0 optimum basis structure of the IfGE 

piobloa for ^ — 2^ wil3i the associated flow x. . . If X is increased, 

—13 

then, the flow on arcs changes due to changing supply and dem^id 
requiremmts, 31ie changed flow must satisfy the flow r^trictions 
(3»7)-(3.1l) of a feasible basis structtire. Since the numbers ir. 

J: 

are uniquely determined for a given basis structure, the optimality 
conditions remain unaffected by increase in X , Hence the hi^^t 
value of X, say X , upto which X can be increased without changing the 


optimum basis structure is determined from the feasibility considerations 
(3.7)-(3.1l). 


It follows from Section 2,5.1 that fbr X e (X. , X) » How in 
any arc (i,i) e X U Y is of the fom ^ + (x- jO vhereas flow 

on other arcs remains unchanged. It follows from Sieoran 2,6 that z. . is 

-m 

givm by 


z. . 

13 


= I 


kef. 




(3.162) 


where I is 1he set of nodes in the reafLting subtree containing 
3 

,node j when (i,3) is dropped from Ihe basis, 3h the n^lhod to 

compute z. . for the problQn desciabed in Section 2.5,1, if we put 

b* = 0 for each (k,l) e TI, it determines z. . for the MCE problem, 
ki ^3 

Having de1;erinined z. .! j is con^^uted £com liie following inequalities? 


o <x_ + (x- A) ^ (3.163) 

+ (^- £ \4 » ^ (i»3) ^ 

X3 ~3.3 13 


(3.164) 



®or each arc (i,4) e X define -Oie nuafcer M. . as fbllows: 

3-^ 


AX = -<|-X--/z.. 

IJ A -I ^ 




~13 ij 


5 if z. . < 0, 

13 


, if z. . = 0 . 


(3,165) 


Similarly, for each arc (1,3) e Y define the number AX. . 

^3 

follows: 


13 


13 -i3 i3 


13 “i3 i3 


if 

z. . > 0, 

13 ’ 

if 

e.j. 

A 

Q 

if 



(3.166) 


Using (3. 165) and (3.166), the conditions (3,163) aad (3.164) 
reduce to 


Hence 


(X_2)<AX ,¥(i,3)eXUY. 


X = X. + min, { AX . . } . 
(i,3)eXUY 


(3.167) 


(3.168) 


Ibr all X e (x, x), (X U Y, L,U, ?) is •fee optimum basis 
structure aid the optimum flow is givai by 


X. . = 4 

13 ^ 


X. . + (X« ^ z. ¥ (i,3) e x U Y, 
-a 3 137 ^ 


(3.169) 


X. . 

-^3 


y (i,3) X U Y. 


Let X. . be the optimum flow for X= X , airther increase 

13 __ 

in X is blocked by an arc (p,<l) £or whi<^ AX „ = min, AX. . } . 

(i,j)eXUY 



If (Pj*!) e thoa - 0 or and if (pjt) £ Y, thoi 

^pq “ ^pq ^q* ^ is incsre^ed furtber wilbout changing the 

opi/iisuiB. basis sl^ructia*©, flow iia arc "violafcs ite bomads. . 3SiiiS|: 

a dual sinplex iterafiou is perfomied to obtain an alt ©mate optimuiE 
basis structure #iich may allow further inerts e in X * 

The dual sxmplea: iteration is performed by dropping the arc (p><l) 
from the basus and selecting an arc in 1 {J ^ U ^ to ®iter liie basiS| 
vfaere L, U and V are the sets of arcs in I* U U (J V U UPta)J 
having x = 0^ b and a respeetiTely* Warn the arc (p|q) is 

-**J xj 

dropped from the basis, two siibtrees are fonaed. Let f f be 

p q 

the resulting subtrees containing node p ^d node q respectively. 

All the arcs which have their one end point in T and another in T 

P <1 

constitute a cocycle Q« Let us define the orLmtation of the coc^cle 
along (p.q) if z > 0, snd opposite to (p,q) if z ^ < 0. In fact, 
the orientation of the cocycle is defined along the direction of flow 
increase with X . Let Q aad j9 be 1±ie sets of arcs in the cocycle 
Q along and opposite to its orientation respectively, (Sbr the sake 
of convenioace, subscripts of Q aad are dropped). 

Since Idae new basis must be a sp waning tree, laae entering axe 
must belong to ^ U .2 • 1 = (Q fl i) U CS fl U) and 

E = ("^ fi ij) U CQ n Y)* If 1 UJE is eii9ty» then Q cV and 

Q cL. In other words, for each (i,j) e Q and x^- 0 

for each (±, 3 ) c Clearly, the flow om not inci*ease along the 

orientation of the cocycle end, hence, the I!fGE problsa is infeasible 



for X > X, 

il, _ for each 

IdL 

(k,l) e U U V 


However, if E is noneaigjlgr, "thoa define a m^toer 
(k,l) e E U U a niMjer for each 


as follows: 



, if (k,l) e L , 


(3.170) 


■^1 - TTjj. » if (kjl) e U, 




"l" 


'KL 



, if 


f if 


(k,l) ^ U , 

(k,l) £ V. 


(3.171) 


It follows from the optimality conditions that 

S Uu and > 0, ¥ (k,l) e U u ¥. (3.172) 

Let 


V = min. [ min. {p_ ) , min. {y.,.}]. (3.173) 

(k,l) e E ^ (k,l) e B ^ 

Let this minimum be achieved for an arc (g,h) ifci(di can belong 
either to 1 or E. Qlie fbllowing theorem provides an optimum basis 
structure for X = x, 

theorem 3.9 : Ihe basis stiucture obtained by r^lacing Ihe arc (p,<l) 

by the arc (g,h) as an unsaturated basic arc if (g,h) e E , or as a 
saturated basic arc if (g,h) e j, is an optimum baais structure for 
X = X . : 

Proof ; It is ea^ to prove using (3.157) end (3.158) that the opti mum 
dual variables, tt*. , with respect to the new basis structure are 



, ’f i E f , 
P 


’ = ^ IT . +U, ^ 3 e5? * if 


z >0, 
Pt ’ 


^ “ P, e if Zp^ <0. 


(3.174) 


It is clear from (3.174) 1iiat if for an arc (i|3) both i md 3 
belong either to or to the optimality conditions (3.159)- 
(3.I6I) ronain unaffected by ehaages in the optimum dual variables, 

!rhus we need to consider the arcs belonging to Q UjQ only. Elrst 
consider any arc (i:,l) e Q . if > 0, then n ^ and 

irj_ = TTi + p ; and if < 0, then - P and »£ = In 

either cage, substituting ir» in (3,159)-(3.l6l), mi thm using 
(3.170) and (3.171), we get 

Pm. i p » ^ % f (3.175) 

Pjj3^+p 10 and if (k,l) e U, (3.176) 

-kl'^ ’^-° * ^ ■ (5.177) 

Bi^e conditions are obviously satisfied in view of (3.172) 
and (3.173). Similarly, it am be shown that if (is:,!) e the 
optimality conditions are satisfied, 

3,8.3 Descro-ption of the ilgoiitbm 

A formal description of the algorithm for the E^CE problem is 
givoa below. 

Step 1 ; let (X (J Y,I>,11»7) be ^ opthotm basds structure of the BfGS 
problem for X= Sjet , ^ Z(^ = 0 mi = 0, 

?ieN. 60 to Stqp 2#, 



Ste£_2; Let 2^^ be the optiiaum flow for ^ , Goiapute 
Define the niMjers &X for each (i,j) e X U Y as 
follows: 


z. 

3-a 


(b. .-X. .)/z. . , if a >0. 


AX. . = <: 

n n T 


13 


_x. ./z. . 

-13 13 


, if z^. < 0, f (i,3) e X, 










, if z. .= 0, 

’ 13 * 


Let AX= min. { AX. . > ^ X = X+AX. a)rallXe(X,X)- 
(i,3)eXUY 

the optimum solution of the MCE pix^blem is giv^ bj 


^ (X- X ) ¥ (i,j) e X u Y, 


X. . = < 

13 ^ 


X. . 

13 


, ¥ (i,j) y X U Y, 


^13 = 


max. {0, x_^} , ¥ (i,i) e A, 


and the cost of capacity expansicm is 


% Z(X) - Z(X) + (X_ 2L) J- °ii\i * 

(i,3)eY 

If X = STOP, otherwise update X. = X and identify the arc (p,(i) e 
X (J Y for which and go to St^ 5, 


3 ! Drop the arc (p,q) from the basis. Let T end I be the 

p q 

resulting subtrees containing node p and node q r^pectively 

Arcs from 1 to 5 ani 3?^ to T constitute a cootcIo* 
p q q p 

Define the oiioitation of the cot^de along (p,q) if z > 0 
and opposite to (p^q) if ^ sad Q be the 

sets of arcs in the cotycle along and opposite ro its oriai- 
tation respectively. Sirther, let 1 = (Q f) L) U (Q fl U) 
and ^ = (Q n U) U (Q n V), Tidiere L, D and ¥ are the 
sets of arcs in 1 yU U¥ ^(pjq)) having x. . = 0, b. . ssad 
a., re^ectively. If 1 U J is empty, go to Step 4; 

•i-J 

otherwise define a number ¥ for eadi (Ie,!) ^ D and a 

Kill 

number for each (k,l) e JE eis follows: 

{tt^-iTi , if (k,l) e L , 


Pld 


=< 


, if (k,l) eu , 


iiki=<: 


V\ °kl» 


Cl"’\ “ °kL» ^ 


Let V = min, f min. { p.,} , min, { } j. Let (g,h) be 
(k,l)eE ^ (k,l) 

the arc vdiich achieves ibis minimum. If CSt^) ® % thoa (g>h) 
enters the bahis as ah unsatui*ated basic arc| otherwise it enters the 
basis as a satumted basic al^. Update the basis structuare. further, 
update the dual variables as 






= ■< 


1T.+ P, ? 3 e 3! . if Z, 


PI 


> 0 , 


v.-p, ¥ 3 e t , if Zp^ < 0 , 



HTCE problQE is infeasible for X > X- STOP. 
3.8.4 M^umerical ExampTfi 


We now solve a numerical example to iilnsfrate toe various steps 
of toe PKOE algorithm, toe nettork is shown in ELg. 3.4. toe nui!d}ers 
r. and are indicated above toe node i. toe nun^eze a. .« b. . 

^ ^ 13* 13 

and c are moationed over the arc (1,3), toe problem is to be 
solved for all feasible values of X ^ 0. 

Solution of toe first six iterations is givm in Table 3.5. 

In the table, the symbol + doiotes the arc mteiing toe basis and 4 
denotes the basic arc leaving the basis. Basis in various iterations 
are shovm in ilg. 3.5. In tois figure, satuarated basic arcs are drawn 
in double lines. 


3,8.5 Capacity Expansion to Maximize Blow 

toe problem of allocating a given bxaiget B to incr^se 
capacities of various arcs so that the flow from source to sink is 
maximized, can be formulate as follows: 

Maximize v, (3.178) 


subj ect to 


1 X _ I 

(3,i) e l(i) ^ (1,3) €0(i) 


X. . == 


fLv, if i=s , 

0, if ii4 ,t, ?i e 

If 1 =^* (3.179) 


° i ’^13- ’ 


I 

(i,3) e A 


13 x3 




(3.180) 

(3.181) 


toere c. . >-0 for each 
X3 


e A. 




Table 3 *3 1 SolutiMi of fhe ne’fcisoii flow problem. 


Ita- 

fa*- X a(X) 


Basic arcs 


Honbasic arcs 


n 0 


0 Sta- 
tus 
2 . . . 

13 


Sta^ 

tus 

2. G,5 0 z. . 


1.5 8 10.5 

X X X 

3 0 1 




4,5 8 11.5 

X X X 

3 0 1 


4.5 1.5 0.5 


(1,3) (2,3) (2,4) (5,4) ; 

8 0 

0 0 

U I. 

I. L 

«» ' 1 

_ 

1 

, .1 

0 0 



(1,3) (2,3) (5,6) (5,4) 

8 0 

9 0 

U L 

U L 

4 

(1,5) (2,3) (5,6) (5,4) 

8 0 

9 0 

U Jj 

U L 


0 


2 

4 




1 Z(X) 


Ba3ic arcs 


Uonbasic arcs 


4. 2 4 


(i,j) (1,2) (3,5) (2,5) (4,6) (5,4) 


8 12 


X X 


0 ^ 1 




(1,3) (3,5) (2,5) (4,6) (5,4) (1,2) (2,3) (5,6) (2,4) 


8 8 12 


U V 


(1,3) (3,5) (2,5) (5,6) (5,4) (l,2) (2,3) (4,6) (2,4) 


11 11 14 
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We refer to the problsa (3.178)_(3..iai) as the Capacity 
Expansion to Maximize mow (CfMP) problem. In this section, we show 
how the PICE algorithm can be used to solve the CMP problem. 

We assume that the value of D is such that the constraint 

(3 .181) is a binding constraint, i.e., in aey optimum solution it is 

Satisfied as an equality. If the network contains a directed path P 

from soizrce to sink for M^Mch I c. . is a finite ntmiber, then 

(i,3) t P 

no matter how lange D is, constraint (3.181) is always a binding 
constraint. 

Using the complementary slacciness coiklitions of the CEMP problem, 
it can be demonstrated that the optiman solution of the OMP problem 
satisfies y. .(b. .+ y. .-x. .) = 0 for eadi (i,3) e A. Using this result, 
the compleoientary slackness conditions of the CEMP problem cen be proved 
to be equivalait to the following conditions s 


0 <x <b . ==> a _a = 0 ^ ^ (i»3) ^ h 

J-iJ * J 

b. . < X. .=> 0 .- 0 . = yc. . affui y. - - x - b , V (i, 3 ) e A, 
i3 13 1 3 13 13 13 13 

0 ,..0, ^ 0 X. . 0 and 3^44 1^(^j3) ^ 

13 13 13 

0 <0— a. < y c. . =5> X. .= b. ^d y . 0, ¥ { 1 , 3 ) e A, 

13 *^13 13 13 13 


V >0 => V^'t = 


y (D - 


(i,3) ^ A 


y c. . y. .) = 0, 


13 "13 


(3.182) 

(3.183) 

(3.184) 

(3.185) 

(3.186) 
(3.187) 


#iere 0 and y > 0 are the dual variables associated with the 

i- 

constraints (3.179) and ( 3 .I 8 O) r^pectively. 
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Substituting iT^= -o^/y for each iel in (3.182)-(3.187), 

we obtain TJhe following equivalent conditions: 



^13=0 , ¥ (i.s) e A, 

(5.188) 



(3.189) 



(3.190) 

0 ir.-'n. < c. . =>x. .= b 

3 1 13 ij 

aod y^^=0, ¥ (i,3) e a. 

(3.191) 

V > 0 => -n-.^TT = l/y , 

t s ^ 


(3.192) 

^ - I c.. y..) = 

(i,3) ^ A 

0. 

(3.193) 


ffe note that the conditions (3,188)-(3.19l) are included in 

the complemoitaiy slackiess ccmditions of the JfCE problem. Ihe 

condition (3.192) can always be satisfied by setting w = o 

s 

y=lA^. Thus, (3.193) is the only aidition^ condition. 

In the HICE problem, let us set 


, ¥ i e U, 



(3.194) 

f-l , if 

i 

= s, 


V 0 » if 

i 

s,t, ¥ i e H, 

(3.195) 

L 1. « 

i 

= 


» , ^ (i,3) 

€ 

A 

(3.196) 


aid consider the optimum solution of the INCE problem when Z(X) = D, 
This solution satisfies (3.193) and haice is ai optimum solution of 
the CEME problem. 
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When the PNCE algorithm is applied to me CMP problem, thm 
at eacn iteration it augments the flow from source to sink on the 
path consisting of basic arcs aad capaciti^ of ail saturated basic 
arcs lying on this path are incre^ed by the aioount of augmented flow, 
!I!he algorithm terminates ■ssh^ either the given budget is ei^austed 
or a path consisting solely of saturated basic ares is discovered 
implying that infinite amount of flow can be at^ment&i on this path, 

Pulkerson [^39 ] solves the OMF problem by solving a seq,uence 
of minimum cost flow problems* At eada iteration, a l£j)elling routine 
is executed to increase the flow from source to sink gjKi when the flow 
ca» not increase, the value of is incremented by one imit. In 
our algorithm, flow is augmented on the loiown path from source to sink 
and in two subsequent iterations, the value of ir . caa increase by 
more than one imit^ Hmce our algorithm requires less it^ations 
in each iteration less confutations are performed. 

Hu [52 ] has also proposed a method to solve the CMP pmhlm.,. 

At each iteration, his method defines arc lengths ^d augmmts the 
flow on the shortest path from source to sinl:. His method is siiidJ,ar 
to our algorithm* The essential difference between ttie two is that 
shortest paths are inflicitly muskrat ed by our algorithm, whereas 
these are obtained by solving shortest path pioblens in Hu*s method. 
Hence our algorithm is computationally superior to Hu*s method* 



CHAPIEB lY 


CONVEX COST NETWOBK FLOW AND mAfffi PHOBLiMS 
4,1 INTRODUCTION 

In this chapter, we stt% ttie ai nimm cost flow pioblem when cost 
of flow over each arc is given by a piecewise linear convex fimction. 

We refer to this problem as the _Goavex (bst Network Blow (OCNP) problem. 

Several practical problens can be foiranlated ^ -^e OCNP problem. 

Cost s-cructure of water transportation in pipeliae netwoiis is convex 
[ 33 ] • fransportation problem, with vsicesTtain dosands a penalty 
for short supply, can be fomulated as the C®E problem [21, 34 ] . 

Capacity expansion of networks [64 3 > end a problem associated with 
transmission of power in electrical netwoitos [51 ] are also CCNP problems. 

The COTE problem has beoi sttMied by several r^eaixhran. The 
following three approaches have be«i r^rted in the literstore to 
handle such problems: 

(i) Primal-Dual Approach i Hiis approach assimes cost fiaictions to 
be piecewise linear convex functions and is essaitially a generalization 
of the out-of-kilter method. An algorithm based on '^s approach is 
suggested by Minty [70,71] . Lawler [64] also discusses this approach. 

(ii) Shortest Path Approach : This approach consists of determining, 

at each iteration, a shortest path from source to sink based on incremental 
costs of arcs and augmentii^ flow over it, Hu’s algorithm [51 ] is based 
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on this approach. She algorLthm s^suaes that aOl cost ftmctions 
lineailzed along scents of ^ common lengih d, such -fliat all 
capacities are multiples of d, 

(iii) ff-^ative Qyole Appioadi ; In this ^pioaoh, the objective flmction 
is improved in successive feasible solutions by gugmenting flows in 
n^ative Qrcles based on incremental costs, lh«i no su(^ impzovemait 
is possible, the solution is optimum. Several authors have used this 
approach and the essential diffei*Qace is in the manner in which negative 
cycles are determined, Bor piecewise linear convex functions, Blonan 
and Bobeits [36] and Billet et, al [35] have offered node l^elling 
methods, and Beale [ll ] and Klein [59] have su^ested matrix manipulation 
methods to identi:^ negative (grcles, Bor continuous convex AmcticKis, 
enumeration is sxjggested by Moaon [6? ] . A moi^ efficient method is 
proposed by Weintraub [93 ] which identifies negative cycl^ by solving a 
finite ntanber of assignment probloiffl. 

It is wellKnowa that the GCNB problem can be traasformed to the 
miniiaum cost flow problem by introducing one varidJle for eas* linear 
segment, GMs traasformation enlaiges tiie network si*st antially. We 

note that in the optimum solution all the additional variable oen. be 
considered implicitly, dhis observation allows us to define the concept 
of optimum basis structure for the CGNB probl^, The optimum basis 
structure is then used to parametrize v, the flow to be transshipped 
from source to sink. Initially, the algorithm obtains an optimum basis 
structure for v=0. The value of v is then increased and an optimum 
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basis structxire is loaintained at each st^ mtil eL-^ier the desired 
flow is established or infeasibilily of C®P problem is discovered. 
Intui'Cively speakirig| the algorithm inplicitly ezxumerates shortest 
paths from source to sink and augments flow over these paths. Hence 
the algorithm belongs to the shortest path approach. Ihe conputatioaatl 
complexity of the algorithm is shorn to be oCmv). 

‘Ihe time-cost tradeoff problem in a CPM network is of consider^le 
practical significance. Cto accoint of its iaportancey the problem has 
been studied by .several researchers [l2, 33, 51, 40, 46 , 58, 63 , 78, 

80, 84 ] . We show that the CCM'F algorithm c^ be used to solve this 
problem. 

We also establish a relationship between the GCSP problem and the 
maximum flow problem with piecewise linear concave gain functions. 

This relationship is a generalization of tiie relationehip betwem the 
minimum cost flow problem and the maximum flow problaa with gains 
recognized by several resea3rchers [66, 74, 90 ] . We suggest modifications 
in the CdTP algorithm in order to solve this problem. 

The problem of optimally sd-locating a given bt»aget to increase -aie 
Capacities of various arcs to maximize the flow in a network is considered 
by sev^al researchers for different capacity ©mansion costs [ 9 , 39, 

50, 52, 79 ] . We show that whoi capacity eapansxon costs are given by 
piecewise linear convex fUactioi:^, this problem can be solved by the 
CCN'F algorithm. 
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A computer program was writtai for the CCIP algoiildan. Bie 
computational peifoimance of the program is reported. fhe pixjgram 
solves GCNP problems with 200 nodes^ 2000 arcs said cost fiaictic^is 
consisting of 5 linear segments in 57 seconds, 

4, 2 PROBLEM STATEMENT 


The mathoEatical statem^t of the C®P problem is as follows: 


Minimize Z = T 0. .(x. .) . 


(4.1) 


sub 3 ect to 


jjecu t:o p 

-V, if i=s, 

(3,i) e I(i) (i,d) e0(i) 1 

’ (4.2) 


Xf. >.0, ¥(i,3) eA, 


(4.5) 


vrfaere 0. .(x. .) is a piecewise linear convex function for eadi 
Ij' 3.3'' 

(i,3) e A, 

It may be noted that capaoily restrictions are not explicitly 
stated because th^jr can be incorporated in C (x .). 

J. J X J 


Let C. .(x..) consist of p(i>3) linear 
3.3 i 3 


. Let 


p = max. {p(i,;}) } . Associated witt each arc 

(i»3) eA 1 ^p(i,d) = t,p(i>3)+1_ 

(1,3) e A are nianbers 0 = ^ \ j < ... < b^^ b_ -...-b^^ 

and slope of (x^. ) in -fee interval 

asstjme that c^^ < c^^ < ... < vhich implies that ^ij(*xj) 

a convex function. lUrther, it is assumed that C^j(o) = 0 for each 

(ifd) e A. 



4,Z OPTIMALITY CONDITIONS 


Let us partition A into the sets 
such that B U B U U B^ is a spanning tree. If we set 

’ (i»3) e ? k = (4.4) 

then there exists a unique flow x. which satisfies the flow 
conseryation equations (4 #2)^ If this flow fuither satisfies 

1 Ic It 

^ij e B , ¥h=1,,..,p, (4.5) 

i O n 

thai we refer to B U B U ...UB*^ a3 a feasible basis aid 
(B^ U B^ U ... U B^,0 °,u\...,dP) as a feasible b^is stiuctiire of 
the CCaSTB problem. We refer to aa arc (ijj) e A as a k-basie arc if 
(i,j) e B^ and a Imonbasic arc if (i,^) e U^, We define m optimm 
basis structure as a feasible basis structure for whitdi the associated 
flow is the optimm soluticai of "fee C5CHB problem. In this section) 

we derive the necessary and suffLcient cojaditLons for a feasible basis 
structure to be an optimum basis structure. 

It is wellhnown [72] that iSie C5CKP problem can be trsnsfoimed to 

a linear program by introducing the variables ^ij*****^'^ 

' B ' Ic 

each (ijj) e A and replacing each by ^ y^. iHie COTP problem 

thoi becomes 

Minimize Z - J ^ y^^i » (4-6) 

(i,i) e A fc=1 



(4.8) 


° "y’- » (i.3) 'i. » k = i,...,P. 

It is also Mown [72] that eny optimm dilution of (4,6)-(4,8) 
satisfies the following conditions: 


y— >■ 0 ==>yA = » b^“^ ¥ & - 1 it-1 


k , k , Ml ^ 5 

y. <b - b, . i=»yr. = 0 
ij ij ij ij 


, ¥ £ = fc4'1,.,.,p. 


(4.9) 

(4.10) 


03ae dual of (4.6)-(4,8) is 


Maximize v(ir^_.iT ) - 
t s' 


I ! 


«?, , (4.11) 


(1,3) j A W ' « « ' W 


subject to 


W ij -°i3» ^(M) e A, k=1,...,p. 


^3 > Of ^ (M) s A, ?h=1, ...,p. 


( 4 . 12 ) 

(4.13) 


and the coii?)leiaentajy slacMess cwMjditicais are 

ic liC! "1 

° <^1.1 - “il ==> w °ii’ ’ ' *> ’ -f. (♦• 


ir.— ir < c. . 

a i 10 


:> yj = 0, ¥ (i, 3 ) e A , ¥fc=1,...,p, (4.15) 


”3-\ ’ “y V ‘i- ‘■5’ ’ ’ 


(4.16) 


Slhese conditions , in iriLsw of (4#9) (4 #10)^ ogn be writtoa. as 


0 <y^ < — b^^=> IT.— ^ (ifj) ^ ^ TsF^fmmmfpf (4*17) 

•'xj ij ij 3 ^ 2.3 

^ - —^1 -=>yL=*>4 4-' ^ =i,2,...,fc. 


07 < c. . 

10 0 1 13 


'ij 13 13 


y£^= 0 , ¥ £ = ^ (ifj) ^ 4, lc= 0 ,...,p. ( 4 . 18 ) 


(Mere = - « 


and c?t^ = » fte ®sB3h (i,i) e A. 

13 , '■ 



are 


The condi-cions (4.17) and (4.I8), in terns of 2. ’ 
lc»»l ic Ic k * 

^ (i»3) e A, ¥k=t,,..,p^ (4.19) 

k ^ kt-1 k 

°ij ^ ^(i»3)e A, ? fc=0,..,,p, (4.2o) 

It follow from (4.19) and (4.20) that the necessazy ^id 
sufficioit conditions for a feasible basis structure (b^ (j U ••• 

UB j U jU ) to be an optimum basis structure are '^at there 

exists numbers ir satisfying Ihe following conditions: 

tl 

’j-’i "" °i3 * ^ ® ^ fc=1»*..,P, (4.21) 

k k+l k 

1 ®ij » ^ (i»j) e IT, ¥ k?c,1,...,p. (4.22) 

We refer to these conditions as the optimality conditions for the 

CCUP problem. We refer to ir satisiying the optimality conditi<»is 

3 

as the optimum dual Tariables associated wi-fc ttie optimum basis 
structure (b"' U B^ U... jU"*,...,!]^). 

It Hay be noted that givaa any feasible basis structure 

(B^ U B^U., JJB^,TJ°,U^,...,D*’), it. cm be lailquely computed by 

3 

setting ir =0 and then using (4.21). If ir. thus obtained satisiy 

S 3 

(4.22), then ir. sire the optimum dual variahles aid (B UBU... U B^, 
3 

U°,Tj\,.,,U^) is the optimum basis strwture. 

4,4 DEVELOBdENT OF THE AWOmTEM 

IhB gjgA rn ±v>tn for 1dti6 CCJIF pioblQQi 1?r0ats as a paraiastor attd 
maintains ^ optimum basis structure at aYeay step. ImtiaUyt an 



optim-um basis structure for v=0 is obtained. me agoiithm thm 
determiiies sudi that this basis structure remains optimm for al 
Y £ A dual sin^lex iteration is then performed to obtain agi 

alt ei:nate optimum basis structu3?e for f #iich allow funiier 
increase in y# 33ae algorithm proceeds in this manner laitil either the 
desired flow is established or infeasibiiity of the problem is indicate* 

4»4*1 Initial Optimum Basis Structure 

We now maJse one eissun^itiwi wiiieh, not veiy reetiictiTe, 

simplifies the con^jutations ojnsider^ly# 

Assumption 4»1 : Hie network does not contain any negative cycle wife 
cl . as fee length of each arc (i, j) e A. 

iO “O 

This asstimption implies that x. .= 0 for eafe (if 3) e A, is fee 
optimum flow for y=0* Let f be the tree of ^ortest paths rooted 
at source with cl. as the length of each arc (if 3) ® A. Stnrther, 

let TT. be the laigth of the short^t path from aiurce to any node 

*3', ' 

jew. Th^ we know [22 3 that 

v.-ir. = c!. , ¥ (i,3) e 2!, (^*23) 

a 1 la 

TT-Vi , ? (i,3) elf (^-24) 

feere T = A^T. ■ 

In view of ( 4 . 23 ) snd ( 4 . 24 ), it follows feat (B^ U B^U... U 
■D°,U^,...fB^) is ah optisam baisis structure for v=Of where B = T, 

U°= T and B^= ... = B^= U^= .*. = • 



4*4.2 Gfaaragteiistic IntervyCL 


Let (B U B U.,, Ub , U be aa optimim basis 

struc Cure for v = t wi1ii the associated flow ^ . . Let Z be the 

-ij - 

cost of this flow, Purther, let be the optimim dual vaiiahl^ 
associated with this basis stmeture, we now deteimdne the moviimm. 


value of V, say v, sudi that the currait basis structure remains 
optimum for aay v e (v,v). lae in-serval (^jt) is Isnom as the 
characteri stip. interval associated with (B^ U B^U ...Ub^,D°,u\...,D^), 


In 1}he basiSy there is a ntiiqno path source to sink* le 

refer to this path as the basic pal^i # Let 1 be tile basic path in 

and P and ^ be the sets of forward and backward arcs in P r^pectiireljr 

Ihe Value of v car be increased hy sending additional flow froii source 

to sink through the basic patii# We assuiiet without any loss of generalityr 

that ^ = 0. using (4#2l) it can be shorn tiiat 

s 


\ = L k\-s* ^ ^ k°ii- 

^ k=1 (i,3)ePnB ^ k=1 (i,3) e P HB ^ 


- ! 


(4.25) 


!Ihus- ir. doaotos the oast of sendir® unit additional flow from 
' t ■ , 

source to sink through the basic path. 


Since it, are miqu^ determined for a given bgeis structure, 

D ■ 

the optimality conditions ranain unaffected by increase in v. However, 
flow on arcs belonging to P ch^ages aid the changed flow must 
satisfy the bound restrictions (4.5) of the basic arcs. We now 
deteuuine the TTi gyi-rmim additional flow, denoted by w, that can be sent 
without violating the bound restrictions. If we define a nuBft»er 
for eadi (i,i) e p as 



w. . = 


» if (ia) e P 


theaa 


Hmce 


1 x...b^\ 

! -13 ±3 * 


if (i,3) e P n B^, 


w = min. { w, } 
(i,3) e P 


(4.26) 


( 4 . 27 ) 


^ = • (4.28) 

Ibr all V e (t,t), (b’ U B^U...UBPy,u\...,lf) is the 

OpfilQUIB. b ssis S*tXliC*tm'6 1/b6 np' fci Trann fiovff is giv6ii by 
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+ 0 W , ?(i,3) e 1 , 
%j. - « w , ¥ (i, 3 ) e p , 
T (1,3) B , 


and the cost of this flow is 


2 = Z + 6 w ■ff . 


where 


e = (v - v) / (t - v). 


(4.29) 


(4.30) 


Lei; X. . doaote the optisajm flow for t = ir. further incre^e 
in T is blocked by 681 arc (a,g) e P for which p = w. If 
there are nore thm one arcs satisi^lng this criteria, 11181 let (a,B) 
be the arc which is nearest to the source. Let (a,B) e B®, If 
(a, 3) e then 5^^ = b®g } and if (o,e) e P, then . 

If V is increased fuiiher wdlhout changing Ihe basis structure, flow 
in arc (a, B) violates its r^pective bound, Hius, a dual sin^lex 
iteration is performed to obteiii^ an alternate optiiaUiii basis structure 
for V = T, which mgy allow ferther increase in t. 



Dual Simplsx Itei-a-ci^n 


The dual sionplex iteration is perfoimed by dropping the arc 
(a,e) from the basis and selecting aa arc belonging to if y U ... 
UU {(a,3)} to enter the basis. It be noted that (a, 6) can 
also ^ter the basis, but at a different status, i.e„ it would either 
be a (e+l)-basic aro or a (e_l)-basic arc. When (a, 3) is dropped 
from the basis, two subtrees are formed, let 3? gpd 1. be the 

S 

resulting subtrees containing source aid sink respectiv^y. Arcs 
from Tg to T^ and to constitute a cutset Q, let Q 
and _Q be the sets of forward and backward arcs in Q respectively. 
Since the new basis must be a spaining tree, the entering are must 
belong to Q. let us define the following sets for convenience: 


= <(i,0) e u° U U^U...Ulf U {(a,3)}: x. .= b^. } , 


ij ij 


¥ 1p=0, 




B _ {(0,6)} , if k=e, 


, ¥ k=t,...,Pa 


(4.51) 

(4.32) 


B 


, if ka^e 


We follow tlie convention that if x. .= b^., liien (i,^) e , iet us 
define the numbers f and for each (h,£) e 0^, ¥ k?=0,1, ...,p 


as follows: 



= IT W. 4- 

h - ‘^hfe * 

(4.33) 

4ji 

II 

1 

I 

(4.34) 


It follows from the optimality conditions that 



Let 


^ 4£-°» E Q* 


(4.35) 


V = min. 


[ min, 
(h,£)e 


T,'V- 


1^'^' 

(h*£)eQ ^ 


(4.36) 


If V = ” , then Q c T? and ^ c D° ^ !IMe implies Tiiat all 
the foOTraid arcs in Q have and all the backward arcs in 

Q have 0. Clearly, flow in the network equgas tbe TnaTiTwim 

flow with as the edacity of each arc (i, 3 ) e a end liie CGEP 

problem is infeasible for v > v. However, if p ^ •» , thai let 1±ie 
mini mu m in (4.36) be achieved at an arc (<lfr) idiich cgc belong either 
to Q or _Q, Let (q., r) e U , !Ehe following theorem provides ^ 
alternate optimum basis structure for v = v. 


Iheorem 4,1 ; The basis structure obtained ty replaclr^ the arc (a,&) 
6y (q»r) as a (f+l)-basic arc if (q,r) e ^ , or ^ a f^basic arc 
if (q,r) e is sn optimum basis structure for v = v. 


Broof : It cai be easily shorn using the optimality condition (4.21 ) 
that the dual variables, ''r! , wi-Si r^pect to the new basis structure 
are 


ir! 

d 



¥a eT 




(4.37) 


Ibr any arc (i,^) e A, if both i and j belcaag either to 
T or to T , the optimsiLity condition (4.22) is obviously satisfied. 
Thus, we need to consider arcs belo ngin g to Q only. Consider any 
arc (h,£) e Clearly, \ and = ir£+ p. Substituting these 

values in (4.22) and th«i using (4.35) £«d (4.34), we get 



(4.38) 

(4.39) 


which are clearly saT:isfied in 7iew of (4.35) and (4.36). Similarly, 

it can be shown that if (h, l) e ^ me optimality conditions are 
satisfied^ 


4«4-*4 Physical lutearpretation 

Initially I the GOTP algoiithia augm^ts tiie flow on the shortest 
path from source to sink# Sie flow blocking arc is dropped frois the 
basis and a cutset results* in arc i^aidi belongs to this cutset aP ^ 
yields inxninium increase in k is sheeted to eater Vtie basis 

tj 

thereby giving a new path from sovirce to sink on itiidi flow o&a. be 
augmented. Since is the cost of tmit alimentation frfwn source to 

sink, the new path is essentially a slwrtest path v?ith respect to the 
changed flow. Therefore, the al®3rLthm in^jlicitly aiumerates shortest 
paths from source to sink and augments flow over these paths. 

Hu*s algorithm [51 ] for liie GCKF problai define appropriate arc 
lengths at each iteration, solves a ^K>rtest path problem aod then 
augmoats unit amount of flow on l^s pa'fc. The proposed algorithm is 
superior to Hu’s algorithm in the following irespects: (i) flow augmenting 
Paths are obtained without solving i^iortest path problems and th^?^y 

Saving computations; (ii) the augmented flow at each iteration can be 

■ . k 

more liisn 1 #iich is a desir^le feature if toe veiues of b are large; 
and (iii) the algorithm can solve problems in toich b^. aire n o n i n t^er 


or evaa irrational 



4.B DESCRIPTION OF THE AISORITm 


A foimal staliQiiQili of fiis CCSUS’ algorithni is givsn bcOLow* 

Stsp 1 * IiSf T be the tree of shortest paths rooted at source with 
°i3 ^ iQigth of each arc ,(i,3) e a. Set B^= T, 

U°= A-T and B^= ... = B^= U^= *.. = if= 4> , (Ehe initial 
optimum basis structure is (B^ U B^ U . . . U B^,U°,U^ , . . .,0^) 

with the associated flow x. .= 0, V(i,j ) e A, Cbn^jute 

■ ^'3 , , 

IT by setting ir = 0 and using Tr.-ir. = c1., ? {ijj) > B 

3 . . 3 ^ 

Set T = 0, Z = 0. Go to Step 2. 


Step 2 : let P end P be the sets of fbrwaid and backward arcs 
in the basic path P respectively. Define the mmiber 
for eadi (i»3) e ^ as foUovras 


w. . = <, 

3-3 


bjj- ^ij » e 1 n B^, 

X..- b^\ if (i,j) e P flB^. 
13 3-3 V _ 


Compute 


w = min. ■f w ^ ^ 
' (i,3) ^ P ^ 


and 


y = min. ( w, (v-v) } . 


Update z. V and Z aS foUovra: 

■ ,^3 .. 

Zi^ + y , V (i,3) ^ P, 

z. . t ^ (i»4) ^ 

13 



V = T + y , 

2 = Z + IT y^ 

If V = T, go to Step 5j olSierwise go St^ 3. 

Let (a, 6) be HJae arc nearest to source ajnong the arcs for 
which = w. Drop the arc (a, g) from the basis. Let 
be the resulting stibtre^ containing source aod 
sink respectively. Let Q be the cutset separating 5^ 

Let 5 _2 be the sets of forward ar>d 

backward arcs in Q r^pectively. Blither, let 

= {(i,j) e if U U... U if U {(o,$)} IX., = bj, } , 

¥ k=0,1,,,,,p. Define the nuiAer for each e Q 

as follows: 


%1 * if (h,A) z’Q nD^ 


^hl 


= <! 


£ 

k 


’fh "■ °h£ * ^ n 


I?^. 


Gon^iute y = min, i V|, }* y = ® , go to Step 4} 

(h,£) e Q ^ 


otherwise identic m arc (q,r) fbr which = P . 

Let (q»r) ^ Obtain the new optimm basis structure 

by replacing (a, 6) by (q,r) as a (f+l)-baaic arc if 
(q,r) ® Q , or as a f-baelc arc if (qjr) ^ _Q, I^ate 


as follows: 



¥ J e fg , 

¥ 3 e f ^ , 


and go to Step 2. 



GGNi' problaa is infe^ible, SOSDP, 

Sfc^ ; The optimum flow of the CdTP pioblem is x. . with the cost 

^ : 

of flow 1 , STOP, 

4.6 CCMPLEXITT OF THE ALGORTTm 


The following theorem ^tablish^ a tight bound on the nuaher of 
coiQputations performed by the algorithm vdaen b^ are int^er, 

k ■ 

Theorem 4,2 : If all b^^ arc integer, the coii?)le3d.ty of the COTP 
algorithm is o(innv). 


lioof t let us caH an iteration of the COrP algorithm as a xmndegenerate 
iteration if positive ajnomt of flow is sMgmait^ on tie basic path in 
that iteration} otherwise it is called a degenerate iteraticai . let 
further define a nond^aaerate subtree as a sihtree containing source 
in which a positive anoint of flow csn be augmented from source to sJay 
of its nodes through the arc belonging to the aibtree. The algoiilhm 
follows the rule that the arc nearest to the source among the flow 
blocking arcs leaves the basis. As a result, T is always a non- 
d^esQ©ra1/© .subtree* Bach dual simplex iteration adds one arc to ihe 
noudegmerate'.siibtreet H^ce, after atmost (n— l) d^enerate iteratious^ 


the nond^eaerate subtree incLudos the sink and there is a nond^oiearate 

iteration. If all are integer. at leant a «at a«»t of 

flow is augmented in a nondegenerate iteratiim. Kie algoii14m, thus, 
terminates in atmost (n-l)T iterations. Since the nrmher of crai^aita- 
tione perfomed at each iteration are o(m), the complad-ty of the 
algorithm is o(mnv). 



. Bg5gLk _ 4.1 : It is obvious from the above pit>of tiiat the algorithm 
terminates xn a finite nuiiher of iterations even if b^. are not 
integers, 

4,7 NUMERICAL EXMiPLE 


In this section, we solve a nm^rleal exacple to illustrate various 
steps of the GCasri* algorithm. Hie network is showa in Sg, 4,1, Data 
of the problem for p=2 is given in fable 4.1, Bbdee 1 aid 6 are the 
source and sink nod^ respectively. ae flow to be ti^isshipped from 
source to sink is 15 units. 


The steps of Ihe algorithm are sMnmaiized in fable 4,2» Hie 

desired flow is establi^ed in six iterations. In -fee tajile, 4 
indicates the arc leaving the ba^e 4 indicates "fee arc entering 


the basis. Hie optimum flow is given in the last row of the table. 

Hie basis obtained at each iteration is shown in KLg, 4,2, The optimum 
dual variables associated wife fee nodes are mentioned above thorn. 







Table 4,1 J Data of the numerical exaai^ile 
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Table 4.2; Soluticaa of tee numerical exanple 


Basic Arcs 


lonb^ic Arcs 


0 0 


(ifj) (1,2) (1,3) (2,4) (3,5) (4,6) 

0 0 0 0 0 

Sta^ B^ B^ B^ 

tus 


(2,3) (2,5) (3,4) (5,4) (5,6) 
0 0 0 0 0 

if if 0° . u° if 







(1,2) (1,3) (5,6) (5,5) (4,6) 

' (2,3) (2,5) (3,4) (5,4) (2,4) 

4 0 

0 

0 

4 

0 

0 

0 

0 4 

b”* 

e' 

B^ 

B^ 

n° 

u 


0 1 

B U 

- 6 

10 

3 


— 



. »» 
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66 

Ite-T ^ ^ ~ — — — 

V ■ 7 "D - ■ ■ 

tion Basic Arcs Honbasic Arcs 

Ho . 




1 (i>3) 

(1,2) (2,5) (5,6) (2,4) (3,5) 
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! X. . 

0 

0 

6 

6 

6 

0 

6 

0 

0 6 
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b' 
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1 
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(i, j ) 
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(2,3) (1,2) (5,4) (5,4) (4,6) 
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13 
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2 

8 

6 

6 

0 

8 

0 
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b' 

B^ 
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B 

if 


if 

O ' ^ 1 ■ 

U U', 
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13 

4 

- 

2 

- 
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- 

- 

- 

■- 



(i,3) 

■ (1,3) (2,5) (5,6) (2,4) (5,5) 

(2,3) (1,2) (3,4) (5,4) (4,6) 

15 

65 

X. . 

7 

2 

9 

6 

7 

0 

8 

0 

0 6 
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4, 8 CONSTRAINED MAXIMIM FLOW FROBIEM 


The problem to maximize flow in a netTOrk subject to a budgetaay 
constraint, can be stated as 


Maximize v, 
subject to 

-T, if i=B, 

' , ^ , ^-ii - ' I X = I 0, if i^,t, ¥i e S, 

(j,i)ei(i) (i,j)eo(i) ^3 ^ ’ 


(4.40) 


V, if i=t, 


"-"ij-V* 

(i,i) ^ A 


(4.41) 

(4.42) 

(4.43) 


We refer to the problem (4.40) - (4.43) as "Sie OMastrained Maxinaaa 
ELow (CME) problem. In this sojtion, we diow how the <30* problem 
can be solved by the COTE algorithm. 

The complementary slackn^s conditions of the CUE problem can be 
shown equivalent to ■^e following emditions: 

(4.44) 

(4.45) 
^^.46) 

(4.47) 


“ Vj "'■v ' *> 


y > 0 => Og-Oj. - 1* 

y (D - I C X ) = 0, 

(i,j) " A 


(4,48) 


where o. and y are -toe dusdt^defcles s^sociated with the constraints 
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Substituting TT. --a^/ y for eadi j e i m (4.44)-(4.48) we 
get the following equivalent oonditicxisi 


0 < X. . 

la 

< b. . = 

=> 

W.-TT. = C. . 

D 1 13 

f ^ 

(if j) 

e A, 

(4.49) 

IT.— TT. < 
1 1 

c. . = 

=> 

X. . = 0 

13 

f ? 

(i»j) 

e A, 

( 4 . 50 ) 

A 

1 

0. . =: 
IJ 

=> 

X. . = b. . 

13 13 

, ^ 

(i,i) 

s A, 

( 4 . 51 ) 

V 

> 0 =: 

=> 

TT -.IT = l/p 

t s 

f 



(4.52) 

p(D _ 

( 

I 

i»i )4 

( 

A 

2. . X. .) = Ow 




(4.53) 


Consider the CCSJF problaa with n=r aid b? . = b and c"* = c 

3-a 3-J ±3 

tor each (ijj) ^ A. Ihe conditions (4.49)— (4.51) are the conplementaiy 
slacloiess conditions of the 0®! problon. (lie comiition (4.52) cm 
d-ways be satisfied by suitably choosing P . Hence the only additional 
condition is (4.53). ' 

Consider the optimusi solution of the COTF problem whoa Z = B. 

This solution satisfies (4.43) and hence is an optimna aalution of the 
CMP problQn, However, if the maxiaam flow is established in the 
network before Z equals B, tboi the constraint (4.43) is iM)t a bindiaag 
constraint and this flow satisfies (4.53) wi1ii V = 0. 

4,9 TIME-CX)ST TRADEOFF PROBIEM IN CPU NEWORKS 

Consider a project consisting of a nianber of activities. She 
precedoace relations between these activities ar© u^icated by idaitiQring 
the activities wilh the arcs of a directed network S = (H,a). I»et 
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aode^aAd t mark the initial and final events of the project 
respectively. 


Each activity (i,j) e 4 i^as ah associated normal comnletion 

time and a crash completion time, !Ehe cost of performing the 

activity betwe® these two extreme times is given by a piecewise linear 

convex function vihich consists of q. . linear segments. Associated 

^ q. +1 

with each activity (i,j) e a are numbers ^0, 

and the cost of short eoing the coipletion tine of '&kB acti¥ity C^ij) 
in the interval by one vmit is b^^ sudi '^at 


13 


0 b^. < b^ <.,.<b^^ 

ij 


The problem of finding tbe mlniTruim cost of rfinrt e rring -ttie project 
to a given duration ^ is knovai as the lime-Oost Tradeoff (fOf) problaa, 
Ihe plot of the least project cost as a function of "file project duration 
is Icnown as‘ the pro.1 ect cost curve . In this section, we show that the 
CGUE algorithm can be used to obtain the project cost curve, 

.Associate a variabl® with each i ^ H, vhich denote the time 
at #iich event i starts. Bet t. , demote the shortaxLBg of activity 
(i,j) in the interval (a^T, a^. Since the cost ftaiction is convex, 
shortening of tbe activity in differait intervals can be treated 


ind^endently. The TOT problem can be stated as follows: 

X ^ 

Minimize L ^ “-li H » 

(i,j) "A 1^1 

subj ect to 

- Pg < X , 


(4.54) 


( 4 . 55 ) 



156 


) D + 

i 10, f (i,3) e A, 

° ^ Ci»3) e Aj ?fc=1,,.,. 




(4.56) 


(4.57) 


Associate the variables y, aid 5 with the coustraints 


(4.55), (4.56) and (4.57) respectively. Ibe dual of the 
is as follows: 


fCT problem 


MaJcunize 


subject to 




y a* X T T ^ 

'i,j) e A (i,3)eA k=t 1 


k 


Xv, 


(4.58) 


pv, if i=s. 


^ ^4,-- I X. . =< 0 , if i/s,t, ¥ieS, 

(j,i)ei(i) (i,j)eo(i) ^ » 


V , if i=t, 


^ij" % -\i* ?k=1,...,q^^ , 


(4.59) 

(4 .60) 


^lOj ¥(i,i)e^ andS^^>0, ¥(i,j)el, ¥k=1,...,q^^ 

(4.61) 

1 2 

Erom a, . > a. . >••• > a. . a^d (4*6o), it follows that optiiioiii 

ij ij ij 

solution of (4.58)-.(4#6l) satisfies *aie followiiag mndition: 


fif. = max* {0| X. b. . ¥ (i,^) e A| ip=1» •••> <L.: • (4 *62} 

Xj 13 xj XJ 

We can substitute (4*62) in (4*58) eliiiinate (4*6o) from iike 
constraints^ We then get the following equival^t problaaj 
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Vi 


Maximize Jalx^ 7 r‘'xk k^-i, t 

^ i '%i”Q--i-i ) aax. {o,x. ._b. .>_xv, 

(i, 3 )eA 1 e =1 ^3 * X 3 ’ 

(4,63) 


subject to 


I x.._ T X = <! 

(d,i)'i(i)^^ (i,3)€o(i)^3 


f~ 

-r t If i=B, 

0 , if l/s,t, ¥iei, (4.64) 

V f if i=t, 


?(i,j)eA, 


(4.65) 


Let 

Vj ^ Vj ^ ’ (i»3) ^ A, (4.66) 

°ij “V* ^^^*3) ® A, ¥k= 1 ,...,p^^ , ( 4 , 67 ) 

b^^. =0 and t)iJ^ = (X), ¥ (i,j) e A. ( 4 . 68 ) 

iUrtfaer^ add an arc (t|S) to the netiiork with imbouoided capacsity 
and cost of flow X ♦ We then get the following eqiaiiral^t flow pioblai; 


Maximize - Y(X) = ^ I CL (^ ) ^ X x. , 

(1.3) cA^5 « 

or, equivaloatly. 

Minimize Y(x) = I 0..(x ) + X x , 

(1.3) £l ^ “ 

subject to 

I X. . _ I X. . = 0, ¥ i e H, 

(j,i)ei(i) 3 f (i,j)eo(i) 

3 ^. ^0,^^(M) e A, 


( 4 . 69 ) 


(4.70) 

(4.71) 

(4.72) 


vdiere C (x ) is the pieeessise linear convex function as defined 
ij il^ 

in Section 4,2. We refer to 1iie problem (4.70)-(4.72) as the Mpdifi^ 
Tim&.Gost a?radeoff (MTCT) problem. 
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Remark 4.2 : The MTCT ^ ■ 

an IS not the dual of the TCf pioblems 
however, both have the same onti™™ 

soiu-cion. !i3ie optimum value 

of tiii0 MTG!r pioblsm iq +ri "hia, t - 

be multiplied by _t to obtain the optimum 

value of the TGT pmblem. 

a.e p„3ect cost ^ ^ 

IS varied from the normal duration of the piojeot to «ie crash 

^113? sfi/XOH O ^ iiii © ■n to i 

project. We now how the aCMF algorithm obtains 
the project cost curve. 

let us partition A into tbs seta b’.b^ 3^,0°, n’ 

and define a feasible basis structure of the IUCT ptoblem as 
(B UB U...UB U{(t,B)>, *ere B’uB^U...UB’’U{(t,s) } 

is a spannine tree, and b'^ and / satisfSr the flow restrictions 
(4.5) and (4.4) respectiveUr. It may be notei that the basis of the 
MTOT problm consists of two subtrees, rtich contain source ^ sink 
respectively, joined by the arc (t,a). It caD be e^iiy shorn that 
the necessary and sufficient conditions for a feasible basis structure 
of -the MTCT problon to be an optinami basis structure are that there 
exist numbers satisfying the following conditions: 


13 


^ “ ^t ’'s » 

(4.73) 

^ ~ ^ (ijj) ^ ® j ¥k=1,,.,,p, 

(4.74) 

- \ » V (ifj) e U^, ¥ fc=0,...,p. 

(4.75) 


We refer to -Mae conditions (4.73 )-(4. 75) as the optimality 
conditions for the MTCT problaa. Observe that the conditions (4.74) 
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and (4.75) are the optimality CMiditions of the OCMF pioblem. We refer 

to as the optimum dual variables of the lEDCf problem. It follows 

from Remark 4.2 that if are the optimum even start times of the 

TCT problan, then p. = _ u. for all j e I. Without any loss of 

generality, we assme in the subsequent discussion that ir ss o* 

s 

We first note that any GJM nstwork is acyclic h en ce 

Assumption 4*1 is satisfied* We also note that for eaA 

ij 

(ijj) ^ -^9 which implies that mspcirngg flow in the network is unbounded* 


(Zi ■^) he one of the characteristic internals of the GCSiF 

problQDa and the associated optinaam basis structui*e be (B^ U 

,U^) with the optimum dual varidales Let 

3 

(b’^U B^U U°, be the alternate optinaaa basis 

structure for v = v with liie associated dual varig^Jles Oj* 
follows £rom (4.37) that 


a. 

3 


» Tf 3 e Tg, 

TT +]i, ¥ j e . 


(4.76) 


Let J = «TT^ ^d 1 = -Oj. . Clearly ^ . 

Let B^ and be defined by (4.32) and (4.3l) respectively. 

It is easy to see that (B^ U B U,,,UB^ U {(t,s)}, U , U ,...,ir) 
is the optimum basis structure of the MIECT problem for X — X with 
the associated optitmam dual varieties It is also easy "to see that 

the same basis structure is the optimum basis structure of KPCT problem 
for X = \_ with the assooLated dual variables Oy 


It is now obvious 



140 


that (S' u S^u...Ub'' U {(t,e)} , {'.....JP) 

structure of the MIOT problem for all X e ( X , I ) 
optimum dual variables ir* = 60 .+ (i.e) ir., ^ere 
Hence, the optimum completion time for eadi (i,j) e 
(_X, I) is (Tr| _ irp and the slope of -T(x) is 


is the optimtna basis 
with the asajciated 

e= (S:_x)/(3Ld. 

A in the intervd. 



How, consider the optimum basis structua« of the CC31P problem for 
v — V, Ihe amount of flow which can be augmentei on the basic path, 
w. Can be finite or infinite. If w is infinite, then it means that 


all the arcs in the basic path are forward arcs ai^ their completion 


times are reduced to the crash completion times. Haice, the proiect 
duration is reduced to the crash duration ^ the fCf probloi is 
infeasible for X < ^ , However, if w is finite, '&ea dropping 

the flow blocking arc yields sa optimm basis structure of the MIGI 
problem for X= ^ ^ which may allow fiirther decrease in X , Thus, we 
see that performing the dual simplex iteration for the GdTP problem 
corresponds to finding the characteristic interval fbr the liEPCT problem 
aid, conversely, finding the characteristic interval for the G(SfP 
problem is equivalent to performing the simples: iteration for the MTCS? 
problsn. 

We observed that for all X e ( ^ , X), the optimum completion 
time of activity (i, 3 ) is given by - €) ihere w* = 9 ( 1 _ 6 ) 

and e= X )• Substitution from (4.76) yields the following 

expression for (it* - irO ; 
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’( - e w, ¥ (i,3) e q ^ 

(Tr.-,rp +0y, v(i,j) . 

( V’^j) » ’ (i,3) i Q. 


(4.77) 


AS X varies from ^ to ^ , 8 varies from 0 to 1 . Accxntiuigly, 
the optimvim completion times of all activitieB belonging to Q decrease 
by an amount 6 v | these tim^ for all activities belonging to 
increase by an amount 6 v • ?^ere^ the optimum conpletion times fbr 
other activities remain unchanged. The algoiitimi, liaus, idaitifies a 
cutset at each iteration and completion tim^ of all the activities 
belonging to this cutset are modified so that the project duration is 
decreased using minimum additional cost. 


Haillips and Dessouky’s algorithm [78] for the TCT probloa defines 
appropriate lower and i;^per bounds for arcs at each iteraticai, identifies 
a cutset by applying a cut search procedure aid modifies the consJletion 
times of all activities belonging to this cutset, ilie proposed algoiirhiH 
is superior to the Biillips ^ D^sou^’s algorithm in the following 
respects; (i) cutsets are implicitly enumemted 1die algorithm egad 
thereby enhacing its efficiency! aid (ii) more general cost functions 
are considered. 


itilkerson [40 ] aid Kelley [58] have proposed labelling procedures 
for the TCT problan which handle piecewise linear convex fmctions by 
introducing one arc for e^h linear segment, Ihe algorithm proposed 
by us avoids this enlarging of "Hie network whicli is highly desirable 
from conputational point of view# 
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4, 10 CONCAVE GAIN MAXIWM PLOf PWBim 


input flow, the maxinam flow with gains probloa 
detez^es the output flow ^en arcs have positive linear 

|^4T 3 « In sPA+*ii^y«4 

on, we consider a gena:^ization of this 

problem where arcs can have oiecewic^ n 

P ecev9xse linear concave gain functions. 

Mathematically, this pinblem can be stated as 


Maximize g ^ 
sub j ect to 


^3^0, V(i,3)sA, 


( 4 . 78 ) 


a» if i=B, 

0 , if ij4,t, ? ieH, 
e , if i=t. 


(4.79: 

( 4 . 80 ) 


where is a piecewise linear concave function consisting 

oi P(i,i) linear segments. We refer to this problem as me 
Concave Q-ain Maximum Jlow (CGMF) problem. 

Let p = max. {p(i,o)} , AsK)ciated with each arc 
(i»3) eA 

(i»j) ^ A are nxmhers 0 = b?. < b^ = hP(i»5)'^'l_ _ 

i3 ij i3 i 

and slope of in the interval is r^^ sudi 

that r! > rj > ... > r?. >0. 

. . T ,1 , ■ T 1. 


ij' ij' 

ij ' "ij ^ ‘ ^ ij 

Let us introduce the variables yL| y.. y: . for eadi 

ij ij ij 

Y • V ^ fe 

Vi,j) e A and replace each x by I 7. . in «ie OSSF probloa. 
We then get the following problem: 
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Maximize $ , 
subject to 


(4.81) 


X ? k k . I ^ F“» 


£ 0(i) k=1 


l^g, if i=t, 

(4.82) 


Jj: t ki-1 

°- ^ij -\j " \ j (ij3) e A, ¥|j=1,...,p. ( 

1 

It is ea^ to see that (4.8l)-(4.83) is equiralent to tke 
OGMB' problem if its eveiy optimum solution satisfi^ the following 
conditions; 


(4.83) 


f = D. , • . 

13 i3 ’ 


id id id 


^ ^ ~ 1j...( k- 1 f 

(4.84) 

¥ £= fc+1,.,,,p. 

(4.85) 


We now show that these conditions are indeed satisfied. Ihe 
dual of (4.81)-(4.83) is 


Minimize -op„ + 

3 


subject to 


(i,3) e A k?=1 


^j ^d ~ *^i “^ij ^ c A, ? k=1,...,p, 

- '^t - 


® xj i Of ^ (if J ) ^ ^ fc=1f Pf 


(4.86) 


(4.87) 

(4.88) 

(4.89) 


where p. and 6 ^. are tbe dual vaiiahles s©a>ciated with (4.82) 
and (4.83) respectively. 
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The complemoitaiy slackness conditions of (4*8l)-(4.83) axe 


f- k /. k , k ^ . 

'A. 

(4.90) 


(4.91) 

6 (- P^-1) = 0* 

(4.92) 


which can be equivalently stated as follows: 


k . k , k-.1 
0 < y. , < b. b_ 

IJ 3l3 

^ Pi~ ^ (1 »d) ^ Aj ? k=1^...]ipi 

(4.93) 

k 
r. . 
ID 

> P . 

^ ^ij” ^ (^»d) e Aj ¥ k=1y ...])p} 

(4.94) 

k 

"'iD 

< Pi 

k .k -k-l „ V , 

^^13 ^ ^3" ^3 » ^ ^ k=1,...,P, 

(4.95) 


6 > 0 

=:^P s? _1, 

t 

( 4 . 96 ) 


Theorem 4«3 : Every optimum solution of (4»8l)-{4.83) satisfies 
(4.84) and (4.85). 

Proof : Suppose that the theorem is not true and there exists an arc 
(i, 3 ) for which y^^ > 0 and y^^ < - b^^ such that i < k. 

Since all ^ G» it follows from (4.93) aid (4.96) that <0 

for all 

and 

Using (4.87) and (4.98) we get 


ID 


j E W« It flirther follows from (4.91 ) ^sd (4.90) that 




k . jf k 

y- > ° t'fy - "i * ‘y 


(4.97) 


±i ‘’ij 13 


(4.98) 
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i3 ~ ij to . 




Hoace 


(4.99) 


since c^ < 0 and > 0 , (4.100) iiipUco that r‘. < r' 


(4.100) 


fc 

ij 


which contradicts our assun^jtion that r^ > r^ 

ij ij* 

aorefore, the problem (4.81 ) - (4.83) is equivalent to the 
oam problem. This equivalence suggests that the OSMP problem can 
be solved as the laaximum now problem with linear gains by introducing 
one arc for each linear segment. This transformation alaiges the 
network substantially. We now suggest modifications in the GGSF 
algorithm in order to solve the CSGMF problem Yhich does not require 
any additional arc to be introduced. le first en^^loy logaillifflic 
transformation to siiiQ)liiy (4.93)-(4.96), 


let 


\ = log (-p^) , Vi e F, 


= log CVr^^), ? (i,j) e A, ¥ k = 1,...,p. 


C4» ’KJl) 


Note that <e^^ <...<e^^ for all (i,3) e 1 . aabstituting 
(4.101) in (4.93)-(4.96), W6 g0t the follovsing cjoii^liioiis; 

° ‘ ‘ ‘ir “S’ ’j-'i ' “S' ’’ 


y^^ * ^ (i>l) 4, ¥ k = 1y ...|)P| 


(4.102) 

(4.103) 


t -ir < 0^ ^ 

1 i °il 

==» ’(M) ®A, ¥k?^,...,p, (4.104) 

e > 0 ==> ir = 0. (4.105) 


r i 
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Ihe ooMltious (4.102)-(4.104) are the eonplementaiy elaotoeee 
conditions of the Oats piobia,. ae additionel conditions (4.105) 
can aiwsys be satisfied, as we show later. Iherefore, we find that 
the complemoataiy slactoiegs conditions of the C®P pioblem and the 
CG-MF problem are same. This is a strong relation^iip bet»em ■Qie two 
problems. 

Assumption 4.1 in the context of the CXJMF problem lilies ttiat 
the network does not contain any flow generating directed <^cle with 

4a ^ ^ ^ ^ ^ ^3rfLnold [47 ] 

has shown that any optimum basis of Ihe problem is a spasming tree. 

Hence, the concepts of feaisible ^d optimum basis structures for Ihe 

1 

CCaTF problem ronain valid for the OGMB' problem. Uhls observation 
permits us to develop the OGMF algorithm, similar to ■file C®F algorithm. 

The OGMP algorithm treats o as a parfijaeter. Bae initial 
optimum basis structure of the OGMF problem is obtained by solving a 
shortest path problem. The associated optimum dual vaxidal^ are 
obtained by setting 1'’,^. = 0 and thga tising (4.21 ). Smiing the 
characteristic interval for "this basis structure is slightly involved 
due to the gain functions. The duelling procedure described by 
Grinold [47 ] to augment the now from source to sink can be modified 
for our problem. Calculations in the dual sin^lex iteratitai remain 
unchanged, except that the dual variables are updated by substracting V 
from TT^ for all 3 «= r . 
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Intuitively speaking, the GGMF alsorLthm 

^ aLgorithffi augmaats the flow 

fiiom source to sink on the oathn • 

paths With maamm gains itech are iaplicitly 

enumerated by the algorithm. The . 

nnitaiess of the CGMP algoritim 

follows from Remark 4*1* 

4 . 21 CONVEX COST CmCITX EXPmsiON MAXmm FLON PumiEM 

problem of optimally allocating a prescribed bixlget D to 
increase the capacities of various arcs, in order to maximize the flow 
in a network, is of considerable practical significance. m Section 3.8, 
problem with linear capacily expaision costs was considered. In 
this section, we show that when capacity expansion costs are given by 
piecewise linear convex functions, this problem can be solved by the 

CCSF algorithm. The mathematical formulation of this problem is as 
follows: 


Maximize v, 


sub;j ect to 


X.. - 


(d,i)ei(i) (i,i)eo(i)^^ 


(4.106) 

jCv, if i=s, 

, 0, if ii4,t, ?i£ir, (4.107) 


*13 


(1,3)L“« 


^0 and y^^ > 0, f (1,3) e A, 


(4.108) 

(4.109) 


(4.110) 


where D^j(x^P is a piecewise linear convex functitaa wi-ai positive 
slopes. We assume that ® 0 for each (i,3) e A. We i^fer 
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to the problem (4,106)-(4^110) as the ^nveg Post flanyitv amnnainn 
Maximum ELow (CCOIME') problem* 

We assume that the constraint (4*109) is a binding constraint. 

If this constraint is not binding, th€n the GGC3£F problem is equivalfint 
to the maximum flow problem where y. . are at their i^per peimissible 
bomds# 


I*t follows from the above assumption and (4e108) that the opti»iis 
solution of the CCCMP problem satisfies the folio wing cooditioass 


y^. = {0, 1, ¥ (i,j) e A. 


(4.111) 


We Can svabstitute (4,1 1l) in (4.109) mA ^iminate y^^. We 


then get the following eq.uivalsit problem: 


Maximize v, 
subject to 

I 


(4.112) 


X. 


(j,i) e I(i) (i,j)eo(i) 


-v, if i=s, 

I = <; 0, if i5^s,t, ?ieH, 

T, if i=t, 


I 0 (X )<D, 
(i,3) A A 


(4.115) 

(4.114) 

(4.115) 


Ay i. 0, ¥ (i,j) e A, 

where 0. .(x. ) is the function obtained by translating 
i3 

by b. . units. • 

3-0 

The following theorem establish^ a dose rdationship betweai 
the optimum solutions of the GCCMF problem and the GOTF problem with 
C. .(x ) as the cost function for arc (i,j) e A. 

10 iO 
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for T-v rtth 2° 

y lo an optlMM solution of the OCOW proU„ fOn d . f ^th 

T as the TOlue Of objectiTefuaotiso. 

o * *^0 the optitaan solution of the CCCHiP problem 

for B = z“ ,4th the value of objective function ,». let 

^ \ ^ ^11 loaow that ^ • 

(ijj) c ^ u l m i mffii cost of 

Shipping v° anomt of flow fnom aourca t» sink. We also inow 

eveiy additional flow ffoffl aource to sink requiree additional cost 

because nonn^atiye slopes for each (i,j) e A. Hms 

> V implies > Z® which contradicts (4,109). Hmce v»<v° 

^ ^ 2 , Since is also a feasible solution of the CGCMP 

piDblem for D = Z°, we must have v» = v°. Surther, Z* = Z° 
because is also a feasible solution of the CCSfP piobloa fcr 

r = v and < Z° contradicts the optimaUt 7 of x° , Iherefore, 

^ o # 

V _ v* and Z - Z and x^^ is sn optimum solution of the CCCMF 
pimblem for D> Z°, 


This theorem shows that the GCSTP alsorithm also solv^ Idie 
CCCIM]? problem, Ihe CCKP algoilthm then can be viewed as treating 
h as a parameter. Xt starts with — 0 £d3d increases its value 
subsequent iterations until Ktther the prescribed value of D is 
attained or the maacimum flow is est^li^^ed in Ihe network. 
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4.12 COMPUTATIONAL BESUms 

A computer pTOgram, named PCCHF, was written in 33HPEAE-I? for 
the CCNF algorithm* The PCCMF stores the network in the form of m 
arc list and computations are performed by using a minor variait of the 
augmoated predecessor index method described in Murly ^72 J . Bie 
listing of the program is given in the Appendix* The progrsBi was 

debugged and tested on DBG-i1090 computer ^stem. The progra® occupies 

I5n + 5m + 2irp words of central memory. 

We solved a nunher of network pioblans by PGOIP for varioxjs 

values of m,n and p and noted down the computational times, le 

also observed the effect on computational times of variations in •tihe 

ic Ic 

process by which b. . and c. . were generated. In this secticai , 
the results of these experiments are presented. 

The computational experiments reguire a large ntm^er of randomly 
generated networks. Initially, we tried to use the program NETGM- a 
wellknown network gomerator [96 ] . However, our code of ITETG-ffi 
could not generate correct networks aid we were laidile to find the 
bugs. Hence we wrote our ovai network goierator NS havii^ similar 
characteristics as those of HETGM, The program NG first generates a 
skeleton network and thoi adds arcs to it until specifications are met. 
In this manner, quite well structured networks are gaierated. 

The data for the eircs was gaierated in the following meaner: 
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‘’y' ”'>^3 (10,30), Vk=1,...,p, 

cl = IRAJJ (10, 20), 

Ic 

°±3 " (2»10)* Vk = 2,...,p, 

where IRAK (a^B) is a fuucliioti 1«hat geaerates a l ytl •foTTnly dis'tributed 
integer number in the interval (a,B). 

Network problems were g^erated by 2JG for various values of "fee 
parameters; n, m and p. !I3ie parameters n^ m and p were respectively 
varied from 20 to 200^ 30 to 2000 and 1 to 5 » S)r each set of parameters 
two problems were solved and average executicoi time to establish 
flow at minimum cost was noted (the execution tiB» is the CPU time 
exclusive of input and output times)* iSieBe times are pr^oated in 
Sable 4.3. 

Ihe computational times demoi^trate the efficiency of the €K2fF 
algorithm for solving large-sized problems. Problems with n = 200* 
m = 2000 and p = 5 are solved by PC®P in less than 37 s^jonds, 
v7e observe that if any two parameteio are fixed, the computati-ohal 
times vary almost linearly with the third pai^eter. In fact, the 
approximate computational time 3? can be calculated by using '&e 
following formula: 

_5- . ■■ ■ 

T = 2nn4> wio seconds, 

■ J |i ^ 

We also varied the manner in which b^^ ar^ were generated 

said observed its effect on the oonputaticnal times. Some of these 



152 


Table 4.3 : Gon^jutational times of KXMF (ac«nition 
times in seconds in Dl(l_1090) 


lo* of 

No . of 



No. of segments 



nodes 

SlTCS 

1 

2 

3 

4 

5 


30 

0.01 

0.02 

0.03 

0.05 

0,06 


50 • 

0.02 

0.05 

0.06 

0.08 

0.12 

20 

100 

0.05 

0.08 

0.11 

0.16 

0,22 


150 

0.06 

0.11 

0.16 

0.20 

0.34 


200 

O.C^ 

0.15 

0.23 

0.33 

0.45 


100 

0.07 

0.16 

0.25 

0.34 

0.42 


200 

0.22 

0.45 

0.70 

0.85 

1,^ 

50 

500 

0.53 

0.91 

1.49 

1.80 

2.55 


800 

0,81 

1.49 

2.41 

3.14 

4.08 


1000 

1.13 

2.15 

3.35 

4.42 

5.83 


200 

0.28 

0.59 

0.87 

1.22 

1.39 


500 

0.95 

1.77 

2.58 

4.22 

4.71 

100 

1000 

1.75 

5.28 

. 6.18 

8,07 

9.79 


1500 

3.56 

5.86 

8.07 

10.23 

14,18 


2000 

4.36 

8.51 

13.49 

17. (X) 

22.03 


500 

1.84 

3.66 

5.70 

8.01 

5.25 

1000 

4.43 

8.73 

9.95 

17.27 

18.89 

2000 

8.21 

14.64 

20.56 

27.16 

36.75 


200 
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observations are given below* 

(i) \hen c. . are generated by using + Ibm (2,X), ^ 

is increased from io to 100, the computation time 
increases by 2C^. 

1*en b^^ are generated by using b^^ = b^^ + im (l0, X), 

the computational times decrease as X is increased. However, 

this change is not as pronounced as in the case of 


(ii) 
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COHaODINS RajAHKS 

In this dissertation, we studied the role of paraiaetric progiransffing 
in solving several goieralizations of some of "tiie weHknowa network problems 
We first studied the parametric analysis of some network problems and 
then considered the following gaier^zations; 

(i) the parametric network feasibility problem; 

(ii) the constrained network capacity expansion problea; 

(iii) the parametric network capaci-ty expansion problem; 

(iv) the convex cost network flow problem; 

(v) the time-cost tradeoff problem in CiM network; 

(vi) the concave gain maximum flow problem; ^d 

(vii) the convex cost capacity expatt&Lon maximum flow problem, 

A coufflion methodology based on parametric prograjmaiE^ is surest ed 
to solve these problems, Ihis methodology utilizes the concept of 
optimum basis structure, Ohe basis s'cructure is a geneiralization of 
the working basis for bounded vaiiable linear progrEffos. Optiiaality is 
maintained because the optimum solution satisfies certain properties 
TO3ich are not satisfied by an aibitsrary feasible solution. Bus is 
particularly true for convex functions vdiere these properties permit us to 
consider several variables implicitly. Since in most of the problaas 
considered by us, convex functicais ^pear either in the constraints or in 
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the objective functioni maintaining the optimality is very us^hl. 
ihepeforej parametric progra mning is used to parsmetzize one parameter 
in order to find the optimum solution for ■fee desired value of the 
parameter. 

As indicated above, the optimum basis striKJture considers seve 3 ?al 
variables implicitly^ Ihis, in certain problems, r^ults in iim>licit 
consideration of several constraints. Haace the size of the optimum 
basis is reduced considerably* On accoiaat of this proper^, basis 
triangularity is preserved in most of the problems. This, in turn, 
yields algorithms which are intuitively ^pealing slid Bass' to isplenmt. 
Integral conputations is also a bypiodvtct of 'Hiis proper^* 

Another advantage of the parametric approach is -aje availabilii^ 
of additional information which may be of ocaisiderable OiSiageiial use. 
Algorithms based on this approach yield the optinaa solution fbr all 
values of the pareuaeter contained in ^ interval. In planning new 
product introduction or exploring expaasion possibilities, this infoxxaaticai 
will hi^p in decision naMng, Also, "Hie paregaetric ^proach easily 
Imds itself to the sensitivity analysis, 

SIsere are several natural generalisations of the problems considered 
by us* A linear programming problem with convex cost fisactlons is one 
such generalisation. The oonc^ts developed fbr the cmvex cost network 
flow piobloa not be directly ^pUcable to this problem. Some 
difficulties are oacountered in sheeting lii© parameter to be parsmetrised 
and -Hnfiine; the initial optimum basis structure. However, the oonc^t of 
■Hie feasible basis structure will stiH be meful. It may be possible 
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to develop a goaerjilization of the siiuplex method vshich moves frcaa 
one feasible basis structure to ttie siaother until m optimum basis 
structure is obtained. 

In most of the problems considered by usj convex functions appear 
in the constraints or in the objective fmction, nevertheless^ our 
approach can be extended to solve problems viiere convex problems ^pear 
in the constraints as well as the objective functi^, linear pxogxssmning 
extensions of these problems can also be handled. Hie following allied 
generalizations can also be solved by suitable modifications and 
generalizations of the algorithms proposed in this dissertation: 

(i) Parametric analysis of the maximum flow problas with gains, "Bie 
maximum capacity path problem end the Tniniimm ratio network 
problems. 

(ii) The parametric network feasibility problem with nonlinear supply, 
demand and capacity functions, 

(iii) The miniTmitn cost flow problem with two eaiditional linear 
constraints, 

(iv) The parametric network capacity expansicaa problem with convex 
expansion costs. 
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